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HYDROLOGIC STUDIES BY ELECTRONIC COMPUTERS IN TVA 


Willard M. Snyder,! M. ASCE 


SYNOPSIS 
Application of electronic computers in TVA watershed studies is pre- 
sented. A comprehensive multiple regression analysis program for the 
IBM 704 computer is described. The use of this program in developing 
and studying hydrologic prediction equations is illustrated by discussion 


of numeric results of fitting two equations to data. 


INTRODUCTION 


Hydrologic studies by the TVA are carried on under a program of research 
and development of watersheds, which are tributary areas of the Tennessee 
River. The program is administered by an Advisory Committee on Tributary 
Watersheds. The purpose of the program and the role of hydrologic investi- 
gations are contained in the first two statements of Understanding and Beliefs, 
as formulated by the Advisory Committee. These are as follows: 


1. The welfare of the people is the clear goal of the tributary watershed 
program, and their full participation is essential to its success. 

2. Research activity on hydrologic withrelated problems and development- 
al watershed activities are broadly distinguishable. Both may be carried on in 
the same watershed area, but each is designed to achieve different but related 
ends. Research activities are designed to produce reliable facts from experi- 
ments with single or defined combinations of controlled practices. These facts 
- Note.—Discussion open until July 1, 1960. To extend the closing date one month, a 
written request must be filed with the Executive Secretary, ASCE. This paper is parte 
ofthe copyrighted Journal of the Hydraulics Division, Proceedings of the American So- 
ciety of Civil Engineers, Vol. 86, No. HY 2, February, 1960. 
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1 Head, Hydrology Sect., Hydr. Data Branch, TVA, Knoxville, Tenn. 
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are basic to the guidance of action programs in developmental watershed 
projects, where practices are not controlled. 


TVA has a variety of programs for development of the resources of the 
watersheds of the Tennessee Valley. These include cooperative research 
projects with state universities, test demonstration farms, fertilizer research, 
and reforestation activities. However, these programs, when undertaken sepa- 
rately, do not achieve maximum development of resources. If the welfare of 
the people be the primary concern, watershed research and development must 
be undertaken on an integrated basis. The problems for which the best possi- 
ble solutions are sought are found in the interaction of such factors as land- 
use, soils-water relationships, and the socialand economic levels of the areas. 

A watershed approach to an integrated solution is necessary if one wishes 
to completely evaluate a program of development and improvement. Starting 
with rainfall as one of the resources of an area, the changing hydrology must 
be considered in relation to the total development of resources of that area. 
Therefore, hydrologic studies and analyses are undertaken as one part ofa 
major program. Causes and effects of many varying factors must be evaluated 
as precisely as possible on watersheds designated for research activities. 
These results form a basis for development programs on other projects. 

Evaluation of hydrologic data in an integrated program requires development 
of ways and means for (1) the adjustment of data to eliminate variability in- 
herent in uncontrolled experiments, (2) detection of hydrologic changes, and 
(3) the synthesis of hydrologic response under changing watershed practices. 
To accomplish these itis necessary to formulate and evaluate effects of vary - 
ing amounts and patterns of rainfall, soil moisture levels and rates of loss, 
among others, and eliminate their effect on volume of runoff, peak rates of 
discharge, and rates of sediment movement so that the residual data, after this 
hydrologic adjustment, are less variable than before adjustment. These re- 
sidual data then form the basis for evaluation and synthesis of external effects, 
such as watershed management, vegetal cover, and possibly soils types. 

Much of the analytic effort involved in hydrologic studies is, of course, re- 
ducible to established computational forms. Where numeric procedures exist 
or canbe devised they can usually be programmed for electronic computation. 
From a practical standpoint, however, only those procedures should be pro- 
grammed in which the volume of repetitive work justifies the initial cost of 
program preparation. 


TVA COMPUTERS 


Electronic computation in hydrologic studies in TVA began in September. 
1955 with use of facilities at the National Laboratory at Oak Ridge, Tennessee. 
In December, 1957 the TVA Division of Water Control Planning obtained a 
small Royal-McBee LGP-30 computer. In June, 1958 installation of a high 
speed IBM 704 computer was completed. Hydrologic computations have been 
programmed for both TVA machines. ie : 

Programming for computers in TVA is on a user basis. That is, programs 
for hydrologic computation are written by engineers engaged in hydrologic 
studies. Programs dealing with system operations are written by engineers 
engaged in river control activities. To this end, selected engineers and engi- 
neering aides from the Hydraulic Data Branch, and every other using Branch, 
have been trained in computer programming. 
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Hydrologic problems programmed for the small LGP-30 computer were 
those which arose during basic processing of information, such as computation 
of mean daily discharge from time-gage height abstracts, computation of daily 
suspended sediment loads from discharge and concentration data, and compu- 
tation of indices of antecedent precipitation. 

Hydrologic problems programmed for the large IBM 704 computer were 
those which contain complex or voluminous computations. Two programs have 
been very useful. One computes daily values of a soil moisture index from 
rainfall, runoff, and evaporation for extended periods of record. The second 
is a program for comprehensive multiple regression analysis based on the 
method of least squares. 


THE MULTIPLE REGRESSION PROGRAM 


The multiple regression program is very definitely the most valuable of 
the current hydrologic computer programs. The method of least squares as a 
technique for evaluating the simultaneous effect of many factors has become 
well established in hydrology. Coupling this technique with regression analy- 
sis makes the scope of analytic results almost unlimited, and appears to offer 
the best approach to hydrologic evaluation whichis a necessary part of the TVA 
Tributary Watershed Studies. 

The multiple regression program used in TVA was written to include five 
major phases of analysis in fitting prediction equations to observational data. 
These will be described in sequential order of computation. 

Product Summation.—Input data for the program consist of individual obser - 
vations and associated parameters which have been computed for each obser - 
vation. Generation of a triangular matrix of product-summations is the first 
hase of computation. This matrix is the setof coefficients of the simultaneous 
10rmal equations for least-squares fitting. A 24 x 24 matrix on 500 obser- 
vations can be generated. The triangular matrix is printed. 

Simultaneous Equations.—A symmetrical square matrix is generated from 
he triangular matrix. This matrix is inverted and the regression coefficients 
are computed. Both are printed. Following this, various statistical quantities 
resulting from fitting the equation are computed and printed. These include 
sums of squares for regression and for residuals, the correlation coefficient, 
he error of prediction and the error of regression, the standard errors of the 
regression coefficients, and the arithmetic mean of each variable. 

Predicted Values and Prediction Errors.—After the statistical results of 
‘itting are printed the observations are again read into the computer. The pre- 
licted value for each observational set is then computed from the regression 
-oefficients and the independent variables. The predicted value is subtracted 
rom the observed value to give the prediction error for each observation. 
Both predicted value and prediction error are printed. 

Predicted Values and Variance of Prediction.—To this point in the re- 
rression program the operator has not been able to exercise any options or 
rogram alternatives on computation or output. The last two phases of compu- 
ation are optional. ; 

It is sometimes desirable to test the statistical variability of the regression 
quation whenused as aprediction device. The prediction errors computed for 
he initial data set do not meet this requirement. However, standard statistical 
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procedures can be found for this computation. 2 This procedure, an operation 
on the inverse matrix, was written into the program. 

To use this optional portion of the program a list of synthetic observations 
must be prepared in advance, and punched oncards or placedon magnetic tape 
along with the list of observations used in fitting. The only requirement in set- 
ting up the synthetic list is to set the independent variables at values at which 
it is desired to test the variability of prediction. The machine computes the 
predicted value and the variance of the predictionfor each synthetic observation 
in the set. For a simple X-Y correlation this portion of the regression pro- 
gram would compute the familiar hyperbolic confidence limits about the line of 
regression. For many independent variables the model is complex and the 
arithmetic becomes quite lengthy. More than 200 product-summations are re- 
quired. For example, to evaluate the variance for a single prediction in a 10- 
term equation, and the computer is a tremendous time saver in this operation. 

Significance of Variables.—The last phase of the computer program was 
written touse in testing the different independent variables for statistical sig- 
nificance. This phase of the program also is optional. 

It is true that the standard errors of the regression coefficients can be used 
to form a test of significance on the statistic “t”. This test, however, is limit- 
ed to one variable. When several variables incombination are to be tested for 
significance, the regression equation must be re-evaluated with this combi- 
nation of terms left out. The resulting change in the sum of ‘squares absorbed 
by the regression forms the basis for analysis of variance, using the statistic 
“F” to determine the significance of the entire combination of variables under 
consideration. 

The TVA multiple regression program allows any number of re-evaluations 
of the regression equation with any number or any combination of variables left 
out. The combinations of terms dropped are independent of the order in which 
the re-evaluations are made. With this complete freedom of choice, the ana- 
lyst specifies beforehand the various re-evaluations which are desired. For 
each reduced equation the program causes computation of the reduced inverse 
matrix and regression coefficients and the other various elements as deter- 
mined for the first evaluation. The predicted values, however, are not comput- 
ed. Comparison of the solutions for all the reduced equations, with the initial 
total equation can then be made on all these elements of the solution. These, 
of course, are in addition to the conventional analysis of variance. 


EQUATION FOR PEAK DISCHARGE 


Some examples of the equations being evaluated by the computer program 
for regression analysis may be of interest. It is, perhaps, a seeming contra- 
diction, that proper utilization of computers requires a better working know- 
ledge of mathematics than when investigations are based purely on empiri- 
cisms and short-cut assumptions. However, computer programs must follow 
precise mathematics if part of the goal of computer utilization is more precise 
analysis. - 

One particularly vexing problem in TVA hydrologic studies has been the 
prediction of peak discharges for some small, single-practice research water- 
sheds in western North Carolina. The watersheds are operated as a coopera- 


2 “STATISTICAL THEORY IN RESEARCH,” 


by R. L. Anderson : 
McGraw-Hill Book Co., New ;York, 1952, p. 202 po Sek esha 
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ive project between TVA and North Carolina State College. Four watersheds, 
ranging in size from 3.7 acres to 5.6 acres are being operated under astatisti- 
sally designed experiment. Results from four different covers on the four 
watersheds will be studied by analysis of variance following a4 x 4 latin square 
Jesign. Preliminary to this analysis it will be necessary to establish the re- 
lationship of the important hydrologic and meteorological variables, so that 
jata can be adjusted for varying meteorological conditions. 

On these small watersheds the normal hydrologic procedure of correlating 
geak discharge with peak-producing rain, maximum intensity, antecedent 
moisture, and so one, would not produce workable peak equations. The runoff 
at the measuring flumes responds so quickly to changes in rainfall intensity 
that there are quite often several peaks during a storm, each peak represent- 
ing a burst of rainfall. After several trials at analysis, it was decided that the 
pattern or time distribution of rainfall would have to beconsidered. The follow- 
ing development leads to such an equation. 


TIME IN 
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FIG. 1.—DEFINITION OF RUNOFF PERIODS WITH RELATION TO 
PEAK DISCHARGE 


Any ordinate of a storm hydrographcan be predicted froma unit hydrograph 
ind incremental volumes of runoff. Peak discharge, being one ordinate, can 
e predicted by an equation which expresses the unit hydrograph principle: 


Q, = U,(RO,) + Up(ROg) + Ug(RO3) + Ug(RO4) . - - + U,(RO,)- - (1) 


n which Qp is peak discharge, U represents the set of unit hydrograph coef- 


icients, and RO is the volume of runoff by periods. iit 7 
The quantities U in Eq. (1) can be considered time -distribution coefficients 


nd they do not have to be at uniform intervals.3 For the western Pa ea 
yatersheds the coefficients were assigned a time definition as shown in Fig. 1. 
n this schematic diagram rainfall and the resulting hydrograph are plotted to 
.common time scale. Zero time is taken as the time of peak, and the time 
lefinition of the coefficients extends backward to the beginning of rainfall. 

3 “HYDROGRAPH ANALYSIS BY THE METHOD OF LEAST SQUARES,” by W. M. 
myder, ASCE Proceedings Paper No. 793, September, 1955. 
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The volumes of rainfall, R, can be obtained from a recording rain gage. It 
is necessary to devise a method to predict runoff from rainfall for each of the 
four time periods. Runoff is an excess of rainfalland can be computed for each 
period from rainfall minus loss: 


ROpaeRe  Lpas sis eee ee (2) 


where loss is considered negative. 

Loss was considered in two parts. The first part was an initial amount 
which varied with season but which was independent of storm rainfall. The 
second part was a progressive amount which accumulated during the storm at 
a rate governed by the accumulation of storm rainfall, but also dependent upon 
initial soil moisture. A soil moisture index, used in this relation, was com- 
puted by modifications to a previously reported budgetary system which com- 
putes a two-level index.4 The deep-level index varies slowly, and was used 
as the expression for season. The upper-level index varies rapidly, and was 
used as a surface moisture index. This two-level index was calculated for 
each day by a separate computer program, using daily rainfall, daily runoff, 
and adjusted monthly pan evaporation. 

The equation used to express the total loss in terms of the two component 
parts was 


Ly = n(N + No) + m(M, + Mot) (Rt + Rot) ......... (3) 


in which Ly is the loss for eachof four time periods, M; represents the surface- 
layer soil moisture index at the beginning of each time period, N is the deep- 
layer soil moisture at the beginning of the storm, R denotes the rainfall during 
each period, No, Mo, and Roy are threshold values of respective variables, and 
n and m are relationship coefficients. 

Eq. 3 was used because it is the simplest possible expression of inter- 
action effect of two causative factors, plus the independent effect of a third 
factor.5 

The surface soil moisture index, M, may be considered as increasing di- 
rectly with rainfall. This is approximately correct when direct runoff is a 
small portion of the rainfall. From this is derived: 


My = Surface-layer moisture index at beginning of storm 
M3 = M4 + Rq 

Mg = M4 + Rqg + Rg 

My, = M4 + Rg + R3 + Re 


Eq. 4 can be used with Eq. 3 to define the loss for each of the four time 
periods. Only rainfall during the preceding periods, initial moisture index of 
the surface layer, and the moisture index of the deep layers are required. Com- 
bining the resulting equations for period losses with Eq. 2, these expressions 
for runoff can be substituted into Eq. 1. The resultant prediction equation for 
peak discharge is algebraically complex in its coefficients, but simple in its 
8 Ss Se See 


4 “COMPUTATION OF EVAPORATION AND EVAPOTRANSPIRATION FROM METE- 
OROLOGICAL OBSERVATIONS,” by M. A. Kohler. Paper prepared for presentation at 
AMS Meeting in Chicago, March 19-21, 1957. 

5 “SOME LINEAR MODELS FOR MULTIPLE REGRESSION ANALYSIS INHYDROLO- 
GY,” by W. M. Snyder. Duplimat Reproduction by Tennessee Valley Authority, Knox- 
ville, Tennessee, May, 1956. 


HY 2 STUDIES BY COMPUTERS 7 


variables. By collectirg coefficient terms the equation can be simplified to the 
following: 


Qp = a + by Ry + bg Rg + bg Rg + bg Rg + cy Ry My + cg Rg Mg 


+ cg Rg Mg + cg Rg Mg + dMyg+eN (5) 


In Eq. 5 the variables R, M, and N to which the peak, Qp, is related have 
been previously defined. The coefficients, however, result from algebraic sim- 
plification and have no real hydrologic interpretation. The numerical value of 
the coefficients can be determined by the methodof least squares. The results 
obtained by thus fitting the equation to data can be analyzed by multiple re- 
gression analysis. 

Eq. 5 has beenevaluated for several sets of watershed cover data from the 
western North Carolina project. Initial results for this interim stage of the 
project appear entirely satisfactory. Some numerical results are shown in 
Table 1. In this table multiple correlation coefficients and the conventional 


TABLE 1.—PEAK DISCHARGE REGRESSION ANALYSES FOR WESTERN NORTH 
CAROLINA COOPERATIVE RESEARCH PROJECT® 


Watershed 


Pasture | Corn Wheat | Pasture| Pasture | Pasture | Wheat 


Number of storms 97 68 65 172 106 62 
Multiple correlation 

coefficient 0.76 0.57 0.42 0.79 0.86 0.88 
Error of regression 0.014 0.256 0.118 0.009 0.018 0.155 


Average peak 0.069 | 0.742 | 0.374] 0.124 | 0.166 


aunits for errors and averages are cubic feet per second. 


statistical parameter, the error of the regression, are tabulated for compari- 
son with the average peaks which were observed under different covers. The 
multiple correlation coefficients show that, in general, the equation is making 
significant adjustments to the data. However, this is a relative measure and 
not too useful. A more absolute check is the error of regression. Table 1 
shows that the errors are small for pasture and over-grazed pasture, but that 
larger errors are associated with corn and wheat. ; 

Another way of studying the equationis by the magnitude of peaks which are 
predicted and the variance, or probable variation, of those predictions. Table 2 
shows predictions for a storm of high intensity with and without rainfall im- 
mediately preceding. The predicted peaks are reasonable for all covers. For 
both the short duration and long duration storms, the peaks are lowest for 
pasture, slightly higher for over-grazed pasture. However, predicted a 
are highest for corn for the short duration or summer intense storm, an 
highest for wheat for the winter intense storm. This follows the known me 
cycle of poor cover under wheat during the winter-time, following fall ok ing 
and poor cover under corn inspring andearly summer following spring p ing. 
: 
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The errors of prediction givenin Table 2 are simply the square roots of the 
respective prediction variances. The probable variations are given in this 
form so that the units are the same as for the predicted peaks. The errors of 
prediction show that, in general, the differences between peaks for various 
cover are real. However, some improvements must still be made. Variances 
for corn and wheat covers must be reduced. Special attention may also have 
to be given to Watershed 5. This watershed, under certain storm conditions, 
has very large amounts of subsurface flow, producing very different peaks than 
if flow were from surface runoff. As a consequence, peak discharge data from 
this watershed are more variable than from the other areas, and the average 
residual errors also tend to be large. In attempting to improve the prediction 
of peak discharges a study will first be made of the prediction errors for the 
individual storms. These data are available from the machine printout sheets, 
and an attempt will be made to associate the errors with cover or watershed 
condition as it varies with the annual cultivation cycle. 


TABLE 2.—PREDICTED PEAK AND ERROR FOR WESTERN NORTH CAROLINA 
COOPERATIVE RESEARCH PROJECT4 


NAETSIES 


Storm Parameters (Inches): R,=0.50, R= 0.20, R 3 RE 0, M=0.50, N=14.0 


Predicted peak 


Standard error of 
prediction 


Storm Parameters (Inches): RE 0.50, R= 0.20, R= 0.50, RE 1.00, M=1.00, N=20.0 


Predicted peak 


Standard error of 
prediction 


8Units are cubic feet per second. 


EQUATIONS FOR VOLUME OF RUNOFF 


An equation for prediction of volume of storm runoff from the small water- 
sheds in western North Carolina is also under development. This equation is 
very similar to that used for prediction of peak discharge. The Similarity lies 
in the definition of rainfall quantities by short periods of time. There is a dif- 
ference, however, in that the volume equation has the most precise time defi- 
ove at the beginning of the storm where loss rates are high and changing - 
rapidly. 

Fig. 2 is a diagrammatic definition of the periods of rainfall used in the — 
volume of runoff equation. As in Fig. 1, rainfall and the resulting hydrograph — 
are plotted to a common time scale. However, zero time in Fig. 2 is taken at 


the beginning of storm rainfall, and time definition extends forward to the end 
of rainfall. 
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The volume of runoff equation, based on the rainfall by periods as in Fig. 2 
and on soil moisture is: , 


RO = a + by Ry + by Ro + bg Rg + by Rg + bp Ry +cyM+cyn.. (6) 
in which RO is surface plus subsurface runoff, Rt are rainfall amounts dis- 


TIME IN 


RAINFALL 


¥ 
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Surface 
plus 
Sub-surface Runoff 
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FIG. 2._DEFINITION OF RUNOFF PERIODS WITH RELATION TO 
STORM RUNOFF VOLUME 


tributed by time periods, and M and N are moisture parameters described pre- 
viously. 

Only a few analyses of this equation can be reported at this time. Results 
of these fittings are encouraging and are shown in Table 3. Multiple corre- 


TABLE 3.—RUNOFF VOLUME REGRESSION ANALYSES FOR WESTERN NORTH 
CAROLINA COOPERATIVE RESEARCH PROJECT?@ 


Watershed 


Pasture 


Number of storms 114 46 
Multiple correlation 
coefficient 0.75 0.90 
0.025 


0.003 
0.029 


Error of regression 


Average runoff 


aunits for errors and averages are area inches. 


gression are low. Some study of 
there may be improvement possi- 
ve some critical intensity, 


lation coefficients are high and errors of re 
the residual errors of these fittings indicates 
ble by taking into account the volume of rain abo 
ag for each rainfall period. 


. 
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COSTS OF COMPUTATION 


It would be impractical to evaluate many equations of the distributed rain- 
fall type by desk calculator. The time required for solution of these lengthy 
equations would allow only a few trials. Many runs of the multiple regression 
program have been possible onthe 704 computer. A cost comparison was made 
after 17 analyses. Computer costs for these 17 runs totalled $663. This does 
not include cost of writing the program but does include costs of checking out 
the program onthe machine. It was estimated that the same analyses using 
desk calculators would have cost about $10,000. This estimate is based on an 
extrapolation of costs for comparatively short equations. There is no allow- 
ance made in this extrapolation for an increasing number of operator errors 
due to the tedium of lengthy calculations, nor to the possibility of increased 
personnel turn-over due to pure boredom. 


OUTLOOK 


The electronic computer will be used at an increasing rate in hydrologic 
studies. The high speed of computation makes possible the accomplishment of 
routine work quickly and efficiently. The working time of hydrologists and 
aides, released from much routine computation, can be used to develop new 
methods, new techniques, and new machine programs. The low unit cost of 
analysis when computed electronically allows rapid trialof these new methods 
over many observations and many data sets. 

Rapid strides are being made in development of techniques of data analysis. 
The great challenge, however, lies in hydrologic synthesis. Hydrologists will 
not be making full use of computer capabilities until techniques are developed 
whereby proposed watershed programs canbe evaluated by integrating the ef- 
fects of changeson many small sub-areas. This has beendone for water yield, 
based on infiltration classifications of the sub-areas. But it is important that 
practical techniques of integration be developed for changes in flood peaks, 
shape of the hydrograph, and transport of sediment. These are factors which 
can unquestionably be improved by good land management and should serve as 
the bases for evaluating proposed programs. 

Allphases of computer applications, basic processing of observations, analy- 
sis of data, and research for new techniques are necessary steps in a continu- 
ing program. The end result of this programis to make more and better hydro- 
logic information available to agronomists, foresters, and hydraulic engineers 
_ who are working to improve the welfare of the people. 
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SCALES OF VISCOUS ANALOGY MODELS 
FOR GROUND WATER STUDIES 


By Jacob Bearl 


SYNOPSIS 
Scales for viscous fluid models used for studying various problems 
of flow through porous media are described. Consideration is given to 
some model techniques and limitations. The discussion covers scales 
for length (including distorted scales), time, discharge, volume, storage 


coefficient, two liquid flows, and capillary effects. 


INTRODUCTION 


The movement of fluid through porous media is a complicated phenomenon, 
and the mathematical equations that describe it can be solved only for a few 
very simple cases. Difficulties are removed to a large extent when model 
studies replace the theoretical approach. Among the various types of models 
which are used for studying ground water movement are the sand model, the 
electric analogy, the heat analogy, the membrane analogy, and the viscous 
flow analogy. This paper is concerned with the viscous flow analogy, also 
known as the Hele-Shaw model. 

The viscous flow analogy model consists of two parallel plates mounted 
together with a capillary interspace (0.5-3.5 mm) between them. The analogy 
is based on the similitude between the differential equations which describe 
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the field of saturated flow of water through porous media and those for lami- 
nar flow of a viscous fluid through the capillary interspace between two 
parallel plates. 

The first model of this kind was developed by H. S. Hele-Shaw2,3 in 1897- 


1899, for the study of flow patterns around variously shaped bodies. 
R. Dachler4 in 1936, was the first to suggest the application of this model to 
ground water flow studies. Since 1940, the model has been much improved, 
developed, and intensively used in the Hydrological Laboratory of the Govern- 
ment Institute for Water Supply, The Hague, Netherlands, for studying various 
problems of ground water movement5,6,7,8,9; in the Israel Institute of Tech- 
nology, Haifa; and by Water Planning for Israel, Tel-Aviv. Other studies and 
applications have been reported by Todd10,11,12, 

The viscous flow analogy can be used for almost any study of two dimen— 
sional flow in porous media whether steady or not. It has been used for seep- 
age through dams and from channels, for tidal fluctuations of the water table, 
for ground water balances, for studies of the movement of the salt water- 
fresh water interface in coastal aquifers, and other problems. In 1956, this 
type of analogy was further improved to include the three dimensional case of 
horizontal flow in an aquifer with a fluctuating water table.13,14 

The principles of the viscous flow analogy are well known and described in 
many articles dealing with this type of model, yet a complete description of 
the computations of scales for the various cases is missing. This paper is 
presented as an effort to fill this gap. 


2 “Experiments on the Nature of Surface Resistance of Water and of Streamline 
Motion Under Certain Experimental Conditions,” by H. S. Hele-Shaw, Trans. Institution 
of Naval Architects, vol. 40, 1898, pp. 21-46. 

3 “Streamline Motion of a Viscous Film, Rept. 68th Meeting of British Assn. for Ad- 
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DERIVATION OF MODEL SCALES 


The two dimensional flow of ground water with a phreatic surface in a non- 
isotropic aquifer is described by: 


2 2 
Kxp Q %p a Kzp o p = Sop ¢p Sue ae) © 6 65 ete sine (1) 
x4, O24, Otp 
For a point on the water table, 
aby db ad, db, 
Kxp Ox, + Kyp Ozp = Sup = Np Oty ¢ Gente neh ak ees (2) 


if the capillary zone is neglected and there is no movement across the free 
surface (See Appendix I). The same two equations can also be written for the 
viscous flow in the capillary interspace between two vertical parallel plates 
of a viscous fluid model. The only difference stems from the fact that the 


model is isotropic, so that K,, = a g a . Hence, 


2 2 
Km(2 ?m 0 =) See d¢m 


+ 
Ox2, dz2, 


and 


= me an ¢ %m r Pac Ca tin ets ar ee 328 (5) 
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into Eq. 2 and comparing the result with Eq. 4, the following relations are ob- 
tained according to the laws of similitude: 


———— 5 oe ae = Nee salts) arolake ene 
Kxr?f K, 62 Ker ?y ny ?y 


Here a can be any number. From the definition of a it follows that %,Z,; 
therefore, the model scales are given by; 


x (2a ee (7) 
Kxyy Ker Zr ~Kzr Kxp 
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and, 


2 
x 
_ Or Zy rae sa shunts (8) 
tT Ker 
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LENGTH SCALES 


Eq. 7 describes the distortion in length scales in cases of nonisotropic 
soils. When Kxp = Kzp the horizontal and vertical scales must necessarily be 
the same. Very often a model must imitate an aquifer having a length much 
longer than its depth (for example, several kilometers as compared to a hun- 
dred meters). It is impossible to design such a model without distortion of 
scales. In such cases, it is necessary to assume that the hydraulic conductiv- 
ity in the horizontal direction is much larger than it is in the vertical direc- 
tion, and then the distortion of the scales will be determined accordingly. The 
error which might thus be introduced (in case of a disagreement between the 
assumed and the existing hydraulic conductivities) is usually small. One has 
only to remember that the flow is essentially horizontal with a relatively 
small vertical component due to the elongated shape of the aquifer. However, 
there are cases, such as in the vicinity of pumping wells or near the outflow 
into the sea, where this assumption might yield large errors, 


TIME SCALES 


The choice of the time scale will be governed by the time interval between 
two successive readings of discharge, levels, etc. Rewriting Eq. 8 ina more 
detailed form, 


2 
Nm Xr Kyp 12 Exp x 


= —_— —— 1 _ = 


t = —— 
Oe PRR Kee Oe ee 


This indicates that an appropriate time scale can be obtained for given length 
scales by manipulating the width of the interspace, b, and the viscosity of the 
liquid, ¥. Consideration of the interrelations between length and time scales 
will enable the most practical laboratory model to be designed for studying a 
given problem. 


DISCHARGE SCALES 


The discharge scale is obtained from Darcy’s law. Starting from the pro- 
totype flow, 


o¢ 
Qxp = - Kxp bp Zp Sz, 
and the model flow, 
O¢ 
Qxm = - Kxm bm 2m oa 
it follows that 
_ Kxr by za 


Oa 2 eee eee (10) 
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Similarly, 
pee ye Dy Kye seers ose sts ss ons (11) 
and, it follows that 


Qxr = Qzr = Qr 


VOLUME SCALES 


The volume scale is obtained from the continuity equation 
Vr = Q; ty = Ny Xy Zr by A aie its atte ay ec oe RO (12) 


Very often areas of open water such as rivers, lakes, or the sea have to be 
imitated by a model. In thesecases the model consists of three plates (Fig. 1) 
instead of two. A narrow interspace is kept between two of the plates while 
the third one, together with cutting the middle plate to the desired shape, 
serves to increase the interspace up to the sum of the capillary interspace 


Model Prototype 


Open water 


Back plate 
Enlarged interspace 


Middle plate 


Capillary interspace 


baa RS 
Front plate 


Qm=l 


FIG. 1 


plus the thickness of the middle plate. When the volume of liquid leaving the 
reservoir and entering the capillary interspace is of importance, this en- 


larged interspace is subject to an additional condition: 
The same volume scale must be applicable to both the narrow and the en- 


larged interspaces. With the notation of Fig. 1, the following equations can be 
written: for the soil 


and for the open water 


' gs Ny Xp Zy by 
Where the bar indicates values for the open water. 

Since 

Vr =Vr, and ny =1 


it follows that nyby = nyby 


a by : 
ee in is Shh agen 13 
or, 7 by ae i Oe ee et ONO Caer) a ; ( ) 


16 February, 1960 HY 2 
REPLENISHMENT SCALES 


The natural replenishment scale can be derived from the equation which 
approximately describes the flow in a phreatic aquifer with a horizontal bot- 
tom. For the case shown in Fig. 3 


6) dh __ og 
& (xn). pon 


gore ( On? ee 

aS Sx\dx) +P =" St 
where h is the thickness of the flow and P is the vertical recharge rate. Since 
hy = Zy it follows that 


or, with K, constant, 


2 
Ny Zy Kxr 2r 
Py = ty = r= eal deletes spots Sees (14) 
r 


NONHOMOGENEOUS AQUIFERS 


In case of a nonhomogeneous aquifer, that is, an aquifer which consists of 
several layers with different permeabilities, animitation by the model is pos- 
sible by varying the width of the interspace. It can be done by introducing thin 
sheets of plastic or any other suitable material with the geometric shape of 
the layers into the interspace, thereby leaving a narrower space free for the 
flow of the viscous liquid. 

Since the discharge scale must remain the same for the various layers it 
follows that K,,b, must remain constant. Hence 


Km bmn il be, bm 
= a = = ES — = a constant 
Kzp bp 12 © VKzp Bp 
or 3 
Dn 
La; CONStANE tac Was ceienes cacurinl anes, owe (15) 
p 
or 5 
bn 
i ae a constant 
xp 


It might be mentioned here that sometimes the condition b2m/ is used 
instead of Eq. 15. This follows from the assumption that Ky, must remain 
constant throughout the model whereas the right thing is to keep Ky b, con- 
stant throughout the layered model. 

A condition which causes some trouble when computing scales for a model 
imitating a layered aquifer is the condition that Np/ by, should be kept con-. 
stant throughout the model. This condition follows from Eq. 9 when ty is kept 
constant throughout the model. Very often this condition cannot be fulfilled, 
but the error introduced is generally small. In most practical cases, the most 
important layers with respect to hydraulic conductivity are the sand or sand- 
stone layers where variations in permeability or in porosity from one layer 
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to the other are relatively small so that the error introduced by neglecting 
this last condition is small. However, a larger error will be introduced in the 
case of more impervious layers. 

In the case of an impervious layer, b,, =0; therefore, the entire inter- 
space is filled. When layers of small permeability are present the computed 
bm for these layers might come out to be very small (for example, 0.1 mm) 
and cause difficulties in constructing and operating the model. These cases 
occur very often when relatively thin semi-pervious layers, in which the flow 
is mainly in the vertical direction, are present. A practical solution to this 
problem would be to use strips which completely seal part of interspace and ~ 
leave a wider interspace in the remaining part as shown in Fig. 2. 

The relationship between the partly sealed and the enlarged interspace is 
determined by the condition that the discharge Q9 through the whole length in 


FIG, 2 


the case of the narrow interspace by;, must equal the discharge Qj through 
the shortened length with the widened interspace bj. Therefore, 


from which it follows that 
Ea Says ee LenS eet a aa (16) 


Experience sfOws that in the case of relatively narrow interspaces the free 
space must be increased by about 15% over the computed values due to 


increased wall friction. The use of these strips was first introduced by the 
Hydrological Laboratory of the Government Institute for Water Supply, The 


| 
| 


_ Hague. 


STORAGE COEFFICIENT CONSIDERA TIONS 


Up to this point only models which do not imitate storage coefficient have 
been considered. It was assumed that specific storage coefficient can be 
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neglected in the case of the phreatic aquifer. In this case the storage coeffi- 
cient of the aquifer is identical to the effective porosity whereas the specific 
storage coefficient which results from the elastic properties of soil and water 
is so small that we can assume Sop = 0 and Som =0 in Eqs. 1 and 2 and they 
cannot be used for the determination of the time scale. 

The problem of imitating a confined aquifer with specific storage coeffi- 
cient, Sg, by a model has not yet been completely solved. Nevertheless, there 
are several possibilities some of which have been successfully tried. 

From Eqs. 1 and 3 it follows that the design conditions for a model of a 
confined aquifer should be: 


Zy Ky Zy 

CP ge ee e's 2S 0 Oe en eens ee (17) 
Xr Zr tr 

or 
2 
S 
th = et ye ee (18) 
Kzr 


If it is desired to consider the specific storage coefficient as well as the 
porosity in the case of anunconfined aquifer, the following condition must also 
be satisfied: 


Sor Zy = Ny © eee © ase 0 4 Cle we 6 «@ 6 (19) 


-An aquifer’s storage coefficient results from the elasticity of the soil and 
the water, but since the pressures in the model are relatively small we can- 


Pesta yt asi 


Phreatic table 


FIG, 3 


not use the elasticity of the plates and the viscous fluid for this purpose. The 
relationship between S and So is given by S = So D, where D is the thickness 
of the confined aquifer. : 
One of the several methods used for imitating the storage coefficient of a 
confined aquifer is to connect to the back plate a vertical pipe the cross sec- 
tion of which will be computed according to the imitated storage coefficient as 
sketched in Fig. 4. The aquifer length is divided into several parts each with 
a length of 1, (for example, 10 cm). The storage coefficient of each of these 
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2 


f 
: 


columns will be represented by a vertical pipe with a cross sectional area a. 
Hence 


and 


r Sp b1Sp 
When the specific storage coefficient, S), varies along the height of the col- 
umn, it is possible to connect to every opening in the back plate a separate 
vertical pipe each having a different diameter. In this case the entire area of 
the plates is divided into rectangles in which the storage coefficient of each is 


FIG. 4 


separately imitated. The same method can be used in order to increase the 
storage coefficient of a model which imitates a phreatic aquifer so that a 
more suitable time scale can be obtained. In this case 


blt+a (22) 


eas i ee Ge: owrertelye, 0) (@_i¢ Te so recneue ae, 


Sm 
TWO-FLUID FLOWS 


An important problem which can be investigated by the model is that of the 
simultaneous flow of two liquids. The problem of salt water intrusion into 
coastal aquifers is an example of one of these problems. In a model two vis - 
cous liquids (distinguished by different colors) can beused. Certain conditions 

ust satisfied for this situation. 
Fe ci tee acaiiie conductivity of the model is the same for the areas 
occupied by the two different liquids, and using the fact that pit sare 
conductivity of the soil is inversely proportional to the viscosity of the water, 


then 


1 oa Cees 
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and Kye = Ky 


v ee 
Hence, =) = (==) tcSabe pil en sis Br atta tetrad (23) 
p/f Vols 


where subscript f indicates the fresh water and subscript s indicates the salt 
water. 


Pie 


Fresh Water, 1; 


Salt Water, Vs 


FIG, 5 


The pressure at any point along the interface is the same on both sides of 
the interface. In Fig. 5, showing a coastal phreatic aquifer, at point A the 
equation for the prototype is 


or ap 5p = hp 
where $2 


Similarly, for the model 


am bm = hm 
and hence ay 6, = hy 
Because Bip Slip * heated tated teniPici soe at ghee 


it follows that the choice of the two liquids will be governed by the condition. 
that: 


1 ee A oe 
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so that 
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CAPILLARY CONSIDERATIONS 


The problem of the capillary rise of the viscous liquid which is flowing in 
the interspace must also be considered. The water table of this liquid does 
not represent a phreatic surface (on which the pressure is the atmospheric 
pressure), but the surface which results from the capillary rise. The height 
of the capillary rise depends upon the width of the interspace, the type of 
liquid, and the material of the plates. The phreatic table in the model is ob- 
tained by deducting the capillary rise between the plates from the liquid sur- 
face elevation, 

In ground water problems where the thickness of the saturated layers is 
much larger than the capillary rise, this phenomenon is of secondary impor- 
tance only. Yet, by choosing a suitable interspace, liquid, and vertical scale, 
the capillary rise in the model can be made to correspond to the capillary 
rise in the soil. 

There are other problems, such as seepage from reservoirs and channels 
in the first stages of the non-steady flow, where the capillary effect (in all 
directions) is of primary importance. Corrections can be made by a proper 
change in the driving head in the reservoir or the channel. 


Prototype 


FIG, 6 


As an example, consider the flow from an artificial replenishment reser- 
voir (two-dimensional case) shown in Fig. 6. The potential at point 1 is: 


@, =a+a=0 


The potential at point 2 on the wetting front is: $9 = (a + £) - he 


and ¢ the ordinates of points on the water 


where denotes the capillary rise 
e aie be written both for the model and the pro- 


surface. The same expressions can 
totype: 
dom = - (a+ he + o)m 
Aom =%im - Pom = (a + he + f)m = (a + ho)m + 5m 


$1p = 0 


bo, = - (a+ he + $)p 


Ady = Pip = OF = (a + he + o)p = (a + help + Sp 


ip pees 
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It follows that 


This means that the water surface in the model reservoir must be raised so 
that 
(a + he)m = (a + he)p 


When raising this level, special attention must be given to the fact that the 
geometric boundaries through which the seepage takes place are not altered. 
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APPENDIX I 


Let F(x,y,z,t) = 0 be the function which describes the phreatic surface. If 
P(x,y,z) is a particle on this surface at time t and P(x + Ax, y + Ay, z + Az) is 
the position of the same particle at time t + At when the phreatic surface is 
moving, then F(x + Ax, y + Ay, z+ Az, t+ At) =0. The change can be ex- 
‘pressed by the total derivative in the form: 


DF _OF OF Ox OF Oy | OF Oz 
Dt Ot * Ox dt * dy dt * dz St 
The potential @ = z + . for the phreatic surface (p = 0) takes the form 
O= 2, or F=f = % = 0. 


Hence, 


Ox O Ox ; 
Now at’ of , and at describe the components of the velocity u, v, w in the x, 


y, and z directions, respectively. 


According to Darcy’s laws 


qa grad: @) anelaeatigsh ortee 44h esi (2) 
hence 
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a = - K, $f =n, 
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where n is the effective porosity and all subscripts indicate directions. 
By inserting (3) into (1), the resulting final equation is 


K.(22)° + x, (32)? + x, [(22)? - (2¢)] «2 9¢ 


APPENDIX II—SYMBOLS 


a = pipe area; 

by; = Width of interspace; 

bp = width of investigated area of prototype; 

g = acceletation of gravity; 

h_ = capillary rise; 

K° = hydraulic conductivity; 

K _ = hydraulic conductivity in the x-direction; 

Kz = hydraulic conductivity in the z-direction; 

m = subscript, indicates model values; 

P = vertical recharge rate; 

p = pressure; 

p = subscript, indicates prototype values; 

Q = discharge; 

r = subscript, indicates the model-prototype ratio; 

So = specific storage coefficient defined as the volume of water released 
from storage in a unit volume of aquifer by a unit decline of head; di- 
mension L; 

S = storage coefficient, defined as the volume of water released from a 
column of aquifer (the height of the whole aquifer) having a unit cross- 
sectional area by a unit decline of head. It is dimensionless; 

t =) time: 

V_ = volume; 

X, y, Z = coordinate axes; 

@ = constant; 

Ys, Yp = Specific weight of salt and fresh water, respectively; 

6 = ratio of the difference between the specific weights of salt and fresh 
water to the specific weight of fresh water; 

py = kinematic viscosity; and 

@ = ground water potential at every point. 


2366 February, 1960 HY 2 


Journal of the 


HYDRAULICS DIVISION 


Proceedings of the American Society of Civil Engineers 


BOUNDARY LAYER STIMULATION IN RECTANGULAR CONDUITS 


By R. G. Cox,! M. ASCE and F. L. Bauer? 


SYNOPSIS 
The effects of artificial stimulation of the turbulent boundary layer 
were studied in a 0.283-ft by 0.5-ft plastic conduit. The length of con- 
duit required for fully developed pipe flow was reduced from 18.5-ft to 


14-ft by 1/8-in. cubes cemented to the conduit entrance curves. 


INTRODUCTION 


This paper describes an investigation of the effects of artificial stimulation 
of the turbulent boundary layer in rectangular conduits and indicates the ulti- 
mate use of boundary layer stimulation in models of hydraulic structures. 

A review of the problem resulting in the study is of interest. For the past 
seven years the Waterways Experiment Station has cooperated with the United 
States Army Engineer Districts in field tests on large hydraulic structures, a 
number of which have been made on sluices and conduits of flood control dams. 
The elements of a typical conduit, consisting of the curved entrance, the con- 
trol gate section, and the conduit proper, are shown on Fig. 1 for a concrete 
gravity dam. Conduits may also contain transitions downstream from the gate 
section, vertical and horizontal bends, and exit portal constrictions. As flow 
passes through each section, it loses energy. For convenience, the energy 

losses occurring upstream from the uniform conduit section are grouped to- 
gether and termed the intake loss. This loss is considered to be the total 


Note.—Discussion open until July 1, 1960. To extend the closing date one month, a 
ive Secretary, ASCE. This paper is part 


written request must be filed with the Execut ‘ r 
of the copyrighted Journal of the Hydraulics Division, Proceedings of the American So- 
piety of Civil Engineers, Vol. 86, No. HY 2, February, 1960, 
‘ Chf, Analysis Section, Hydr. Analysis Branch, Hydr. Div., U.S. Army Engr., Water- 
ways Experiment Sta., Vicksburg,’ Miss. 
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FIG, 1.—ELEMENT OF A TYPICAL CONDUIT 


POOL ELEVATION 


FRICTION PRESSURE 
GRADIENT 


FIG, 2.—ENERGY LOSSES IN A TYPICAL CONDUIT 
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available energy head H minus the conduit velocity head and the friction loss 
hg and miscellaneous losses in the uniform conduit section (Fig. 2). The in- 
take loss, as defined, includes the losses attributable to the entrance curves, 
fo the gate section and slots, and to the tunnel transition. It may be expressed 
as a function of a dimensionless coefficient by 


Ke = V2/2e a tS Ree es eae ee (1) 


4 , 2 
in which the energy losses (intake) he = H - hg - a and V2/2g is the velocity 


head (conduit proper). In every case where comparison between model and 
prototype results has been possible, the prototype loss coefficient is consider- 
ably greater than that observed in the model. For example, the observed Pine 
Flat Dam sluice prototype coefficient is 0.16. A comparable laboratory study 
indicates a coefficient of 0.07. Similar results were noted for Denison Dam 
flood control conduits where coefficients of 0.12 in the model and 0.19 in the 
prototype were observed. 

Models of hydraulic structures are usually designed on the basis of geo- 
metric similitude and Froude’s law. The latter requires that the inertial and 
gravitational forces have the same relation in the model as in the prototype. 
Consequently, the model will have a smaller Reynolds number than the proto- 
type and relatively greater friction losses. Therefore, the model is made as 
smooth as possible to obtain a closer friction relation to the prototype. It is 
believed that the turbulent boundary layer may develop more rapidly at the 
very high Reynolds numbers that occur in the prototype than at the low Reynolds 
numbers that occur in the model. Further, the entrance, gate slots and tran- 
sition of the prototype couldestablish a higher turbulent energy level than could 
be sustained by wall roughness alone. In the case of the model, the boundary 
layer may develop more slowly and may not reach a turbulent energy level 
comparable to that of the prototype. If such is the case, artificial stimulation 
of turbulence in models may be necessary to reproduce the intake loss coef- 
ficient of the prototype. 

In an attempt to understand the reasons for the apparent differences between 
modeland prototype intake losses, the Waterways Experiment Station has under- 
taken three approaches to the problem. First, theoretical studies of boundary 
layer development and wall shear stress on curved surfaces were initiated, 
and are still in progress. Second, in cooperation with the Omaha District, a 
horizontalstrut was installed in one of the 22-ft-diameter flood control tunnels 
of Ft. Randall Dam. This strut contained numerous instruments for measuring 
velocity distribution and turbulence fluctuations. Tests were made with flow 
up to 50 fps. In addition, a specially designed probe was used to investigate 
velocity distribution within one foot of the tunnel wall. And third, a laboratory 
study of the effects of artificial stimulation of the turbulent boundary layer on 
flow in a rectangular conduit was made. The latter study forms the basis for 


this paper. 


PREVIOUS STUDIES OF ARTIFICIAL STIMULATION 


s engaged intesting the drag on ship models 
have artificially stimulated the 


difference in boundary layer 


_ For a number of years, engineer 
at the David Taylor Model Basin, U. S. Navy, 
boundary layer to compensate for the relative 


j 


+ 
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growth in the model and prototype. Breslin and Macovsky3 report that without 
stimulation, this difference results inappreciable error in the predicted power 
requirements of the prototype at design speeds. This error has been attribut- 
ed to the fact that the slower development of the turbulent boundary layer at 
the low model Reynolds number results in smaller shearing stresses than does 
the more rapid development at the high prototype Reynolds number. Prior to 
1950 there was but little published literature describing attempts to investigate 
the effects of artificial stimulation on the boundary layer itself. The work of 
Breslin and Macovsky was a pioneer study in this field. 


THEORY OF TURBULENT BOUNDARY LAYER 


The theory of the turbulent boundary layer has been presented in numerous 
publications. Much work has been done on boundary layer development on flat 
plates and in circular conduits, and a limited amount on open channels of large 


cae caaee. 
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FIG, 3.—BOUNDARY LAYER DEVELOPMENT 


width-depth ratios. An extensive search of published literature revealed no 
record of any studies conducted on rectangular conduits of varying ratios of 
height to width. 

A brief summary of the theory of boundary layer development in a conduit 
is pertinent to the investigation to be described. It is generally known that in 
pipes the flow is laminar with Reynolds numbers less than 2,000; that a tran-_ 
sition range occurs between 2,000 and 3,000; and that above a Reynolds number 
of 3,000 the flow is turbulent. 

Boundary Layer Parameters.—Turbulent flow in a conduit develops by a 
progressive thickening of the boundary layer in the region close to the wall 
until this thickness is equal to the radius of the conduit at which point complete 
turbulence is established. In the case of the rectangular conduit this thickness 


eee ee 
3 “Effects of Turbulence Stimulators on the Boundary Layer and Resistance of a Ship 
Model as Detected by Hot Wires,” by J. P. Breslin and M. S. Macovsky, without stimu- 


lation, this difference results in appreciable error in the Navy Department, David Tay- 
lor Model Basin Report 724, August 1950. 
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would be one-half of the minimum conduit dimension. This phenomenon, for 
two-dimensional flow, is illustrated by Fig. 3, in which uy is the velocity in 
the undisturbed central core; u represents the longitudinal component of the 
velocity within the turbulent boundary layer;and 6 denotes the thickness of the 
turbulent boundary layer. 

Rouse? defines the nominal thickness of the turbulent boundary layer 6 as 
the thickness of the zone in which the major portion of the velocity variation 
takes place. Velocity distribution in the nominal boundary layer is illustrated 
by Fig. 4. 

Two important parameters used in the study of the turbulent boundary layer 
are the displacement thickness¢ 6 and the momentum thickness 9. The dis- 
placement thickness is usually defined as the distance which the streamlines 


FIG, 4.—VELOCITY DISTRIBUTION IN THE 
BOUNDARY LAYER 


have been displaced from the boundary as a result of the turbulent boundary 
layer. The displacement thickness4 can be expressed as: 


6 
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Aclearer concept of this parameter has been presented by Lighthill5 who states 


«, . . between the surface and any streamline just outside the boundary 
layer, there must be a constant volume of flow per unit span. This bor: 
be so if the flow reduction inside the layer is compensated for by an ou ‘ 
ward displacement of such a streamline through a distance 64 (whic 


. g) ; Inc., N. ¥., 1950. 
RG ing Hydraulics,” by H. Rouse, John Wiley and Sons, Inc., 
5 SE archi een Thickness,” by M. a Lighthill, Cambridge University Press, 


. ournal of Fluid Mechanics, Vol. 4, Part 4, August 1948. 
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produces a flow increase U 6 4, since the velocity is U in the region of 
streamline displacement).” 


In the foregoing, 6 1 is the displacement thickness and U is the streamline ve- 
locity. 

The momentum thickness @ is generally defined as a measure of the defect 
of the momentum flux or decrease of energy within the boundary layer. The 
momentum thickness can be expressed as: 


6 
o- [ a Ce es RE or wn (3) 
oO uy uy 


Fig. 5 illustrates the physical and mathematical significance of these two pa- 
rameters. 

Other useful parameters in the study of the boundary layer are the local 
momentum Reynolds number: 


inwhich 7 is the kinematic viscosity of the fluid; and the local resistance coef- 
ficient as defined by Ross :6 


Cy = (4.4 + 3.81 log Ree snc ctee sa ae (5) 


DESCRIPTION OF TEST APPARATUS 


Tests were conducted on a rectangular plastic conduit 0.283 ft wide, 0.5 ft 
high and 30 ft long. The axis of the conduit was mounted perpendicular to the 
vertical face of a pressure tank. The transition from the pressure tank to the 
conduit was made by four elliptical surfaces formed to the equation: 


2 2 

X Y 
PR tie oe Th a 
(0.75D)2(0.25D) ¥ 


in which D is equal to the dimension of the conduit in the direction concerned. 

Flow was controlled by regulating valves and measured by venturi meters 
located between the pressure tank and the water supply system. The pressure 
tank head was measured by a mercurial U-tube manometer and referenced to 
the center line of the conduit. The conduit pressure gradient was determined 
by side center-line piezometers in the conduit proper connected to a manome- 
ter bank. Fig. 6 is a schematic plan of the general test installation shown in 
Fig. 7. 

Horizontal traverses at selected stations were made by means of a stag- 
nation tube. This tube consisted of a traversing rod 0.25-in. in diameter to 
which was fixed a 0.5 in.-long hypodermic needle, 0.028-in. O.D. and 0.016 in. 
I.D. The traversing rod and hypodermic needle were inserted in the conduit 
at the selected station prior to beginning of operation. Leakage around the 
traversing rod was controlled by means of packing glands. Fig. 8 shows the 
traversing rod and hypodermic needle inserted in the conduit. 


6 “Turbulence Flow in the Entrance Region of a Pipe, by D. Ross, Paper No. 54-A-89, 
ASME Preprint, 1955. 
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FIG. 8.—STAGNATION TUBE 
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The boundary layer was stimulated by blocks cemented to the curved en- 
trance as shown by Fig. 9. The blocks were 0.02-ft square with 0.02-ft clear- 
ance between each block. Separate tests were made on 1/16-in. and 1/8-in. 
high blocks. In each test the stimulators were fastened to all four surfaces 
of the entrance curves. 


TEST PROCEDURE 


The procedure followed was identical for each of the three series of tests 
(no stimulation; 1/16-in. stimulation, and with 1/8-in. stimulation). Tests were 
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FIG. 9.—ARTIFICIAL STIMULATION 


made with heads of 10 ftand 13 ftof water in the pressure tank for each degree 
of stimulation. The resulting discharge through the conduit was determined. 
The conduit side center-line piezometric pressures were read on the manome- 
ter bank to the nearest 0.005 ft. Velocity distribution data in the boundary 
layer were obtained by traversing the stagnation tube across the conduit. Lee 
first two readings were taken at 0.005-ft intervals from the conduit wall. Suc- 
‘cessive readings were made at 0.01-ft intervals to the centerline of the con- 
duit. The tube was always oriented to permit determination of the x-velocity 
component. Fig. 10 shows the piezometer locations. The velocity traverse 


stations are shown on Fig. 11. 
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DISCUSSION OF DATA 


In the interest of brevity, the discussion of data is limited to tests on flows 
having a Reynolds number of 6.02 X 105. This Reynolds number was computed 
using the diameter of a pipe equal to four times the hydraulic radius of the test 
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FIG, 11.—VELOCITY TRAVERSE STATIONS 


conduit. The results of tests on flows with other Reynolds numbers were simi- 
lar to the results presented in this paper. 


Velocity Data.—The local static pressure at each traversing station was 
computed from the pressure data measured at the side center-line piezometers. 
The stagnation tube and the local static pressure data were used to compute the 
velocity distribution across the conduit. A semi-log plot of the velocity data 
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with respect to the conduit wall was made. Two straight lines could be drawn 
through each set of plotted data: a sloping line starting at the conduit wall and 
a line perpendicular to the conduit wall. The sloping line defines the logarith- 
mic law of the velocity distribution in the turbulent boundary layer and the 
perpendicular line the velocity in the undisturbed core. The boundary layer 
thickness was determined graphically by the intersection of these two lines. 
Velocity traverses for unstimulated flow are shown on Fig. 12. The pro- 
gressive thickening of the boundary layer from station 1 to station 4 is shown 
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FIG. 13.—BOUNDARY LAYER STIMULATED 


iti i i i file lines. The absence 
by the position of the intersection of the velocity pro 
of the constant velocity line at station 6 indicates that the boundary layer is 
fully developed. This station is located about 50 equivalent diameters down- 
stream from the entrance. t 
Velocity traverses for stimulated flow are shown on Fig. 13. With sys ene 
‘block stimulation, the sloping line at station 5 extends to the center line fc) oe 
conduit. An almost similar condition exists at station 4 with 1/8-in. bloc 


H stimulation. 
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The boundary layer thicknesses at each station was determined from the 
plots of the velocity profiles. These data, plotted on Fig. 14, show the pro- 
gressive thickening of the turbulent boundary layers with and without stimu- 
lation. The effects of successive increases in the degree of stimulation are 
shown by the inward movement of the boundary layer and the decrease in con- 
duit length required for fully developed turbulence. 

Pressure Data.—The observed conduit pressure profile without boundary 
layer stimulation is shown on Fig. 15 (pressures referenced to side center line). 
This profile indicates that the friction pressure gradient extends to within11 ft 
of the conduit entrance. The steeper portionof the pressure gradient upstream 
from this point occurs in the regionof the conduit where the boundary layer is 
rapidly developing. The effect of stimulationon the pressure gradient is shown 
on Fig. 16 (pressures referenced to side center line). With each increase in 
the degree of stimulation, a corresponding decrease occurs in the elevation of 
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FIG, 14.—STIMULATION AND BOUNDARY THICKNESS 


the pressure gradient in the region of boundary layer development together with 

a slight upstream extension of the friction pressure gradient line. 
Friction.—The observed pressure gradient in the last 19 ftof the conduit was 

used to compute a friction factor by the Darcy-Weisbach formula. For flow 

with a Reynolds number of 6.02 X 105 an “f” value of 0.011 was obtained. This 

value plots below the smooth pipe curve on the Moody-type friction diagram. 
The equation for the wall-shear-stress coefficient, developed by Ross 
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for circular pipes was used to show qualitative effects of stimulation on the — 
wall shear stress. Application of this equation to the data from the unstimu- 
lated boundary layer resulted in a decrease in the wall-shear-stress coefficient 
from 0.0027 at station 2 to 0.0022 at station 6 where the boundary layer was 
fully developed. With maximum boundary layer stimulation the coefficients 
were found to decrease from 0.0026 at station 2 to 0.0023 at station 4, again 


HY 2 


PRESSURE, FT OF WATER 


RECTANGULAR CONDUITS 


FRICTION GRADIENT 


EXTENSION 


DISTANCE FROM ENTRANCE, FT 


FIG, 15.—PRESSURE PROFILE - NO STIMULATION 


LEGEND 


© NO STIMULATION 
oOo i716” BLOCK STIMULATION 
178" BLOCK STIMULATION 


FRICTION GRADIENT ~~ 


EXTENSION 


PRESSURE, FT OF WATER 


DISTANCE FROM ENTRANCE, FT 


FIG, 16.—EFFECTS OF STIMULATION ON 
PRESSURE PROFILES 


37 


38 February, 1960 HY 2 


where the boundary layer was assumed fully developed. While the Ross equation 
was developed from data for circular pipes, its application to a rectangular 
conduit is considered valid. 

Intake Losses.—The empirical intake loss coefficient defined earlier inthis 
paper is dependent upon extension of the friction gradient upstream to the in- 
take structure. A slight decrease in the slope and elevation of the friction pres- 
sure gradient was noted with maximum stimulation. It is believed that this re- 
sulted from an increase in the intake loss and a decrease in discharge too small 
to measure with available equipment. It is hoped that other tests can be made 
with greater stimulation and more refined equipment. 


CONCLUSIONS 


The following conclusions can be drawn from the investigation: 


a. Unstimulated fully developed pipe flow will occur at about 50 equivalent 
diameters from the entrance of rectangular conduits with a height-width ratio 
of 1.77 when operated in the range of Reynolds number of 108. 

b. Stimulation of the boundary layer will effectively reduce the conduit 
length required for fully developed turbulent flow to about 40 equivalent diame- 
ters. 

c. The friction pressure gradient indicates fully developed turbulence at 
about 30 diameters from the entrance. 

d. Computation of intake losses based on pressure gradients within 30 di- 
ameters of the entrance will result in erroneous indications of the intake losses. 
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SCOUR AT BRIDGE CROSSINGS 


By Emmett M. Laursen, 1 M. ASCE 


SYNOPSIS 
Relationships are proposed for the prediction of scour at piers and 
abutments for the case in which sediment is supplied to the scour hole. 
The relationships were obtained from a combination of an approximate 


analysis and laboratory experiments, and depend on knowledge of the 


flow conditions at the bridge site. 


GENERAL ASPECTS OF THE PROBLEM 


In order to design the foundations of a bridge over an alluvial stream, it 
is necessary to know the elevation of the stream bed in the vicinity of the 
piers and abutments. The elevation in question, however, is not the elevation 
at the time a survey happened to be made, but the lowest elevation which will 
occur during the anticipated life of the bridge. Unfortunately, engineers are 
equipped with slide rules rather than crystal balls; therefore this lowest 
elevation must be coupled with a probability of occurence. On this basis an 
economic analysis is possible with the cost of construction for one stream 
bed elevation lower than another considered as an insurance premium for the 
decreased chance of loss (including interruption of traffic). 

Causes of Scour at Bridge Piers and Abutments —A lowering of the stream 
bed in the vicinity of the piers and abutments can occur from a variety of 
causes. A useful distinction can be made by separating the various causes 
into two general categories; (1) those characteristic of the stream itself, and 


Note.—Discussion open until July 1, 1960. To extend the closing date one month, a 
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(2) those due to the modification of the flow by the bridge crossing. For the 
past ten years scour due to the modification of the flow has been studied at the 
Iowa Institute of Hydraulic Research under the sponsorship of the Iowa State 
Highway Commission and the Bureau of Public Roads. As a result of that in- 
vestigation means for predicting the scour at bridge piers have been evolved. 
The relationships which were obtained will be presented in skeleton form in 
a later section; for a full discussion of the investigation and the design rela- 
tionships, reference should be made to the final reports which have been pub- 
lished in bulletin form by the Iowa Highway Research Board.2,3 

Even without the complicating feature of a bridge crossing the bed and 
banks of an alluvial stream cannot be considered as fixed. Superposed on the 
very slow changes measured in geological time, there are changes which may 
occur suddenly in a single flood, periodically during a series of water years, 
or slowly during the life of a bridge. Past observations at the chosen site, or 
at similar sites, must be relied upon to estimate the scour which may occur 
in the future. Each river and each reach must be studied to understand its in- 
dividual, almost personalized, characteristics. 

The reach which includes the chosen bridge site may be degrading: 
naturally, as the erosional agent of the geological cycle, as a consequence of 
a dam some distance upstream, or as the result of stream straightening 
downstream. Old surveys and other records should yield evidence of natural 
degradation which can be used to predict the future. Degradation caused by 
the works of man can be analyzed, at least approximately, on the basis of 
sediment-transport relationships, that is, if the plans of man can be antici- 
pated. 

Meandering streams will tend to have the greatest depth at the outside of 
the bend and shallows between the bends. During high water, the bends will 
scour and the crossings fill; during low water, the roles will reverse. As the 
loops of the meanders grow, there may be natural cutoffs with resultant deg- 
radation of the reach above and aggradation of the reach below. Usually the 
meander grows, scouring the outside of the bend, filling the inside; it is not 
unusual during a flood, however, for a chute to develop across the inside of a 
bend and even replace the old channel. 

Braided streams are characteristically wide and shallow, with deeper 
channels which may shift erratically. Even if it is not believed that “if any- 
thing can go wrong, it will” is a natural law, it must be admitted that there is 
a fair chance that during the life of a bridge a deep channel may develop at a 
pier at the time of a flood. 

Any contraction, whether of the main channel or of the overbank flow, will 
result in scour during high stages. As shown by L. G. Straub,4 the principle 
of continuity applied to both the discharge and the sediment load permit this 
case to be solved analytically since the flow conditions in both the contracted 
and uncontracted reaches can be considered uniform. It is likely that this 
case of scour has given rise to the old riverman’s belief that during a flood a 
plains river will scour its bed as much as the water surface rises. That such 


2 “Scour Around Bridge Piers and Abutments,” by E, M. Laursen, and A. Toch, Iowa 
Highway Research Board Bulletin No. 4, May 1956. 


3 “Scour at Bridge Crossings,” by E, M. Laursen, Iowa Highway Research Board 
Bulletin No. 8, August 1958. : 


4 “Approaches to the Study of Mechanics of Bed Movement,” by L. G. Straub, Pro- 
ceedings (First) Hydraulics Conference, State University of Iowa, Iowa City, Iowa, 1940. 
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scour cannot be general throughout the stream has been well demonstrated by 
£. W. Lane and W. M. Borland.9 

The Long Contraction.—By describing the conditions in the uniform 
reaches above and in a long contraction by a discharge equation and a 
sediment-transport equation, it is possible to solve for the ratio of the depths 
of flow in the two reaches—and, therefore, the depth of scour. The solution 
will depend, in detail, upon the equations selected to describe the flow and 
transport. Fig. 1 is a definition sketch of the long contraction in which there 
is both channel and overbank contraction. The Manning formula can be used 
to describe the flow conditions 
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An approximate form of the total-sediment-load relationship recently pro- 
posed6 can be used to describe the sediment concentration 
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That portion of the shear associated with the sediment particle of size D is 
termed T.' and 7, is the critical shear for sediment movement. The coef- 
ficient b and the exponent a depend on the mode of sediment movement and 
are functions of the ratio of shear -velocity to fall-velocity. Only the exponent 
a is of importance in the final solution and the functional relationship can be 


expressed as: 
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For flood flow conditions, the assumption can be made without serious sudo 
that To'/T,. - 1= To'/Tc- At equilibrium Q = Qc + Qo and fy Q. = Cg At; 
algebraic manipulation then permits the depth ratio to be expressed as a 
function of the contraction parameters 
’ ” by E. W. Lane and W. M. Borland, Transactions 
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FIG, 1.—DEFINITION SKETCH OF LONG CONTRACTION 
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FIG, 2,.—SCOUR IN A LONG CONTRACTION 
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FIG. 3.—EQUILIBRIUM DEPTH OF SCOUR AT A PIER 
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Since the ratio of the n values should not be too different from unity and the 
power is at the most 0.37, this factor can be safely neglected. Assuming that 
the material scoured out during a flood is redeposited over a large area so 
that the depth of scour dg = yg - yy, the foregoing equation for an overbank 
construction reduces to 


and for a channel contraction to 
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These relationships are presented graphically in Fig. 2. 

The case of the long contraction is of interest not only for its own sake, 
but for the light its solution sheds on the nature of scour in rivers. The 
shear-velocity in the contracted and uncontracted reaches of a particular 
stream should not be very different; so that for a given mode of sediment 
movement the depth of scour for a channel contraction is, as a first approxi- 
mation, dependent only on the geometry of the contraction as described by the 
width ratio and the approach depth. The assumptions which give rise to the 
approximation are that the n value is the same in the two reaches, and that the 
ratio To'/ T, is large compared to unity. Although reasonable and justifiable, 
these assumptions result in the neglect of any secondary effect of velocity 
and sediment size. 

In the case of the overbank constriction, the ratio of the total discharge to 
the channel discharge and the depth of the approach flow are sufficient to de- 
termine the depth of scour, granted the two assumptions previously mentioned. 
In this case, the mode of sediment movement does not affect the depth of 
scour. The detailed geometry of the approach conditions will, of course, de- 
termine the flow distribution by determining the velocity and cross-sectional 
‘area of the flow in the channel and on the floodplain. However, it matters not 
whether the velocities are high or low, or whether the flow sections are large 
or small; but only what discharge ratio results therefrom. 
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A bridge crossing is in effect a long contraction foreshortened to such an 
extreme that it has only a beginning and an end. The flow at the crossing can- 
not be considered uniform, but the solutions for the long contraction can be 
modified to describe the scour at bridge piers and abutments with the use of 
experimentally determined coefficients. 


LOCAL SCOUR AT PIERS AND ABUTMENTS 


Experimental Investigation.—Following the floods of 1947 in which the 
State of Iowa experienced a considerable monetary loss because of thedam- 
age to or the failure of bridges, the Iowa Institute of Hydraulic Research em- 
barked on an investigation of scour around bridge piers and abutments. For 
the preliminary phase of the study in which the effect of the geometry of the 
pier or abutment was studied, a flume 10.5 ft wide and 35 ft long was con- 
structed. Two supply lines and two tailgates permitted the wide flume to be 
operated as two 5 ft flumes by the addition of a center wall. Since neither a 
sand elevator nor a trap was built into the flume, the time period of the runs 
was limited and only qualitative results could be obtained. The velocity and 
depth of flow of the standard run, however, resulted in general movement of 
the bed, and the observations which were made established the nature of the 
phenomenon that was occurring. 

Albeit qualitatively, the importance of the length-width ratio of the piers, 
the angle of attack of the stream against the piers, and the length of encroach- 
ment of the abutment was established in this flume, as well as the unimpor- 
tance of small details of the geometry of pier or abutment. It was also dem- 
onstrated experimentally that there was an equilibrium or limiting depth of 
scour. In the preliminary work, piers and abutments representative of cur- 
rent Iowa designs were studied. The study of geometry was continued’ with 
more generalized pier forms and scour arrestors. The effect of debris was 
was also investigated (debris can be considered as an enlargement of the pier 
somewhat less than the size of the debris mass depending on the permeability 
of the mass). 

For the second phase of the program in which the effects of velocity and 
depth of flow and of sediment size were to be assessed, a transport flume was 
constructed which was 5 ft wide and 35 ft long and equipped with a sand ele- 
vator and trap. Most of the runs were made with a representative Iowa pier 
set at an angle of 30° to the flow. The depth of scour was measured by an 
electrical scour meter which could sense the position of the sand by the dif- 
ference in conductivity between the water-sand mixture and plain water. As 
wide a range of velocity, depth, and sediment size was utilized as the charac- 
teristics of the flume permitted. Froude numbers approaching the critical 
and sand jumping the trap limited the maximum value of the velocity and the 
minimum size of the sand, respectively. The minimum value of velocity and 
the maximum size of the sand were limited by the requirement that there be 
general movement of the sediment as bed load. Nevertheless, a sixty-fold 
change in the rate of sediment transport was obtained—probably the most sig- 
nificant indicator of the range of the tests. The depth of scour at the front of 
_ the pier is shown in Fig. 3. No systematic scatter could be detected in the 
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measurements —indicating that the depth of scour was uniquely determined by 
the geometry. 

This transport flume was also utilized for a model study of a bridge pier 
on which field measurements were made. The measurements in the field were 
obtained with an adaptation of the laboratory scour meter. Comparison of the 
model and prototype data indicated that the depth of scour could be treated as 
simply another length and that the equilibrium depth of scour obtained in the 
field. 

Measurements on a system of multiple cylinders were also made. The 
most significant findings of these tests were that the depth of scour did not 
depend on the degree of contraction (or proximity of an adjacent cylinder) 
until the scour holes overlapped, and that the minimum depth of scour when 
interference did occur was given by the solution for the long contraction. 

For the third phase of the study, the first flume was modified by being 
lengthened 7 ft, and equipped with a sand elevator and traps. The sand was 
supplied to five ft on one side of the flume, simulating the river channel. 
Clear water was supplied to the other side of the flume which simulated the 
floodplain, or overbank area. The experimental program for the investigation 
of the scour at the abutment included a variation of the ratio of discharge of 
overbank and channel, the depth of flow, the position and angle of the pier, and 
the presence of vegetal screening along the bank line. A few experiments 
were also conducted on the effect of spur dikes off the end of the abutment 
parallel to the bank, and of screening on the floodplain in the approach to the 
abutment. 

Design Relationships For Abutments.—The observation previously men- 
tioned, that the local depth of scour does not depend on the degree of contrac- 
tion until scour holes around neighboring obstructions overlap, suggests an 
extension of the solution for the long contraction to the case of local scour. 
For a normal sand, the width of a scour hole at right angles to the flow is 
about 2.75 times the depth of scour. Fig. 4 is a definition sketch of a long 
contraction, which it will be assumed approximates the case of the abutment. 
The scour in the long contraction is assumed to be a fraction 1/r of the scour 
at the abutment dg. Rewritten for this special case Eq. 7 becomes 
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Since the depth of scour is unknown, Q, is unknown. Therefore, the discharge 
Qy Over an arbitrary width w is substituted where Qy/w = Q,/2.75 dg. (The 
width w should be approximately equal to 2.75 dg and in practice a trial-and- 
error procedure will be necessary). Algebraic manipulation will now result 


in 
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in which yg is the average depth of flow in the width w. Eq. 10, with a value 


of r of 4.1, is plotted in Fig. 5, together with the experimental data for runs 


with the approach fill and abutment normal to the direction of flow. Various 


set-back distances of the abutment and clear distances between the abutment 
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FIG, 5,—SCOUR AT AN OVERBANK CONSTRICTION 
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and the vegetal screen are represented among the data. The solid points are 
irom the last, and best, series of runs in which the abutments was not set back 
from the bank line. The high points are from early runs in which the velocity 
of the approach on the overbank was so great as to be unrealistic. The floor 
in the overbank area was later lowered to obtain a smaller velocity of ap- 
proach. 

For large values of the parameter Qow/Qwyo, it is apparent that the ap- 
proximate solution describes the scour remarkably well. For small values 
of the parameter, the discrepancy is due to a cross flow from the channel to 
the overbank in the approach to the constriction. This cross flow could be 
observed by means of dye streaks. Since the water-surface slope in the 
channel was greater than in the overbank area, the water surface at the head 
end of the flume was higher in the channel than in the overbank. Correcting 
the overbank discharge by an approximate calculation of the cross flow would 
shift the experimental points over to the analytic curve. This has not been 
done, to emphasize the importance of fully determining the overbank flow. In 
practice, the rate of flow on the floodplain will be determined by considering 
the slope of the water surface, the depth of flow and the area of the flow sec- 
tion, and estimating an n value. If the overbank flow is small, cross flow in 
the area immediately upstream of the constriction can give rise to a condition 
similar to that in the laboratory flume. If this cross flow is possible, but can- 
not be evaluated, the upper curve in Fig. 5 should be used. 

A similar approximate solution can be obtained for the case of an en- 
croaching abutment. Using Eq. 8a for the condition of bed-load movement 
(VgyS /w < 1/2), this relationship is 


d 1.70 
1 ons as ® +1) = blah ete eee dy 
Vo 7 Yo 


where & is the effective length of the abutment such that (denoting the flow 
obstructed by the abutment as Q,) 
aw bats 
LY 

Eq. 11 with r equal to 11.5 is plotted in Fig. 6 together with some experi- 
mental points. The square points were obtained during an investigation of the 
relief bridge problem and are for the case of a vertical, blunt-ended, normal 
wall. The round points are data from the investigation of multiple cylinders 
in which there was no interference effect. 

Among the variations in geometry tested was the angle of incidence be- 
tween the direction of the approach fill and the direction of the flow in the 
channel. The results are presented in Fig. 7 as a multiplying factor, or co- 
efficient, Kg to be applied to the depth of scour calculated for the normal 
crossing (6 = 90°). Although only determined for the overbank case, these 
coefficients should also be applicable to the encroaching abutment. A check 
on the points for the cylinders plotted in Fig. 6 would indicate that the effec - 
‘tive angle of the curved wall of the cylinder is about 45°. 

The effect of setting the abutment back from the normal bankof the stream 
is difficult to assess. In the laboratory experiments no measurable effect 
could be noted. The long contraction solution can be further modified to con- 
sider that a part of the overbank flow remains on the floodplain, and that the 


we 
WY 


q 


48 February, 1960 HY 2 


e) 
o 
0 
0 


5 


10 15 20 
£/yo 


FIG, 6.—SCOUR AT AN ENCROACHING ABUTMENT 
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FIG. 8.—EFFECT OF THE RATIO OF SHEAR-VELOCITY TO FALL-VELOCITY 
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scour hole moves back with the abutment. This approximate solution also in- 
dicates that the effect of set back is small. Both the experiments and the 
analysis, however, do not model a real set back in certain respects. Because 
of the bank, the depth of flow in the vicinity of the abutment will be less. If 
the set back is very large, the presence of the river nearby is immaterial. 
The condition is then a special case of a relief bridge and the scour is due to 
clear water. If the bank material contains considerable clay or is well 
covered with a good, well-rooted, tight turf the material near the abutment 
might be very resistant to scour. 

For the case of the encroaching abutment, the mode of sediment movement 
will affect the depth of scour. An approximate evaluation of this effect can 
be obtained by comparing Eq. 8a, b andc. For the same values of By/ Bo the 
depth of scour for the different modes of movement were calculated. Defin- 
ing K- as the ratio of the depth of scour under suspended-load conditions to 


depth of scour under bed-load conditions, Fig. 8 was prepared. Note that the 
fall velocity should be that of the bed material being scoured out. The so- 
called wash load should be ignored. The effect of a change in the mode of 
movement should be applicable also to the local scour at a pier. Experiments 
made with a circular cylinder and a bed material 0.04 mm in diameter con- 
firmed the order of magnitude of the effect. The depth of scour—very difficult 
to measure and extremely sensitive to unsteadiness of the flow—was about 
50% greater than for bed-load conditions. The ratio of shear-velocity to fall- 
velocity in this experiment was about 20, and the concentration between 5 and 
10% by weight, that is, 50,000 and 100,000 parts per million. The ratio be- 
tween the rate of suspended-load movement and the bed-load movement was 
about 500 to 1. 

Design Relationships For Piers.—All the available data on piers were ad- 
justed to scour around a rectangular pier aligned with the flow, and the design 
curve shown in Fig. 9 was drawn with conservatism and with due regard for 
the reliability of the various data. Also shown in Fig. 9 in Eq. 11 with r = 
11.5 (note that b = 22). The most important aspect of the geometry of the pier 
was the angle of attack between the pier and the flow, coupled with the length- 
width ratio of the pier. A family of curves is shown in Fig. 10 as a multiply- 
ing factor Kgyj, to be applied to the depth of scour obtained from the basic 
curve. If the pier is set at an angle to the flow, it is recommended that no al- 
lowance be made for shape. If it is certain that the pier is aligned with the 
flow and will remain so during the life of the bridge, a shape factor can be 
applied as shown in Table I. 

All the experimental work on which the design curve was based involved 
bed-load movement. If the mode of movement is different, the ratio of shear - 
velocity to fall-velocity factor K- in Fig. 8 should be used. 

Besides the local scour at the pier, it is possible that the scour centering 
at the abutment may extend its area of influence out to piers near the bank. 
Thus, a fraction of the depth of scour at the abutment may have to be added to 
the local depth of scour. Typical scour holes have been shown whereby this 
effect can be estimated.3 et 

In the bulletins of the Iowa Highway Research Board the application to de- 
sign is fully discussed, including schematic design examples illustrating the 
use of all of the design relationships. On the basis of these relationships the 

prediction of local scour at piers and abutments is simple and straightfor- 
‘ward. In order to use the relationships, however, the flow conditions at the 
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TABLE I 


Shape coefficients K, for nose forms. 
(To be used only for piers aligned with flow.) 


Nose form Length-width Ratio US 
Rectangular E 1.00 
Semicircular & 0.90 
eral 0.80 
Elliptic 
rye, il 0.75 
27a 0.80 


Lenticular 
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bridge crossing must be known. As any civil-hydraulic engineer knows, the 
forecasting of the river flow of the future is not easy. Few gaging sites have 
been operated for 100 years, and 100 years of record are not sufficient to de- 
fine precisely the magnitude of a 100-year flood. Unfortunately, even rarer 
floods during the anticipated life, say 50 years, of a bridge, have a probability 
of occurence which is large enough to be reckoned with. This problem, of 
course, is Simply the old basic problem of hydrology. 

In addition to magnitude, other characteristics of the flow to be expected 
must be foreseen; the stage, the direction of flow (which may vary with stage), 
the division of flow between channel and overbank (which will vary with stage 
and the stage of the floodplain). Although these details can never be foretold 
with complete confidence, they are essential to the prediction of the scour 
which will occur. The evaluation of the flow characteristics must be judicious, 
and a calculated risk must be accepted; otherwise, the worst conditions must 
be assumed and overdesign will usually result. If sufficient control can be 
exercised over the valley upstream from the bridge crossing, the risk can be 
reduced to the chance of the flood occurrences. Each bridge crossing will 
have to be treated individually, and in the final analysis economic considera- 
tions will dictate the solution. 

Backwater at Bridge Crossings.—When scour occurs at a bridge crossing, 
the flow pattern which then obtains does not resemble that of a two- 
dimensional constriction. The flow approaching the obstruction (pier or abut- 
ment and embankment) ducks beneath the surface and passes through the con- 
striction as a spiral roller. The remainder of the flow passes over this roller 
and through the constriction barely noticing the piers and abutments, or that 
there is a constriction. There is no contracting jet as in a two-dimensional 
slot; or as in an open-channel constriction with inerodible boundaries. Asa 
result, the velocity increases only as it must to transport the sediment load 
downstream from the bridge; the energy losses at and due to the bridge 
crossing are minimal. Since the flow is disturbed, there are, of course, some 
losses. Fig. 11 shows the measured rise in water surface in the abutment 
flume. This measured rise is partly due to the velocity-head change between 
the approach flow in the channel and the flow in the contracted section down- 
stream from the bridge. The average depth in the downstream section can be 
obtained from Eq. 7 and the curve for the velocity-head differential has been 
evaluated on this basis. The backwater at the bridge is the difference be- 
tween the velocity-head differential and the measured rise of the head of the 
flume. The backwater in a comparable inerodible flume according to Liu, 
Bradley, and Plate8 is also plotted for comparison. It is readily apparent 
that the increase in the flow section due to the scour reduces the backwater 
to a fraction of the amount that otherwise would occur. i 

It is likely that a major cause of backwater at a bridge crossing on an al- 
luvial stream might be the resistance to flow of the contracted stream in the 
reach below the bridge. Depending on bank height and cover on the floodplain 
the cross flow from the stream back to the floodplain might take place very 
slowly. If the resistance to flow is greater in the contracted reach than in the 
normal stream, backwater will result. Moreover, any backwater which does 
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occur in an-alluvial river tends to cause aggradation upstream, and a con- 
comitant lengthening of the backwater curve. 


SUGGESTIONS FOR FURTHER STUDIES 


Field Measurements.—In order to apply with confidence the relationships 
proposed for predicting the scour at bridge piers and abutments, field meas- 
urements are needed to verify the conformity of model and prototype. Some 
of these measurements should be made at sites of extremely simple geometry 
like that of the Skunk River pier near Ames, Iowa. Such sites can be easily 
modeled in the laboratory, and no questions should arise as to extraneous ef- 
fects. Other measurements should be made at sites of complex geometry 
which perhaps cannot be modeled in the laboratory. The scour can be pre- 
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FIG. 11.—BACKWATER IN AN ALLUVIAL CHANNEL 


dicted with the aid of the relationships presented, and any discrepancies may 
indicate factors that should receive further study. 

The experience of the past can be used to substitute for field measure- 
ments. By hindcasting, the scour which occurred around piers and abutments 
at selected crossing during major floods can be evaluated. If bridges that 
failed and that did not fail during the same flood are checked, a correlation 
should be possible among hindcasted depth of scour, permissible depth of 
scour (based on foundation design), and failure or non-failure. 

Laboratory Investigations.—The one general aspect of scour at bridge 
crossings which has not been discussed is the scour at relief bridges. Relief 
bridges are openings placed in the highway embankment as it crosses the 
floodplain of the valley. For design purposes the flow through the relief bridge 
must be considered to consist of clear water. If there is a sediment load it 
will tend to be much finer in size than the material which will be scoured out. 
The depth of scour due to clear water is most assuredly a function of the 
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velocity of flow and the sediment size as well as the geometry of the opening 
and its environs. The limiting condition for this type of scour isa boundary 
shear that is equal to the critical for the material of which the boundary is 
composed. Since there is no supply of sediment coming into the scour hole, 
at the limit there can be no sediment going out of the scour hole. This limit 
is finite, but is reached asymptotically in infinite time. It is very difficult in 
the laboratory to scale all factors controlling clear water scour and a combi- 
nation of experiment and analysis is indicated. 

For an abutment which is set back from the normal bank less than the 
depth of scour, the relationships proposed herein should be applicable. An 
abutment set back over three times the depth of scour should be considered a 
relief bridge. Following a successful investigation of the relief bridge prob- 
lem, it would seem in order to study the effect of set backs between these 
limits. 

One could, of course, go on studying the effect of various geometries of 
piers, abutments and overall crossing. However, the results of this investiga- 
tion indicate that the effects of these details are minor, and of less importance 
than the error that can be expected in the evaluation of the flow conditions at 
the site. Further studies of geometry should probably be deferred until field 
measurements indicate they are of sufficient importance. 
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SEDIMENT TRANSPORT AND DELTA FORMATION 


By E. Kuiper,! F. ASCE 


SYNOPSIS 

This paper deals with sedimentation computations that were carried 
out in connection with dike design for the Saskatchewan Delta Reclama- 
tion Project. The first part of the paper gives a description of the re- 
gime of the Saskatchewan River in general and of the Saskatchewan Delta 
in particular. The second part presents some observations on delta 
formation. The third part discusses the sedimentation problems that 
were encountered in the design of dikes. The fourth and final part pre- 
sents the method that was adopted to determine the future rate of sedi- 


mentation on the floodplains adjacent to the projected dikes. 


REGIME OF SASKATCHEWAN RIVER 


The Saskatchewan River system, Fig.1, originates in the Rocky Mountains, 
where some of the peaks exceed an elevation of 10,000 ft. above sea level and 
where the permanent creeks begin to flow at an elevation of 5,000 ft. to 7,000 ft. 
Most of the runoff in the mountains comes from the annual snowmelt in June, 
which may last for several weeks. During a normal snowmelt, there is no 
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sediment transport to speak of. The water inthe streams is clear andthe peb— 
‘bles and boulders on the stream bed do not move. Most of the sediment trans- 
port takes place during short, concentrated rainstorms. At such a time, the 
¢reeks become torrents and they bring down the debris from the mountain 
slopes in tremendous quantities. The majority of the mountain streams flow in 
V-shaped valleys with a steep gradient. They cascade over rapids and water- 
falls and are apparently in a state of active erosion. Some of the streams, 
however, flow in valleys with a mild gradient. These streams have a braided 
character. The valley floor is filled with gravel and small boulders and gives 
the impression of being a horizontal plane between the mountain sides, over 
which the numerous stream channels wander back and forth at will. Ofcourse, 
this valley plane is not quite horizontal, but it has a longitudinal slope which is 
sufficient to carry discharge and sediment load. 

Numerous small headwater streams leave the foothills of the mountains at 
an elevation of about 4,000 ft. Their bed consists of boulders and gravel; their 
slope is about 10 ft. per mile; and during normal snowmelt conditions, the wa- 
ter has a visibility from 6 ft. to 10 ft. 

At this place, one would normally expect the beginning of the alluvial fans 
that are often found at the foot of mountain ranges. However, in the Saskatche- 
wan Basin, it appears that the slope of the land east of the mountains is steep 
enough to carry the streams with their entire sediment load so that no deposi- 
tion takes place. In fact, the slope is so steep that active erosion takes place. 
Near the foothills, the stream valleys are small, but in an easterly direction, 
they increase in depth and width and unite into valleys of still greater dimen- 
tions. The visibility of the water rapidly decreases from practically clear to 
a visibility of about four ins. over a distance of about 50 miles, indicating that 
sediment is derived from local river-valley erosion. 

Between Macleod and Lethbridge, a distance of about 60 miles by river, the 
following observations were made during a normal snowmelt flood stage in 
June. The river valley had a width of about 1 mile, at places narrowing down 
to a quarter of a mile. The depth of the valley increased in downstream direc- 
tion from about 75 ft. to 125 ft., and occasionally to 175 ft. The river mean- 
dered from one valley bank to another, and where it touched, huge active land- 
slides could be seen. The valley banks consisted of shale overlain by glacial 
till. There were no signs of ancient alluvial fan deposits. The river had a 
varying width of several hundred feet and a varying maximum depth from 6 ft. 
to 15 ft. The bottom consisted of coarse gravel and boulders, well-sorted and 
well-rounded, indicating that occasional transport of this bed material takes 
place. While the boat drifted on the stream, the clear tinkling sound of coarse 
sand and fine gravel moving over the bed could be heard, the bottom of the boat 
acting as a soundboard. The slope of the river was about five ft. per mile. The 
suspended-sediment concentration was found to range from 400 ppm to 500 ppm. 
The low floodplain banks of the river channel were composed of layers of sand 
to coarse gravel, each layer well-sorted. At one place, an old oxbow was ob- 
served at a higher elevation than the present river channel, indicating that the 
river is in a state of degradation. Except for a few minor rapids, the river 
had a rather uniform maximum surface velocity of about 6 fps. This uniform 
condition may be due to the fact that the river valley is excavated in soft shale, 
the coarse bed material being transported from upstream. Thus, there is no 
cause for major rapids or waterfalls. The gradient of the river is probably 
governed on the one hand by the amount of bed-material load that comes from 
the mountains, and on the other hand by the rate of erosion that takes place in 
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the valley. It would seem that the eroding shale banks form an important 
source of the suspended-sediment load that is found in the river. 

Between Lethbridge and Medicine Hat, the river valley is similar to the 
foregoing description. The river channel and river valley gradually become 
wider and deeper, but the bed material remains composed of coarse gravel 
and small boulders, and the actively eroding valley banks are still composed 
of shale with an overburden of glacial drift. 

Between Medicine Hat and Empress, an interesting change takes place. The 
river banks increase in height until they are more than 400 ft. above the river 
level. There are no floodplains, and immediately adjacent to the river channel, 
the valley banks, composed of soft sandstone overlain by shale and glacial drift, 
rise up in weird formations that resemble badland topography. The slope of 
the river decreases from about 3 ft. per mile to 2 ft. per mile, and the maxi- 
mum surface velocity from 4 fps to 3 fps. The maximum depth ranges from 5 
ft. to 15 ft., and the width of the river from 700 ft. to 800 ft. The interesting 
feature of this stretch is the rather abrupt change inthe composition of the bed 
material. Near Medicine Hat, the river bed is composed of pebbles and small 
boulders. Halfway between Medicine Hat and Empress, the river bed is com- 
posed of small boulders, with intermittent stretches of sand and fine gravel. 
Near Empress, the river bed is composed of sand, with an occasional patch of 
gravel or a boulder sticking out of the sand. These last observations were 
made by dragging a sampler over the river bottom from a drifting boat, each 
time over a length of several hundred feet. The type of noise that the sampler 
made, plus the contents of the sampler, provided a reliable indication of the 
composition of the river bottom. 

Between Empress and The Forks, the river valley is from one to several 
miles wide and from one to several hundred ft. deep. There are still no flood- 
plains to speak of. The rolling valley sides, covered with grass and brush- 
wood, come down to the river edge. The river channel has a slope of about 1 ft. 
per mile and is from one to several thousand ft. wide. Within the rigid bound- 
aries of the valley banks, the river tends to braid. Even with a small boat and 
during the snowmelt flood, one frequently gets stuck in shallow water. The bed 
material consists of rather uniform sand with a mean diameter of 0.30 mm. 
There are no signs that erosion is very active in this stretch of the river. 

In the foregoing discussion, the course of the Old Man River and of the South 
Saskatchewan River was followed. The North Saskatchewan River has very 
similar characteristics. The Bow River carries less sediment, probably due 
to its milder gradient and lesser valley bank erosion. The Red Deer River has 
its drainage area mostly in the rolling plains east of the Rocky Mountains. 
Normally, the runoff from these plains is low, but during heavy rainstorms, the 
river may carry a considerable discharge with a relatively high sediment con- 
centration. This sediment is partly produced by the plains and partly by the 
river valley, which has in its middle part a very pronounced badland topography. 

Between The Forks and the delta, the gradient of the river increases to 
several feet per mile. Occasional rapids are found, and the bed material con- 
sists of gravel and boulders. The valley banks are from one to a few hundred 
feet high and are composed of glacial drift. At places, where the river chan- 
nel nudges into the valley banks, huge active bank slides over the total valley 
height, and several thousand feet long, can be seen. At the end of this stretch, 
the river channel descends about 50 ft. over Tobin and Squaw Rapids and en— 
ters the Saskatchewan Delta, shown in Fig. 2 oer 

About 10,000 yr. ago, receding glaciers allowed the Saskatchewan River to 
take its present course. Near its lower end, where the river entered glacial 
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Lake Agassiz, most of the sediment was deposited and the formation of the del- 
‘ta area, as it is shown at the present time, was commenced. Deposits of 
‘coarse material were confined to channel locations, while the fine silts and 
‘¢lays settled down farther away from the river. After the channels were ex- 
tended over some distance into the lake, while possibly being split up into 
several branches, a situation developed whereby the river would somewhere 
break through its own formed banks and choose an entirely new course. Inthis 
fashion, the Saskatchewan River changed from one group ofchannels to another 
for thousands of years. Some of the oldest channels are completely covered by 
subsequent sedimentation. Some of the more recently abandoned channels are 
still plainly visible on aerial photographs and show up in the landscape as 
densely wooded ridges, enclosing low-lying marsh areas. The whole delta 
area, disregarding local ridges and depressions, is a gently sloping plain 
about 30 miles wide and 120 miles long, with a gradient of about 1 ft. per mile 
near the upper end and 1/2 ft. per mile near the lower end. The single river 
channel that enters the delta cone has a slope of about 1 ft. per mile and may 
carry during normal flood stages a peak flow of 100,000 cfs, containing 0.3% by 
weight of sediment. The several river channels that leave the delta and flow 
into Cedar Lake, one of the remainders of the ancient glacial lake, have a slope 
of about 1/4 ft. per mile and may carry during normal flood stages acombined 
peak flow of 50,000 cfs, containing 0.1% of sediment. The difference in sediment 
load between the upper and lower end is used in raising the delta surface. The 
remainder of the sediment load is used in river-mouth extension into Cedar 
Lake. 

The average annual sediment transport into the Saskatchewan Delta is esti- 
mated at 12,000 acre-ft. About one third of this volume is sand with an average 
diameter of 0.20 mm., and the remainder is wash load (silt and clay). Prac- 
tically all of this sediment is transported in a suspended state and is deposited 
somewhere in the delta. During low and medium flood stages, the flow of wa- 
ter and sediment remains within the natural river banks and deposition takes 
place in the river channel and at the mouth of the river channels. During high 
and extreme flood stages, the river banks are overtopped and deposition takes 
place throughout the entire delta area. These three modes of sediment depo- 
sition: in the river channels, on the floodplains, and in the river mouth, will 
be discussed in more detail in the following paragraphs. 

The river channels at the apex of the delta cone have a slope of about 1 ft. 
per mile, which is in accordance with the observed regime slope on the lower 
North and South Saskatchewan Rivers, as noted earlier. The river channels 
near the foot of the delta cone at Cedar Lake have a slope of about 1/4 ft. per 
mile. This reduction in river slope corresponds with a reduction in sediment- 
transporting capacity. It is estimated from observations and supporting com- 
putations that the capacity to transport the washload is hardly affected, but 
that the capacity to transport the sand is reduced by about 80%. Asa result, 
sand is deposited in the river channels; the coarsest fractions upstream, and 
the finer fractions downstream. This is reflected in the composition of the 
bed material, which shows a mean diameter of 0.30 mm., 0.24 mm., and 0.19 
mm., respectively, near the apex, the middle, and the foot of the delta. 

_ During high flood stages, the river banks are overtopped and a flow of wa- 
ter and sediment is directed towards the marshes. Immediate ly beyond the 
edge of the river bank, the coarsest sediment begins to precipitate. Soil sam- 
ple analyses show that the natural levees are composed of the coarsest frac- 
tion of the suspended sediment that is carried by the river. Away from the 
river, the deposited sediment becomes finer and finer, until only the wash load 
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FIG. 2.—SASKATCHEWAN 


remains which is deposited in the quiet marshes and lakes adjacent to the river 
channels. After a delay, which may range from a few days to several weeks, 
the overflow water, now deprived of its sediment load, returns some distance 
downstream via the outlet creeks of the marshes into the river channels and 
helps to reduce the average suspended-sediment concentration of the river 
flow. 

At Cumberland Lake and at Cedar Lake, river channels emerge into alarge 
body of water. Due to the inertia of the moving water, the river flows forward 
as a homogenous body of water into the lake and consequently, it retains its 
silt load for a while. After a certain distance, however, friction overcomes 
inertia, the water slows down, and sediment precipitates. This deposition of 
sediment takes place where velocities decelerate, which happens to be mostly 
at the flanks and in front of the retarding river current. Thus, two longitudi- 
nal bars and a number of shoals are formed. From aerial photographs, taken 
25 yr. apart, it appears that the bars and shoals keep growing in elevation and 
that the openings in between develop into river channels. Some of the less im- 
portant channels silt up until only one or a few main channels remain. All of 
the sand and most of the silt carried by the river is deposited near the river 


HY 2 DELTA FORMATION 61 


| 
| 
1h 
| 
| 
li 
4p 


DELTA 


mouth. Some silt and most of the clay is deposited further away in the lake. 
The bottom of Cumberland Lake and the western part of Cedar Lake, for in- 
stance, are almost entirely covered with soft mud. Since Cumberland Lake is 
only 10 ft. deep at flood stages, occasional winds keep stirring up the bottom 
deposits and as a result, the water that flows out of Cumberland Lake is always 
discolored. Cedar Lake, however, is larger and deeper, causing all sediment 
to permanently settle down, so that the outlet of Cedar Lake has clear water at 


all times. 
DELTA FORMATION 


ru 
“i The essential process in the formation of a river delta may be summarized 
as follows: A river brings water and sediment into a lake or sea. Through de- 


position of the sediment, the river mouth is extended. Since the river needs a 
ae to carry its water and sediment, the lengthening of the river channel will 
result in a rise in river stages at any given point on the lower reach of the 
river. There is a limit to how much a river can rise above the adjacent land, 


% after this limit is reached, the river breaks through its own banks, chooses 
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a new course, and beginsto repeat the former procedure. The manner in which 
this process takes place will be describedin more detail in the following para- 
graphs. 

When a river enters a lake the flow of water is not dispersed in all direc- 
tions, but due to its momentum, it maintains for some distance its original di- 
rection and moves as a well-defined stream into the lake. An appreciable 
amount of sediment, carried by the river, will be deposited at the sides of this 
stream, where moving and still water meet, and will form two banks that pro- 
vide further guidance to the river current. When after recession of the flood 
the lake and river levels drop, these banks may become exposed. This will 
enable them to consolidate and possibly to grow some vegetation. When the 
next flood comes, the same cause for further deposition on the banks is still 
present and may even be increased due to the vegetation, while the causes for 
erosion are lessened through the consolidation andthe presence of the vegeta- 
tive cover. When this goes on from flood to flood and no avulsion through the 
natural river bank takes place, the river mouth will extend like a finger in the 
lake. Sometimes, this form of river-mouth extension is very pronounced, and 
at other times, it is hardly noticeable. There are several factors which affect 
this phenomenon. First of all, the type of sediment carried by the river is im- 
portant. A large proportion of fine sediment will produce a relatively high 
sediment concentration in the overflow water and hence, a rapid sedimentation 
of the banks. Moreover, fine sediment will result in cohesive banks that may 
delay the occurrence of an avulsion and thus promote the extension of a single 
mouth into the lake or sea. Secondly, a rapid progress of the river mouth will 
result in an excessive concave upwardcurvature of the streamprofile. Inother 
words, the slope nearthe mouthbecomes relatively low. Hence, the mouth may 
protrude over a large distance into the lake or seabefore the height of the river 
banks becomes critical for an avulsion. There are several possible causes for 
rapid progress of the river mouth, like a shallow lake or a large concentra- 
tion of sediment. Since the slope of a graded river is determined by the bed- 
material load and not by the wash load, the circumstances would be ideal for 
finger-like extension when a river witha low concentration of bed-material load 
and high concentration of wash load would enter a shallow lake. The mouth 
would progress rapidly, the banks would be cohesive, and the slope would be 
small. 

After an avulsion takes place, the river will be using boththe oldandthe new 
channel for a certain length of time. However, there is a tendency for the old 
channel to silt up and to be abandoned for the following reasons. First of all, 
an avulsion often takes place in the outside of a river bend where the water is 
piled up during flood stages. This meansthat the new channel assumes a more 
or less straight course, whereas the old channel turns off to the left or to the 
right. Due to the tendency of the bed load to move towards the inside of a bend, 
the old channel will receive proportionally more sediment than water. In or- 
der to carry this sediment, the river needs a steeper slope, which it can obtain 
only by silting up its upper end. Even when the distribution of water and sedi- 
ment over the old and the new channel is in proportion, then the new slope on 
the old channel still needs to be steeper because two separated channels are 
less efficient in transporting sediment than one combined channel. 

In addition to the above causes for the actual silting up of the old channel, 
there is the tendency of the new channel to produce a local steep gradient at 
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the site of the avulsion. This is first of all because of the local difference in 
elevation between the river channel and adjacent land, and secondly, because 
the distance from the avulsion to the point where the sea or lake level is 
reached is usually less on the new channel, especially if deposits of the old 
river project far into the sea or lake. This steeper gradient in the new chan- 
nel near the avulsion will cause a draw-down water profile on the river up- 
stream. This will result in local scour and a further decrease in river stages 
near the entrance of the old channel, thus depriving this channel of more and 
more of its original flow. 

The deposition of sediment in the old channel will take place in such a way 
that the remaining cross section is in equilibrium, as discussed previously. 
That means that deposition will not only take place on the bottom of the chan- 
nel but also at the sides. The channel will become shallower and narrower. 
Finally, the entrance of the channel may become so narrow that it plugs up 
with debris or it may become so high that it becomes grown over with vege- 
tation during low-water periods. After this has occurred, the channel will be 
abandoned and be filled up with marsh vegetation. 

The rapidity of this silting-up process depends upon the proportion of ex- 

‘cess bed-material load that enters the channel. If this is low, for instance, due 
to the entrance being located nearthe convex side of the main current, the river 
will silt up slowly. In the meantime, the new channel also silts up due to nor- 
mal aggradation. Hence, the initial degradation upstream of the avulsion is 
checked and is reversed to aggradation. It may even be that the aggradation of 
the new channel catches up withthe aggradation of the old channel and that both 
rivers become equal partners. The aspect that these two rivers are less ef- 
ficient in their sediment transport than the former one river channel causes 
their slopes to be steeper. Consequently, the river channel upstream of the 
avulsion will aggrade above the level that it had at the time of the avulsion. 
After this has continued for awhile, another avulsion may take place, elimi- 
nating the two old channels and replacing them by one new channel, 

The formation and aggradation of the new channel, which initially flows 
through marshes and shallow lakes adjacent to the old channel, is somewhat 
similar to the formation of the old channel, as discussed in the above. How- 
ever, one point of difference is of interest. In the old channel, it was assumed 
that the river bed consisted of normal river bed-material. The new channel, 
flowing through low areas where formerly fine sediment was deposited, may 
find a bottom of cohesive material that resists scour. Hence, the formation of 
the cross section may be different from what it would be if the river flowed 
over its own bed-material deposits. It could be expected that in the beginning 
the channel will be relatively wide. After the river digs into the clay deposits 
part way, it could be expected that normal meandering or casual shifting will 
be less pronounced. It could also be expected that the river slope can be 
steeper than would be in accordance with the transported bed-material load, 
due to its resistance against scour. 

From the foregoing discussion, it follows that an avulsion may lead to the 
presence of two more or less stable and more or less active river branches. 

If this is repeated a number of times, the familiar picture of river arms in a 
delta, comparable to the veins in a leaf, is obtained. In downstream direction, 
the river channel will become increasingly smaller due to repeated bifurcat- 
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There is another reason why river channels in a delta decrease in size in 
downstream direction. During flood stages, when the river carries most of its 
sediment and when the formation of the channels takes place, water consistent- 
ly overflows the river banks towards the adjacent marshes and lakes. Asa 
result, the flow that remains in the channel will become increasingly smalle1 
in downstream direction. Since a channel is formed in accordance with the 
volume of its flow, it follows that the channel capacity will also decrease ir 
downstream direction. 

The consequence of this decrease in channel capacity is that only one flooc 
flow with a certain magnitude will cause uniform overflow over the delta. 
Floods of greater magnitude, which occur less frequently but which have a 
higher concentration of sediment and are of a longer duration, will flood more 
severely the upper part of the delta area; while floods of relatively small 
magnitude will overflow more extensively in the lower part of the delta area. 
In this way, the sedimentation is kept in balance. 


SEDIMENTATION PROBLEMS 


It was found from soil surveys that practically all the alluvial deposits ir 
the Saskatchewan Delta, covering an area of about 2,000,000 acres, are suit- 
able for agricultural development provided the water table can be controlled. 
An area of 150,000 acres, near The Pas, has recently been reclaimed, while 
plans are being made for the future reclamation of an additional 1,000,000 
acres of delta land. 

Reclamation of part of the delta surface will result in a smaller area being 
available for deposition of sediment and hence, in a faster rise of the surface 
level. This in turn will increase the flood elevations and thus endanger the 
projected dikes. The problem arises to determine the future sedimentation of 
the floodplain in quantitive terms and to find possible remedial measures. 

Generally speaking, there are three principal methods of coping with the 
present sedimentation problem. The first method is by reduction of the sedi- 
ment load through soil conservation. The second method would aim at increas- 
ing the sediment-transporting capacity of the river channel by river-training 
works. The third method includes all plans for controlled storage of sediment 
by reservoirs, or in the delta area. In the following paragraphs, the applica- 
bility of these methods will be discussed. 

In cases where the sediment, carried by the river, finds it origin in surface 
erosion of the drainage basin, soil conservation methods could be effective. On 
the one hand, they improve the usefulness of the land, and on the other hand, 
they may provide a solution to the sediment problem. However, the sediment 
carried by the Saskatchewan River finds its origin not only onthe prairie lands, 
but also in stream-valley erosion in the upper part of the watershed. To re- 
duce this type of sediment supply would call for extensive channel-protection 
works. Without entering into detail, it can be stated that the benefits in the 
delta area would be too small to justify the cost of such works. 

The method of increasing the sediment-transporting capacity by river- 
training works is applied extensively on rivers in Europe and occasionally on 
rivers in North America. In most cases, these works also serve the purpose 
of maintaining a navigation channel and protecting valuable shoreline property. 
In the Saskatchewan Delta, such works, which may cost anywhere from $50,000 
to $500,000 per mile of river channel, would be too costly to be justified by the 
benefits, at least for the time being. 
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Hence, if the sediment cannot be retained at the place of its origin, nor car- 
ried through to the mouth of the river, then it must be deposited somewhere 
‘along the course of the river, and the problem is to find the location that will 
do the least damage. Storage in an upstream sediment-control reservoir would 

provide an adequate solution. For the first few hundred years, the sedimenta- 
tion problem in the Saskatchewan Delta will not be so severe as to justify the 
construction ofan exclusive sediment-control reservoir. However, if upstream 
reservoirs are going to be built for other purposes, the feature of sediment 
control will be added automatically. It could even be expected that the conse- 
quent reduction in sediment load will cause degradation of the river channels 
in the delta, and thenceforth, will lower flood stages and improve drainage. 

If no upstream reservoirs are going to be built during the first few hundred 
years, the Saskatchewan River will deposit its sediment load in the delta. At 
the present time, the larger part of the sediment load is deposited in the Cum- 
berland Lake Delta. It was from a quantitative sedimentation study that this 
situation can be maintained for another 150 yr. From then on, increasing 
quantities of sediment will pass the Cumberland Lake Delta and will be de- 
posited in stream channels and on floodplains downstream. For some time to 
come, the only way to cope with the consequent rise in flood levels will be to 
raise the adjacent dikes accordingly. 

There is undoubtedly a limit to the possibility of raising dikes to keep pace 
with the natural sedimentation outside of the dikes. The ultimate slope of the 
river whereby it can carry the complete sediment load throughout the delta 
area is about 1 ft. per mile. A line under this slope, beginning at the future 
mouth of the river near the outlet of Cedar Lake, will reach a point near 
Cumberland House more than 100 ft. above its present elevation. It appears, 
therefore, that the area which is presently going to be reclaimed will have to 
be sacrificed sometime in the far future. However, it would seem entirely 
feasible to compensate for this loss by reclaiming the areas that have been 
silted up in the meantime. Thus, a situation would develop whereby every 
thousand years or so, another section of the delta area is made available for 
reclamation. When spread over periods of a thousand years, the cost of mov- 
ing the inhabitants from the area, which is going back to the river system into 
the area that is going to be reclaimed, appears very small and should be no 
objection to present reclamation. 

After having decided that the most appropriate remedial measure to cope 
with future sedimentation of the floodplain would be to raise the dikes as need- 
ed, the remaining problem was to determine the rate of sedimentation adjacent 


to the dikes. 


SEDIMENTATION COMPUTATIONS 


A serious attempt was made to compute in detail how the future deposition 
of sediment in the Saskatchewan Delta would take place. Precise sediment 
sampling was carried out at a dozen stations. Analytical computations were 
applied to determine the existing mode of transport and to estimate the rate of 
deposition of all sediment fractions for new equilibrium conditions. However, 
in the end, so many assumptions had to be made that the reliability of the an- ~ 


swer became highly problematical, particularly because it was found that slight 
changes in the assumptions, well within reason, could produce significant 
changes in the end result. It was, therefore, decided to determine the rate of 
future sedimentation, not on the basis of precise equilibrium caiculation, but 
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FIG. 3.—CROSS SECTION A-A 


on the broad basis of an estimated ave rage annual sediment load coming into 
the delta and an estimated distribution of deposition over the delta surface. 
Before these two items are discussed in some detail in the following para- 
graphs, it may be appropriate to point out at this place that detailed analysis 
and precise calculations were definitely helpful in exploring the nature of the 
sediment problem, and providing insight into the proper solution of the prob- 
lem. However, the present sediment problems appeared too complicated to in 
sist on precision to the very last answer. 

The average annual sediment load was estimated on the basis of sediment 
sampling carried out during 2 yrs. A relationship was established between 
river flow and total sediment transport. From this relationship and the flow 
duration curve, a sediment-transport duration curve was prepared, which in 
turn yielded the average annual sediment-transport figure of 8,000 acre-ft. 
per yr. Since the 2 yr. of observation were preceded by about 10 yr. of above-» 
normal rainfall conditions over the drainage basin, resulting in a better-than- 
normal vegetative cover over the sediment-producing prairies, it was decided 
to increase this figure by 50%, resulting in a total average annual sediment 
transport of 12,000 acre-ft. per yr. ~ 

The average distribution of sediment deposition over the delta area was 
estimated from a study of profiles and cross sections over the delta cone. One 
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such a cross section is shown on Fig. 3. Topographic surveys were carried 
out to determine the elevation of river banks and marsh lands. Hydrometric 
surveys were carried out to determine the distribution of flow during low, 
normal, and flood stages of the river system. Sediment sampling was carried 
out at various locations in the delta to determine where sediment was deposit- 
ed for certain flow conditions. Historical information was collected from 
journals of early explorers and from the memory of older residents. 

Based on the foregoing studies, the following conclusions were reached: 
At the present time, about 2,500 acre-ft. per year. of sediment load is carried 
by the Old Channel and remains in the river channel well past The Pas. The 
remaining 9,500 acre-ft. is carried by the New Channel and is nearly all de- 
posited in the Cumberland Lake Delta. Only 1,000 acre-ft., mostly wash load, 
leaves Cumberland Lake and is carried to Cedar Lake. Full development of 
the delta area includes reclamation of all land south of the Old Channel and 
Saskatchewan River, and will divert the Old Channel into the New Channel. By 
then, the New Channel will carry the full 12,000 acre-ft. of sediment load. 
Twenty years later, the bed-material deposits will have reached the outlet of 
Cumberland Lake and increasing quantities of bed material will begin to flow 
down the Saskatchewan River. Thirty years later, the Saskatchewan River, be- 
tween Cumberland House and The Pas, will carry its full capacity of bed ma- 
terial load and deposition will begin afterwards, starting at Cumberland House 
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and gradually extending in downstream direction. The foregoing sedimentation 
is graphically illustrated in Fig. 4. It was estimated that 100 yr. from now the 
extreme flood levels in the Cumberland House would be raised by about 2 ft., 
in The Pas area by 1/2 ft., and near the apex of the delta area by a negligible 
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amount. From then on, the extreme flood levels throughout the whole delta 
area upstream of The Pas will continue to rise from 1 to 2 ft. per 100 yr. 
Temporary and local deviations from this average rise may occur due to ab- 


normal climatic conditions or due to avulsions of the river channels. 
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TRAP EFFICIENCY OF RESERVOIRS, DEBRIS BASINS, AND DEBRIS DAMS2 


By Charles M. Moore,! M. ASCE, Walter J. Wood,2 F. ASCE 
and Graham W. Renfro? 


SYNOPSIS 
An attempt has been made to review and summarize available infor- 
mation pertaining totrap efficiency and to make itavailable to engineers 
and the general public in a unified form. As additional data are collected 
on these and other projects, it is hoped that it will be added to that pre- 
sented herein so that the subject of trap efficiency may be advanced as 


much as its importance warrants. 


INTRODUCTION 


Trap efficiency of reservoirs is a subjecton which there is a lackof readily 
available information for use by engineers in the design of water supply reser- 
voirs, flood control structures, and other kinds of impoundment or diversion 
structures. Trap efficiency is expressed as the ratio between sediment ac- 
cumulation and sediment outflow. One objective of the Task Force on Reser- 
voir Sedimentation, ASCE, has been to make a summary review of the knowledge 
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now accumulated, but not available in a unified form to the profession at large. 
The authors of this paper have attempted to do this. 

It appears that the most authentic information available in print is a paper 
by Gunnar M. Brune.4 This paper covers the normal ponded reservoirs very 
thoroughly as of this date. A summary of the paper is made herein, including 
the more important findings and results. In addition, a brief review is made of 
the cooperative studies of trap efficiencies of upstream floodwater retarding 
structures currently being made cooperatively by the U. S. Geological Survey 
and the U. S. Soil Conservation Service; and of the work being done by the Los 
Angeles County Flood Control District. There may be other published infor- 
mation of which the authors are not aware. 


SUMMARY OF PAPER BY BRUNE 


As pointed out by Brune, the trap efficiency of a reservoir depends upon a 
number of factors. Among these are the ratio between the storage capacity 
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From Transactions, American Geophysical Union, Vol. 34, No. 3, June 1953. 


° 25 50 7 100 125 150 75 200 
CAPACITY - WATERSHED RATIO ( STORAGE CAPACITY PER SQUARE MILE OF DRAINAGE AREA IN ACRE-FEET) 


FIG, 1.—TRAP EFFICIENCY AS RELATED TO CAPACITY-WATERSHED RATIO 


and inflow, age of the reservoir, shape of the reservoir basin, the type of out- 
lets and method of operation, the size grading of the sediment, and the behavior 
of the finer sediment fractions under various conditions. A number of attempts 
have been made tocorrelate trap efficiency with oneor more of these factors. 

One of the earlier studies of trap efficiency was made by Brune and R. E. 
Allen,® who developed a curve relating percentage of eroded soil caught ina 
reservoir with capacity per square mile of drainage area. The trap-efficiency 
values given by this method, however, are necessarily low because they are 
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based upon rates of erosion which, as pointed out by L. M. Glymph, Jr.,6 are 
higher than rates of sediment production to the reservoir. Glymph found that 
the rate of sediment productionis inversely proportional to drainage area and 
is usually lower than the rate of erosion because of the deposition of some of 
the eroded material such as colluvium, on flood plains and in stream channels. 
The ratio between storage capacity and inflow has been expressed in a gener- 
al way by the capacity-watershed (C/W) ratio. C. B. Brown” first developed 
acurve relating C/W and true trap efficiency. The curve shown in Fig. 1, 
with some additional records, is represented by 


Cr = 100 [1 - 1/(1 + 0.1 C/w)] 


where C7 is the reservoir trap efficiency in percentage, and C/W is the reser- 
voir capacity in acre feet per square mile of drainage area. There is a con- 
siderable spread in points, as may be noticed. 

The main reason for the spread of points in the curve is that reservoirs 
having the same C/W ratio may have very different capacity-inflow ratios, as 
previously stated by L. C. Gottschalk:8 


“A reservoir in an aridor semi-aridregion may have a low capacity- 
watershed ratio, yet not receive enough inflow in any one year to cause 
water to be discharged over the spillway. In contrast, the volume of 
mean annual flow from a watershed of equalsize in a humidarea may be 
equivalent to 25 times that of a reservoir having the same capacity- 
watershed ratio. In the dried region 100% of the incoming sediment load 
is trapped, whereas in a humid area possibly only 70% is trapped.” 


Allen Hazen? first introduced the storage, or capacity-inflow, ratio in 1914. 
He used it for determining reservoir requirements, however, and not as an in- 
dex of sediment trap efficiency. 

For individual reservoirs, curves can be drawn correlating trap efficiency 
with detention time in days. Whitney L. Borland has prepared1l0 such a curve 
for Imperial Dam Reservoir, shown in Fig. 2. Such curves are quite satis- 
factory for specific reservoirs, since other factors such as sediment charac- 
teristics, shape of reservoir, and method of operation tend to remain constant. 
Detention time may be defined as the period that the inflow remains in the 
reservoir before being released downstream. 

M. A. Churchill!! took into account both detention time and velocity of flow 
through the reservoir. He has developed a “sedimentation index” which repre- 


6 “Relation of Sedimentation to Accelerated Erosion in the Missouri River Basin,” 
by L. M. Glymph Jr., SCS-TP-102, U. S. Soil Conservation Service, Lincoln, Nebr., 
July 1951. thsi 

7 Brown, C. B., Discussion of “Sedimentation in Reservoirs” by B. J. Witzig, Proc., 
Amer. Soc. Civ. Eng., V. 69, No. 6, pp. 793-815, 1493-1499, 1943. 

8 «Analysis and Use of Reservoir Sedimentation Data,” by L. C. Gottschalk, pp. 131- 
141, Proc., Federal Inter-Agency Sedimentation Conference, Washington, D. C., Jan. 


1948, ; 
2 © t ided in Impounding Reservoirs for Municipal Water Supply, 
ee A : i 77, pp. 1539-1640, 1914. 


by Allen Hazen, Trans., Amer. Soc. Cinna ve 1 : - 
: 10 Unpublished data in files of Bureau of Reclamation, compiled by Whitney L. Bor- 


land, Denver, Colo., 1951. ‘ ; L 
11 Churchill, M. A., Discussion of “Analysis and Use of Reservoir Sedimentation 
Data” by L. C. Gottschalk, pp. 139-140, Proc., Federal Inter-Agency Sedimentation 


Conference, Washington, D. C., Jan. 1948. 
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sents the period of detention divided by mean velocity. His curve, relating 
trap efficiency to the sedimentation index for several TVA reservoirs, is shown 
in Brune’s paper. While this curve is very satisfactory where the data are 
available, such information as period of detention and mean velocity is not 
readily available for most reservoirs. 

In Brune’s study a thorough search was made for all reliable records of 
reservoir trap efficiency, canvassing all agencies of the Federal Inter-Agency 
River Basin Committee known tobe interested in the subject. Information was 
also gathered on capacity, annual inflow, shape of reservoir basin, type of out- 
lets and method of operation, observed gravity underflows or density” cur- 
rents, and any other pertinent data. 

Some forty-four records were gathered and are summarized in Table 1, in- 
cluded herein. Of these, forty are for normal ponded reservoirs, with infor- 
mation on the effect of sluicing and venting operations on three. Two records 
are for desilting basins, and two for semi-dry reservoirs. 


TRAP EFFICIENCY 


DETENTION TIME IN DAYS 


FIG, 2.—TRAP EFFICIENCY AS RELATED TO DETENTION TIME, IMPERIAL 
DAM RESERVOIR, YUMA, ARIZONA 


As pointed out by Brune, it is possible to study trap efficiency of reservoirs 
by a number of different methods, and he found it necessary to use nearly all 
of them. Those who are interested in the details are referred to Brune’s paper. 

In taking suspended-load measurements downstream from a reservoir it is 
necessary totake themas near the dam as possible, because downstream from 
reservoirs streams tend to pick up a new sediment load quickly by scouring 
and degrading their channels. For example, as C. B. Brown has pointed out,12 
the average sediment concentration in the Colorado River below Hoover Dam 
as of 1947 increases from 0.006% 12.8 miles below the dam to 0.091% 98.6 
miles below the dam, or 15 times. Sucha change can cause grave errors in 
computations of trap efficiency if suspended-load measurements far down- 
stream from the dam are used. : 

In computing trap efficiency it is also necessary to take dredging into ac- 
count. Since dredged areas in a reservoir usually fill with sediment again in ~ 


12 “Sediment Transportation,” by C. B. Brown, Chai “ i 
- ziti 3s ’ pter XII of “Engineering Hydrau- 
lics,” by Hunter Rouse, John Wiley and Sons, New York, N. Y., Aug, 1850b. Rants 
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a short time, dredging tends to increase the amount of sediment trapped in a 
reservoir. 

In Table 1, reservoirs and ponds of allsizes were used, ranging froma farm 
pond with 0.038 sq mile drainage area to Imperial Dam Reservoir, draining 
184,600 sq miles. The chief criteria in selecting reliable records were accu- 
racy with which trap efficiency could be determined and the availability of in- 
formation on the various factors which affect trap efficiency. As Gottschalk 
stated,8 the laws of sediment deposition are the same for all types of reser- 
voirs, including stock ponds. Likewise, the trap efficiency of reservoirs is af- 
fected by the same factors, regardless of the size of the reservoir. 

It has long been recognized1l2 that capacity-inflow (C/I) ratio would be a 
more accurate indexof trap efficiency than the capacity-watershed (C/W) ratio 
which has heretofore been so widely used. The C/W ratio must be used within 
definite hydrologic regions, and notas an indexof comparisonover the country 
as a whole, for, as is obvious, with the same C/W ratio the trap efficiency will 
increase as the runoff per unit area decreases. 

Table 1 of Brune’s paper is reproduced herein and the records have been 
plotted as shown in Fig. 3. Brune found that correlation for normal ponded 
reservoirs (conventional reservoirs as distinguished from desilting basins and 
dry reservoirs) operated without any special efforts at sluicing sediment, is 
much better using the C/I ratio than using the C/W ratio. 

Records 32-37ain Table 1, and in Fig. 3, represent reservoirs suchas Fort 
Peck Reservoir and Lake Mead which have C/I ratios greater than 1.0. In 
other words, these reservoirs have a capacity larger than the annual runoff 
from their watersheds. When the effects of evaporation and seepage are also 
considered, it is obvious that water will rarely be discharged over the spill- 
ways of such reservoirs, and that trap efficiency must therefore be close to 
100%. 

The C/I ratio thus provides a means of differentiating between the “hold- 
over storage reservoirs” and “seasonal storage reservoirs” described by 
Thomas Maddock, Jr.13 Reservoirs with a C/I ratio of 1.0 or less may be 
classed as seasonal storage reservoirs, and those with a ratio greater than 
1.0 as holdover storage reservoirs. 

Reservoirs may have very different C/W ratios and annual inflows, yet have 
the same C/I ratio and the same trap efficiency, other factors being equal. 
For instance, a reservoir in the humid part of the country may have a C/W 
ratio of 100, an annual inflow of 1,000 acre ft per sq-mile, and a C/I ratio of 
0.1. Another reservoir in an arid region may have a C/W ratio of 10, an an- 
nual inflow of 100 acre ft per sq mile, and the same C/I ratio of 0.1. Both 
reservoirs, assuming they are of the normal ponded type, will have roughly the 
same trap efficiency (around 87%, from Fig. 3). 

Theoretically, a reservoir can have a trap efficiency of zero only when its 
capacity is reduced to zero or the runoff approaches an infinite quantity. Also, 
a trap efficiency of 100% is theoretically possible only when the reservoir can 
contain the inflow long enough to allow all sediment to settle out. An inspection 
of Table 1 and Fig. 3 shows that there are several instances on record of trap 
efficiencies of zero and 100%. 


SS a a ee eee 
13 “Reservoir Problems with Respect to Sedimentation,” by Thomas Maddock, Jr., 


pp. 9-14, Proc., Federal Inter-Agency Sedimentation Conference, Washington, D. C., 
Jan. 1948. 
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Williams Reservoir, a power reservoir near Williams, Inc., had a trap ef- 
ficiency of less than zero during the period of measurement between 1930 and 
1939, thus revealing that scouring actually had occurred. Three instances of 
trap efficiency equal to 100% were measured at H. Lage Pond, Fort Peck Reser- 
voir, and Denison Reservoir. Details on these records may be foundin Brune’s 

aper. 
z "Table 1 and Fig. 3 show the effects of sluicing and venting. At Arrowrock 
Reservoir, Idaho, sluicing is estimated to have reduced the trap efficiency 
from 93.0% to 90.3%, and to have increased the sediment outflow from the reser- 
voir by 39%. Likewise, at Conchas Dam, N. M., sluicing operations decreased 
the estimated trap efficiency from 97.3% to 95.8%, and increased the sediment 
outflow from the reservoir by 56%. Those who are interested in the details of 
these and other sluicing and venting operations are referred to the original 

aper. 
As pointed out by Brune, it appears that sluicing and venting operations have 
much better chances of accomplishing their purpose when they are timed to 
intercept gravity underflows as they reach the dam. If gravity underflows are 
not present, such operations should at least be timed to meet the higher sedi- 
ment concentrations brought in by flood flows. It is likely that accurate timing 
to intercept gravity underflows can treble or quadruple the amount of sediment 
lost from the reservoir. 

All of the foregoing pertains to normal ponded reservoirs. A word is now 
in order about desilting basins and dry reservoirs as studied by Brune. The 
trap efficiency of desilting basins is governed largely by their shape. Hadley 
Creek New Desilting Basin, Ill., (Table 1 and Fig. 3) covers 280 acres but has 
an average effective depth of less than 4 ft. A normal ponded reservoir with 
the same C/I ratio would have a trap efficiency of about 85%. This basin, how- 
ever, traps 98.8% of the sediment reaching it, largely because of its shallow- 
ness and the small interval of time required for the sediment to reach the 
bottom. Other factors, such as flocculation, the coarse nature of the sediment, 
and the location of the basin at a point where the stream gradient decreases 
abruptly, have their effect on increasing the trap efficiency. 

An even more efficient type of desilting basin is exemplified by the All- 
American Canal Desilting Basins, Ariz. The method of operation, in addition 
to the factors previously mentioned, causes these basins to have a high trap 
efficiency. A normal ponded reservoir with a C/I ratio as low as these basins 
could be expected tohave a very low trap efficiency, yet the basins trap 91.7% 
of the sediment reaching them. 

Thus, according to Brune, desilting basins may have trap efficiencies above 
90% in a much lower range of C/I ratios than for normal ponded reservoirs. 

Although no information is available on the trap efficiency of truly dry reser- 
voirs, two reservoirs which have been operated in a semi-dry condition pro- 
vide indicative data. 

The Senecaville Reservoir, Ohio, a flood-control reservoir built by the Corps 
of Engineers, and the John Martin Reservoir, Colo., are examples of reservoirs 
being operated in a semi-dry condition. The studies on the Senecaville Reser- 
voir revealed the trap efficiency to be 48.4% when operated with the gates wide 
open and 94.3% with the gates closed. The trap efficiency of the John Martin 
Reservoir was 62.2%, an unusually low figure for a reservoir with such a high 
C/I ratio (1.71). | 

Thus it is apparent that even holdover reservoirs like John Martin, if oper- 
ated so as to allow large flows of water to pass through the dam, may have 


age in the detention pool at a contr 
is determined from consideration o 
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trap efficiencies in the 60% range rather than about 90% as would be expected 


with normal operation. Truly dry reservoirs, such as the Miami Conservancy 
District reservoirs, undoubtedly have much lower trap efficiencies, probably 
in the 10% to 40% range, depending chiefly upon the C/I ratio. 

The more important conclusions reached in the study by Brune have been 
Summarized in the foregoing portions of the paper. However, those who are 
interested in more details are referred to the original paper. There are nu- 
merous references listed at the end of his paper. 


UPSTREAM FLOODWATER RETARDING STRUCTURES 


As pointed out previously, the only known published information available 
on trap efficiency is that by Brune, which has been summarized. Since much 
of this information applied to storage reservoirs, it is not Known whether it 
applies also to the type of floodwater retarding structures (particularly with 
dry pools) being constructed by the Soil Conservation Service inthe watershed 
protection program. For this and other reasons, the trap efficiency studies to 
obtain additional information for design were initiated by the Soil Conservation 
Service under authorization of Public Law 566 in cooperation with the U. S. 
Geological Survey. The studies are being executed by the Geological Survey 
with funds transferred by the Soil Conservation Service. 

The effectiveness of floodwater retarding structures in trapping sediment 
from storm runoff has been observed only in a general way. The objective of 
the trap efficiency studies is to provide planning data for the design of such 
structures in the range of watershed conditions exhibited in the United States, 
to delineate the principal factors that influence trap efficiency, and to study 
the sediment yield and characteristics of small watersheds. The selection of 
individual projects has been made jointly by representatives of the Geological 
Survey and the Soil Conservation Service and has been guided by the relative 
need for design data, by the availability of supplementary hydrologic data and 
by the diversity of watershed characteristics. 

Before proceeding further with the discussion of the cooperative trap- 
efficiency studies it might be well to discuss briefly the nature and general 
design of an upstream floodwater retarding structure. Fig. 4 is a section of a 
typical structure. The structures are rolled earth embankments with earthen 
vegetated spillways. They are designed without conservation pools. Their de- 
sign does provide storage space for sediment, the needed space being deter- 
mined by a sediment source survey of the watershed of each structure. In the 
design for this storage, sediment rates based on present conditions of the water - 
shed for ten years and on treated conditions for 40 yr are used where the 
structure is designed to operate at full effectiveness for a 50-yr period. De- 
pending on the character of incoming sediment, from 10% to 30% of the sedi- 
ment storage is allocated in the flood pool, the remainder being provided be- 
low the inlet of the principal spillway in the sediment pool. 

The principal spillway of a floodwater retarding structure usually consists 
of a reinforcedconcrete drop inlet with a heavily reinforced concrete pressure 
pipe, having a minimum diameter of 18 in., laid on a concrete cradle through 
the embankment. 


The principal spillway automatically releases water from temporary stor- 
olled release rate. Its discharge capacity 


f the various physical and economic factors 
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associated with the planning of a single dam or a system of retarding dams. 
Variations in discharge capacity influence detention storage time, and corre- 
latively the trap efficiency. 

A sluice gate may be constructed in the inlet of the principal spillway to 
permit the release of water stored below the crest pending the accumulation 
of sediment. This is essential in some areas to satisfy downstream water 
rights. 

The capacity of the detention pool is such, that flow through the emergency 
spillway occurs not more often than once in 25 yr for any structure. Many 
such structures are designed with sufficient detention storage to accommodate 
a flood of 100-yr or more estimated frequency before the emergency spillway 
is brought into operation. 

The sediment and flood pool volumes vary with the size of drainage area, 
runoff characteristics, sediment yield and other factors. In general, however, 
the sediment pool storage is less than 25% of the total available storage behind 
the structure. In the Southwest area the average volume of sediment space is 
209 acre ft as compared to 1,069 acre ft of total storage, an average for 146 
structures. The character of sediment, character of native water, and detention 
storage time appear tobe the most significant parameters influencing trap ef- 
ficiency. 


Reinforced Concrete 
Gantilever Section 


FIG. 4.-TYPICAL SECTION OF A FLOODWATER RETARDING STRUCTURE 


So much for the general design of the floodwater retarding structures as it 
might affect trap efficiency. 

The operating procedures are adjusted to the particular characteristics of 
each structure. Twelve structures in eleven states were selected as repre- 
sentative for trap-efficiency studies.14 Information on the operation of the 
structures as of June 30, 1956 covers only the first ten structures. The last 
two structures were added after that date. 

The station name, location, date established, area of watershed, capacity 
of the detention pool (from sediment pool elevation to the earth spillway ele- 
vation), and the detention time of each project are shown in the report in the 
order of their installation. The procedures at each project location are de- . 
signed to provide: (1) qualitative data on the character of water and sediment 
flowing into the reservoir; (2) quantitative data on the amount and properties 
of sediment retained in the pool area; and (3) quantitative data on the amount, 
rate, and properties of sediment discharged from the structure. Occasional — 
samples of inflowing storm runoff, in the order of about 10 per yr, are collect- 


14 “Progress Report, Investigation of Trap Efficiency of Floodwater Retarding Struc- 
ores July 1955 to June 1956,” U.S. Geological Survey to Soil Conservation Service, 
Feb. 11, 1957. 
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@d to describe the concentration and size distribution of inflowing sediment and 
pertinent chemical properties of the inflowing water. The Soil Conservation 
Service has in each instance the responsibility of making periodic pool sedi- 
“mentation surveys to determine the amount, distribution and characteristics 
of the deposits. The frequency of pool surveys is determined by the rate of 
accumulation and the occurrence of unusual flood events. The amount of sedi- 
ment discharged from the structure is measured by sampling at a frequency 
sufficient to delineate a continuous record of sediment concentration with time. 
Sampling frequency ranges from hourly intervals during critical periods of 
storm discharge to weekly intervals during periods of base flow discharge. 
The size distribution and flocculation characteristics of the sediment, and perti- 
nent chemical properties of the water, are determined about ten times per 
year. 

Supplemental data relating to precipitation, runoff, inflow-outflow relation- 
Ships, physical and cultural characteristics of the watershed and physical 
characteristics of the structure are available in each project watershed from 
concurrent investigations by the Soil Conservation Service and Geological 
Survey. 

Operation of five of the projects was begun in the interval between January 
and June 1956. Study of the two reservoirs in Nebraska is unique in that struc- 
ture No. 2 is within the watershed of structure No.1. At the present time only 
brief information covering the results of the studies is available. The estab- 
lishment and testing of relationships used in computing water and sediment 
discharge customarily cause more lag between observation and computation of 
data for the initial period than for subsequent periods of operation. 

Since results of measurements as of June 30, 1956, the latest available in- 
formation on some of the structures, are meager due to the shortness of the 
period and the few runoffevents since initiation of the studies, they are not be- 
ing included. 

However, the numerous storm events that occurred during March, April and 
May 1957 in Texas, Oklahoma and Arkansas provided considerable information 
and data relative to sediment inflow and outflow for the floodwater retarding 
structures on which trap-efficiency studies are being made.15 These data are 
being analyzed and tabulated by the Geological Survey. 

Although sufficient time has not elapsed and the storm events above the 
structures have not been numerous enough to give suitable data on which to de- 
sign, it appears that a goodstart has been made to obtain reliable and accurate 
data needed for evaluation of trap efficiency for structures of this type. Until 
such time as the data does become available the Soil Conservation Service, 
among others, will of necessity have to adapt the presently available information 
for other types of reservoirs to the type of structure being considered. In this 
respect, indications are that, for the floodwater retarding structures of the 
type being designed by the Soil Conservation Service and having sediment pools 
such as those constructed in the Southwest area, 85% to 95% of the incoming 
sediment is trapped. This appears to be a reasonably conservative estimate, 
but will be revised as more accurate data become available. 


Conservation Service in the woe Gulf 
National ASCE Conven- 

States,” by G. W. Renfro and C. M. Moore, Paper presented at Nat: : 

tion, Eouraulids Division, Jackson, Miss., Feb. 20, 1957, and published in the Journal 

of the Hydraulics Division, ASCE, Oct. 1958, Paper No. 1806. 
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RESERVOIRS, DEBRIS BASINS AND DEBRIS DAMS 


Before reviewing and summarizing the work being done by the Los Angeles 
County Flood Control District in the control of debris-laden streams, it might 
be helpful to know something of the background of the debris basin program in 
the Los Angeles area. The following information is taken from a letter report 
prepared in December 1952 by a member of the District’s staff on the subject 
of “Control of Debris Laden Streams.” 

Southern California has experienced severalsevere floods in the past which 
have produced heavy debris flows from the mountains. One notable example is 
the flood resulting from the storm of December 31, 1933 - January 1, 1934. 

A brush fire in November 1933 burned over an area of 4,500 acres in the 
mountains north of La Crescenta and Montrose about 15 miles north of Los 
Angeles. The mountain drainage in this area has a very rugged topography 
rising some 2,500 ft ina horizontal distance of 1-3/4 miles. The average 
“Slope of the debris cones at the foot of the mountains is between 8% and 10%. 
Individual watersheds in these mountains range from 0.2 to 1.6 sq miles in 
area. 

During the storm mentioned, more than 12 in. of rain fell. This amount of 
rain falling on steep, recently burned-over areas produced great quantities of 
debris. It is estimated that as much as 87,000 cu yd of debris persq mile was 
produced from the 4,500 acres during this storm. Although velocities of mud 
flows were sufficient to move many large boulders, 10 tons or more in weight, 
the bulk of the debris produced was probably less than 1/4 in. in size. 

While a heavy rainfall on a recently denuded surface will produce great 
quantities of debris,a major storm canproduce similar quantities from an un- 
burned saturated area. For example, the storm of March 2, 1938 fell on fairly 
well saturated ground and produced about 72,000 cu yd of debris per sq mile 
from the same area. One area of 0.27 sq miles produced debris at a rate of 
120,000 cu yd per sq mile. 

At the time of the flood of January 1, 1934 a gravel pit existed near the 
mouth of Haines Canyon. This pit had been provided with a temporary spill- 
way to investigate its action as a debris basin,collecting the debris and al- 
lowing desilted water to pass on down the channel. The debris damage down- 
stream from this pit was not as severe as in the areas not Similarly protected. 

With this example before it, the Los Angeles County Flood Control District 
in 1934 embarked on a program of debris basin construction. Later the U. S. 
Army Engineers undertook the construction of additional basins in the Montrose- 
La Crescenta area. 

Essentially, this system consists of basins located at the mouths of potenti- 
ally hazardous canyons andof concrete-lined lateral channels leading from the 
basins to a major collector channel. The lateral channels being located on 
debris cones, are naturally on steep grades and consequently carry the flow at 
high (super-critical) velocities. Careful attention must therefore be given to 
super-elevation at curves and to oscillation of the water prism thereby creat- 
ed. Rectangular reinforced concrete channels have been found to be the most 
satisfactory for these conditions of flow. q 

As mentioned, the New Year’s storm of January 1, 1934, produced about 
87,000 cu yd of debris per sq mile of watershed. This led to the adoption of a 
design capacity value of 100,000 cu yd per sq mile for debris basins in the 


Montrose-La Crescenta area. This has been the accepted design value when 
a permanent installation is considered. 
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However, since debris production is dependent upon a number of variables 
such as rainfall intensity, slope of watershed, vegetation and soil conditions, 
it has been the experience of the District that the rate of debris production per 
square mile for areas other than the Montrose-La Crescenta area can range 
between 40,000 and 10,000 cu yds per sq mile. In view of this wide variation, 
some means of estimating previous amounts of debris flows from a given 
watershed provides the most satisfactory means of establishing production 
rates. At present there is no completely reliable means of using known de- 
bris production rates for a given watershed in determining the production rates 
for a different watershed, except by a general comparison of topographic fea- 
tures and rainfall characteristics. 

The basins themselves have been of two types, (1) excavated, and (2) non- 
excavated. 

The excavated basin is obtained by digging a pit of the required capacity and 
using the material removed to form a dam of suitable curvature around the 
downstream partof the basin. In this dam an adequate concrete spillway struc- 
ture 6 ft to 10 ft in height is constructed, generally in line with and having its 
foundation at the existing streambed. Several uncontrolled discharge ports 
are generally incorporated in the spillway to allow for some release during 
storms, and partial post-storm drainage. It is essential that an intake struc- 
ture be provided from the normal streambed to the bottom of the excavated pit. 
Otherwise, channelerosion will occur at the inlet tothe basinand progress up- 
stream. A line connecting the lips of the inletand spillway structures is gener- 

_ally at a slope that is about one-half to two thirds of the natural slope of the 
stream channel. These basins are particularly appropriate where complete 
desilting, at least up to the ports, is desirable. 

The unexcavated debris basinis similar to a small storage reservoir, with 
a few variations. The dam is generally a compacted earthfill with heavy rock 
armoror reinforced concrete on its upstream slope. For purposes of estimat- 
ing the storage capacities, the spillway lip is at an elevation that will create 
the desired debris storage capacity. The debris slope is assumed to be about 
one-half to two-thirds of the natural channel slope upstream from the spillway 
lip. No inlet structure is required for an unexcavated basin. A provision is 
generally made to allow for some release during storms. This may be done 
by means of ports through the spillway or some other outlet structure. These 
basins lend themselves to partial desilting, that is, to the retention of debris 
of such size only as would be detrimental to channelimprovements downstream 
from the basins if allowed to be discharged therefrom. This may be accom - 
plished by designing the width and slope of the basin soas to cause the desired 
segregation of debris to take place automatically. se 

It has been the policy of the Los Angeles Flood Control District to limit 
debris accumulation to 25% of the original basin capacity before debris re- 
moval is considered necessary. There is no fixed rule for removal, however, 
for a rather flexible basis has been used in the operation and maintenance of 
the debris control system. The allowable debris accumulation is considered 
in the light of the probable debris production and the condition of the watershed 
relative to vegetation and geological characteristics. 

Debris disposal presents a major problem in rapidly growing southern 

j lated debris has no agricultural 

California. Generally speaking, the accumul wale 

value in spite of the very high percentage of fine material. After a ee i 4 

years of operation of a debris control system with the ees Ae poe eS 
moval, the nearby disposal areas are exhausted and more remote 


82 February, 1960 HY 2 


be found. This makes mechanical debris removal a costly procedure. It would 
be desirable to segregate the bulkof the fine debris and cause it to be discharged 
with the water, as previously discussed. This type of debris removal elimi- 
nates the excavation cost of the accumulated debris and minimizes interference 
with the normal replenishment of sand for the beaches of southern California. 

In general, floating debris does not present a great problem for the debris 
basins. The watersheds are usually quite small and the volume of this type of 
debris is not large. However, until the water surface in the basin rises above 
the ports, there is always the possibility of clogging. Patrolmen have been 
able to keep the ports clear more or less satisfactorily during storms. The 
District has been experimenting with inclined trash racks made from railroad 
rails and spaced roughly one-half of the minimum dimension of the ports. The 
results.have been satisfactory so far, though incomplete. When the water over- 
flows the spillway lip the floating debris is allowed to pass on down the chan- 
nel. At the end of the runoff season when the basin has been allowed to drain, 
the remainder of the floating debris accumulates on the sides and bottom. It 
is then either collected and burned or removed from the basin by trucking. 

At times itis necessary to construct temporary debris basins and channels 
before sufficient funds are available to build permanent structures. Under 
these conditions the design capacity, 100,000 cu yd per sq mile, is sometimes 
reduced to 50,000 cu yd per sq mile, orless. Such a basin might have asphalt- 
lined or wooden inlet and spillway structures, as well as an asphalt-lined dis- 
charge channel. The use of asphalt to replace concrete requires caution, par- 
ticularly where impact due to coarse debris and boulders is concerned. Asphalt 
will not stand up under such punishment. It will, however, resist abrasion due 
to ordinary drag of a bed load of debris, the maximum size of which is of the 
order of three inches. 

Necessarily, any flood protection must be based on a definite plan of con- 
trol. Portions of an area might require storage reservoirs for debris control 
and flood flow regulation while smaller watersheds may needonly a system of 
debris basins for satisfactory protection. 

Since December 1952 the most important modification of the early debris 
basin concept and design has been the elimination of dead storage and of the 
open port outlets which were provided in addition to the spillway in many of 
the debris basins, and the substitution of an outlet tower. This is located at 
the low point in the bottom of the debris basin offset from the line of flow be- 
tween the inlet and the spillway and serves as the entrance to a conduit which 
passes under the dam and discharges in the channel downstream. This outlet 
design accomplishes two things: (1) It permits storm flows which are rela- 
tively free of debris to pass on through the basin with a minimum of ponding, 
thus leaving basin capacity for debris storage, and (2) it serves as a trash 
rack. These towers have been constructed in circular form of 5-ft inside di- 
ameter of reinforced concrete or reinforced gunite. They have a 6-in. wall 
thickness with a series of round perforations four inches in diameter on the 
outside and five inches in diameter on the inside face of the tower. The tower 
is of sufficient height to project at least one foot above theoretical debris slope, 
equal to 75% of the original streambed slope. The outlet pipe normally is a 
36-in. reinforced concrete pipe with a slope of not less than 5% and a capacity 
of at least 150 cfs for the maximum expected head of water in the basin. By 
means of a controlling orifice at the base of the tower, flowunder pressure in 
the outlet is prevented. Other modifications in design have no specific bearing 
on the subject of trap efficiency and need not be discussed here. 
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In presenting information relative to the reservoirs, debris basins and de- 
bris dams of the Los Angeles County Flood Control District, it appeared that 
perhaps the best way would be to follow the format of Table 1 of Brune’s paper.4 
This procedure was followed insofar as the data permitted. Tables entitled 
“Reservoir Debris Production Data,” Table 2, and “Debris Dams and Debris 
Basins—Debris Production Data,” Table 3, are included herein. 

From the data available it was not possible to include three columns which 
were included in Brune’s paper. These were (1) average annual inflow, (2) C/Ib, 
and (3) trap efficiency. This omission may be explained by stating that relative 
to debris basins and the majority of the debris dams, continuous inflow and 
outflow records are not available. In fact, very few records of inflow and out- 
flow are available. Due to the lack of inflow records, the values for C/I> and 
trap efficiency could not be computed. However, it is believed that in nearly 
all cases the trap efficiency, at least of course-grained debris, approaches 
100%. With the exception of five of the basins, namely, Aliso Creek, Scholl, 
Snover, Ward Canyon, and Verdugo, the material is trapped, as will be later ex- 
plained. In the abovementioned basins, the silt has been allowed to accumulate 
to the level of outlets. Much of the fine material and that which is in suspen- 
sion has been allowed to go on through. This has been somewhat in the nature 
of an experiment, but so far appears to be a successful plan of operation for 
these basins, at least during seasons where the watersheds are not burned off 
and no subsequent floods occur. In the case of Scholl and Ward Canyon, flash- 
boards are placed in the outlet notches to regulate to some extent the level of 
debris in the basin. At the other basins, except in the case of sizable floods 
which fill the basins with sediment, all material and most of the inflow is 
trapped during a given storm. If, subsequent to a storm, it is found that less 
than 75% storage capacity remains, the material is removed by shovel and 
trucks to a nearby disposal area. 

The omission of the three-above-mentioned columns in the tabulation for 
flood control dams is due to the fact that although the District does have in- 
flow and outflow records, sediment generally goes through the outlets only 
during periods of controlled sluicing. Under certain other conditions such as 
during flood releases, small quantities of extremely fine silt and sand may pass 
through the valves. Normally flood releases are made only when an appreciable 
head exists on the valve. Other releases are made subsequent to storm peri- 
ods and therefore after most of the suspended material has settled out. As 
stated above, the District believes that trap efficiency for the reservoirs 
generally approaches 100%, even during periods when both inflow and outflow 
take place. As the reservoirs serve the dual purpose of flood controland water 
conservation, water is stored to limits compatible with flood operation. 

In some of the reservoirs where conditions permit, capacity is regained by 
the hauling of some of the silt accumulation by permittees. In general, this 
does not constitute as large a storage reclamation as is gained by sluicing, as 
may be noted from the table. Sluicing, of course, is dependent upon sufficient 

‘residual flows. This is often a problem, as the reservoirs may not be drawn 
down due to conservation commitments until late summer or fall when inflows 


are low. : 
A study has recently been comple 
ing Division, relative to the quality o 
and removed to disposal areas; the latter have 
_ the District exclusively for disposal of debris. 
terial trapped by the debris basins is of interest 
-_jectof trap efficiency, a study was made of grading an 


ted by the District’s Foundation and Test- 
f debris excavated fromthe debris basins 
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LOS ANGELES COUNTY FLOOD CONTROL DISTRICT 


Project Dunsmuir Debris Basin _ 
Hole No. Composite of 5 Pits 
Sample No. Depth 

Loc. Elev. 


Disposal Area 


Comparison of Debris Material with 
80 Pit Material from San Fernando Valley 
—— Rock Producing Companies 


FOUNDATION AND TESTING DIVISION 
GRADING ANALYSIS 
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Percent Rock 
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LOS ANGELES COUNTY FLOOD CONTROL DISTRICT 
FOUNDATION AND TESTING DIVISION 


Project Hails Debris Basin 
Hole No, Composite of 3 Borings 


Sample No, Depth 
Loc. Disposal Area Elev. 
100 200 100 


Comparison of Debris Material with 
90 Pit Material from San Fernando Valley 
— Rock Producing Companies 


Sampled by C.L.&A.L.R. Date 5/57 
Tested by P.P.&C.L. Date 5/57 

45.3 Percent Sand 48.4 

Rock 7.70  F. M. Sand 2.17 _ 


GRADING ANALYSIS Sampled by ek, Date s/s7 
lab. Tested by P.P.&C.L. Date S/57 
Percent Rock 24.0 Percent Sand 64.4 
FLOM. Rock 8.20 F. Me. Sand 2.36 
ASTM SIEVE NUMBER 
50 30 16 8 4 3/8" 3/4" 13" 3" 


(LIMITS PER A.S.T.M. © 125) 


FIG. 6, 
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basins including Cooks, Dunsmuir, Halls, Pickens, Ward, and West Ravine. 


| The grading analyses for the various basins were compared with pit material 
_Trom the San Fernando Valley rock producing companies. Figs. 5 and 6 are 


typical gradation curves for material from the Dunsmuir and Halls debris 


basins. 


The rock content of the debris material found in the various disposal sites 


-yaried from 10% to 45%, with an average of 29%. This is considerably less 


than the average of about 58% rock found in the commercial material. 

In general, the gradation of the rock was considerably finer than that of the 
commercial sources. Excluding the material from West Ravine, Halls, and 
Aliso-Wilbur debris basins, the debris rock compared favorably with that of 
the commercial sources relative to abrasive loss. 

The debris material from the various disposal sites contained between 
48% and 74% sand, with an average of 63%. This is greater than the average 
of 40% sand found in the commercial material. The quality of the debris sand 
compared favorably with the commercial material. 

In order to determine more accurately the degree of trap efficiency of the 
Los Angeles County Flood Control District structures, it will be necessary to 
obtain more data directly related to the subject. These data should include 
continuous inflow and outflow records and measurements of sediment coming 
into and passing through the structures. If such data were available, the trap 
efficiency could be computed on a C/I ratio basis and perhaps a correlation 
could be made for these structures as has been done for normal ponded reser- 
voirs. 
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HYDRAULICS DIVISION 


Proceedings of the American Society of Civil Engineers 


STRATEGIC ASPECTS OF URBAN FLOOD PLAIN OCCUPANCE? 


By Gilbert F. White! 


SYNOPSIS 
The paradox of continuing flood control and rising flood losses is ex- 
plained in part by elements entering into decisions to change urban 
flood-plain occupance, 1936-1957. Curbing the mounting toll calls for 
public action to broaden the range of choice open to managers of flood- 


plain properties in adjusting to flood hazard. 


INTRODUCTION 


A striking feature of national efforts at flood control in the United States is 
that in the more than two decades since a national policy was launched in 1936, 
the mean annual toll of flood losses has continued to rise. While at least 4 bil- 
lion dollars have been spent for engineering works to reduce and control 
floods, the economic losses from occasional flood disaster have mounted. To 
understand this paradox it is necessary to examine what has been happening in 
the human occupance of flood plains. To find out the full import of changes in 
flood-plain occupance it is essential to identify those factors which enter into 
both public and private decisions as to flood-plain use. 

This paper outlines briefly the broad problem of flood-loss reduction in the 
United States, reviews the findings from a recent study of changes in urban oc- 
cupance of flood plains, and then suggests three strategic aspects of such oc- 
cupance that seem likely to affect the results of any further engineering efforts 


at flood control. 


 Note.—Discussion open until July 1, 1960. To extend the closing date one month, a 
written request must be filed with the Executive Secretary, ASCE. This paper is part 
of the copyrighted Journal of the Hydraulics Division, Proceedings of the American So- 
ciety of Civil Engineers, Vol. 86, No. HY 2, February, 1960. 

a Presented at the July, 1959 ASCE Hydr. Conf. in Fort Collins, Col. 

1 Dept. of Geography, The Univ. of Chicago, Chicago, ll. 
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The broad problem of flood-loss reduction is that the rate at which floor 
losses are being eliminated by construction of engineering or land-treatmen 
works is of about the same magnitude as the rate at which new property is be 
ing subjected to damage. Even though the data on which estimates of nationa 
flood- loss are based seemtoo inaccurate to warrantany precise comparison 0 
the two rates, it is clear, both from the aggregated damage statistics and fron 
the record of selected flood plains, that the heavy investments in flood protec. 


THE NATIONWIDE FLOOD 


y PROBLEM 


RESIOUAL DAMAGE AFTER aef, 
AUTHORIZED PROGRAM pee® 
1S COMPLETE C 

eee ee 

~8' ¢' 


= 


WORKS COM — 
PLETED OR IN 
OPERATION 


TO BE PREVENTED 
BY AUTHORIZED 

BUT NOT STARTED 
WORKS 


TO BE PREVENTED BY WORKS 
UNDER CONSTRUCTION 


WOTE: EXCLUSIVE OF FLOOD DAMAGES IN SMALL UPSTREAM 
TRIBUTARIES ESTIMATED BY DEPT. OF AGRICULTURE 
AT #300 MILLION ANNUALLY. 


FIG. 1.—GRAPH PRESENTED BY THE CORPS OF ENGINEERS TO THE HOUSE 
COMMITTEE ON PUBLIC WORKS IN 1957 SHOWING AN ESTIMATE OF 
MEAN ANNUAL FLOOD-DAMAGE POTENTIAL AND OF DAMAGES 
PREVENTED. (FROM 85TH CONGRESS, 1ST SESSION, HOUSE 
COMMITTEE ON PUBLIC WORKS, PRINT NO. 1, 1957). 


tion have effectively curbed the losses in many areas but that new damage po 
tention is being built up at the same time. 

Studies of seventeen selected urban areas having flood problems reveal 
general and persistent encroachment of urban structures upon the flood plain 
during the period 1936-1957, even in areas in which there was net decrease i 
total population. They show certain distinctive patterns of encroachment, an 
indicate highway construction and flood-control works astwo major stimulant 
to growth. The situations at Boulder, Colo. and Denver, Colo., illustrate som 
of the major findings. : 
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_ In the light of that experience it is argued that at least three aspectsor ur- 
ban occupance not ordinarily considered in flood-control plans in the past must 
be taken into account in the future if the tide of rising flood losses is to be 
turned. First, itmust be recognized that engineering works are only one of the 
possible human adjustments to flood hazard. Second, it must be seen that the 
complex of elements entering into decisions as to future occupance of flood 
plains includes many considerations in addition to the traditional one of cost- 
benefit evaluation. Third, the range of choice now permitted property man- 
agers in dealing with the flood hazard is so restricted that radical changes 
must be made in public policy to broaden the range of choice among possible 
adjustments and to assure a full appraisal of each choice. 


TRENDS IN FLOOD LOSSES 


There are two major sources for national estimates of flood losses. The 
Corps of Engineers issues, from time to time, an estimate of the total poten- 
tial flood losses and of the losses that have been prevented or will be pre- 
vented by works constructed or planned by the Corps. A sample graphical 
statement of such an estimate for 1954 is given in Fig. 1. As of that year the 
Corps found a potential mean annual loss of 964 million dollars.2 (All esti- 
mates are adjusted to 1957 price levels.) After all works then authorized had 
been completed there would remain a balance of approximately 210 million dol- 
lars that would not be prevented. These losses are computed separately for 
each of the major hydrographic areas of the country (Fig. 2). 

The Weather Bureau issues the other series of flood-loss statistics, having 
published each year since 1903, an estimate of reported losses for each of its 
districts. The national totals for 1902-1955 are shown in the lower part of the 
bars in the graph in Fig. 3. These totals have been adjusted to 1957 price 
levels, and the full height of the bar shows the totals on a comparable basis. 

No attempt will be made here to appraise the validity and the discrepancies 
among these two sets of estimates and other less comprehensive series. This 
has been done in a separate publication.? Aggregate data alone are unsatis- 
factory and need to be checked in more detailed fashion. It is sufficient here 
to point out one characteristic which the two national estimates have in com- 
mon. Bothsuggesta pronounced upward trend in the size of annual flood losses 
since 1936, when the first national flood control legislation was enacted. The 
estimated mean annual losses in 1936, were about 212 million dollars. The 
Corps estimate of remaining annual losses after allowing for works completed 
was 444 million in 1954, and 700 million in 1959, (at 1959 prices).4 The Weath- 
er Bureau estimates of mean annual losses for 1924-1953 were 25% lower than 


those for 1944-1955. 


2 “House Committee Print,” U. S. Congress, House Committee on Pub, Works, 1957, 


No. 1, 85th Congress, 1st Sess. p. 2. ; , 
3 «Changes te Urban Occupance of Flood Plains in the United States,” Gilbert F. 


White, Wesley C. Calef, James W. Hudson, Harold M. Mayer, John R. Sheaffer, ae 
Donald J. Volk, Univ. of Chicago Geography Research Paper, 1958, No. 57, pp. 1-19. 
_ 4 Hearings on Civil Works Appropriations Bill,” U. S. Congress, House Appropria- 
tions Committee, 1959, Pt. I, p. 18. 
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REASONS FOR MOUNTING FLOOD LOSSES 


One partial explanation for the apparently rising toll of recorded losses is 
that the inflation in dollar values has exaggerated the size of recent estimates. 
As shown in Fig. 3, even when values are adjusted to a common year there 
still remains a clear and sizeable increase. 

A second reason advanced for the apparently mounting losses is improve- 
ment in enumeration and estimating methods. This unquestionably has played 
a role as Weather Bureau procedures have been made more uniform and as 
Corps of Engineers studies, have reached areas not previously covered. Some 
who _ have been studying this problem believe that changes in enumeration meth- 
ods may account for as much as 10% to 15% of the increase. 

Perhaps of greater importance is the difference in number of large andin— 
frequent floods. William G. Hoyt, F., ASCE and Walter B. Langbein, F., ASCE 
have shown a 35% increase_in their “flood index” from the first 25 yr. of the 
century to the second 25 yr.’ There was a remarkable bunching of rare flows 
in the lower Missouri and Northeastern basins during the 1950’s. This hydro- 
logic record may account for as much as 25% of the increase in losses. 

Probably the most important reason for the rising trend in flood losses is 
to be found in the continuing encroachment of human occupance upon flood 
plains. This takes the form of new structures, of changes in the intensity of 
existing structures, and of structures which so reduce the hydraulic efficiency 
of valley sections as to increase the hazard in affected reaches of the stream. 
Although there have been numerous studies of flood-control projects there has 
been no comprehensive investigation of the changes that have been occurring in 
the areas subject to flood. One recent study conducted by University of Chicago 
geographers dealt with urban occupance because it has a high and increasing 
proportion of flood losses, is compact and is more susceptible to change. It 
shows that there are at least 1,020 urban places with a populationof more than 
1,000 that have well defined flood problems. That study also gives a precise 
picture of the changes in seventeen urban areas selected for their diversity in 
valley section, flood frequency and height, size, population growth, and types of 
land use. A few findings from those field examinations will be summarized, 
and then illustrated from the Boulder and Denver situations. 


GROWTH PATTERNS 


The most evident and widespread trend that is to be observed in the urban 
flood plains is one of growth. In every place studied, including several in 
which the total population declined during the 21-yr. period, the number of 
structural units in the flood plain increased.6 The growth rate ranged from 
less than 2% in a stable city like Wheeling, W. Va., to more than 600% ina 
rapidly expanding place such as Dallas, Tex. 


ee eee 
5 “Floods,” by William G. Hoyt, and Walter B. Langbein, Princeton Univ. Press, 1955. 
6 Changes in occupance were measured in “structural units.” A unit was defined for 
different occupance classes as: Residential, A-a single or double-family dwelling, B-a 
multi-family dwelling for 3-6 families, and C-a multi-family dwelling for more than 7 
families; Commercial, Industrial and Transport, separate structures in multiples of 


aunt sq. ft, and open working space in multiples of 25,000 sq. ft; Public, each separate 
uilding. 
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In most places the commercial and industrial structures grew at a more 
-apid rate than residential structures, although in a few of the places such as 
in sectors of Los Angeles, Calif., the residential expansion was substantial. 
Public structures in most places showed substantial growth rates. 

Certain areal patterns of change seem to repeat themselves in urban flood 
dlains (Fig. 4). One typical pattern is found where residential areas in or bor- 
ering on the flood plain enlarge by moving down nearer the channel. A second, 
.§ that in which industries already established in a flood-hazard zone expand 
aterally along the stream bank and at right angles away from it. A third oc- 
zurs where a central business district spreads out into the flood plain, follow- 
ing the major traffic arteries. 

There is no readily discernable relation between rate of population growth 
and rate of change in structures on the flood plain: Too many other factors 


River Flood Plain ou 2 Occupance 
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FIG. 4:—SCHEMATIC MAPS OF THREE STRATEGIC LAND-USE ASSOCIATIONS, 
SHOWING PATTERNS OF GROWTH IN THE FLOOD PLAIN. 


such as urban function and availability of landplay a part. Itis clear, however, 
hat two forms of public action have had a powerful stimulating effect upon in- 
fasion of flood plains. These are the highway program and the flood-control 
rogram. New highway construction in urban areas tends to follow the low 
rradients of stream valleys and the less densely settled sectors of some flood 
lains. It has caused substantial removal of low quality residential structures 
wut it also has been an inducement for commercial and industrial establish- 
ents to follow along into the flood hazard zones. 
i" As might be sare the construction of new flood-protection works fre- 
juently has been the signal for accelerated movement into the flood ene 
Thus, the completion of reliable works along the Trinity River in pee eX. 
ssured development of a large-scale commercial district behind the levees. 
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But there are certain implications of flood-control works that are not as obvi 
ous. The lack of protection is not necessarily a deterrent to urban invasion ¢ 
flood plains: In cities which have had no serious flood in 50 yr. as well as i 
cities that had unprecedented flooding only a few years before, there has bee 
extensive building without any immediate prospect of protection. Moreover, i 
certain valleys, such as in the Tennessee at Chattanooga, the completion c 
river-control works upstream has been followed by further movement into th 
flood-plain where frequency has been reduced but where high flows still ar 
possible. Although most Federal flood-control works are built to protect agains 
a project flood and conceivably will one day, however infrequently, be exceede 
by a larger flow, there is a universal disposition to believe that the rare flo 
will never come. This means that the number of situations in which a catas 
trophic disaster may follow one of those rare flows is increasingas new levee 
are completed. 

These generalizations are illustrated in the near-by cities of Boulder an 
Denver. 


FLOOD-PLAIN OCCUPANCE IN BOULDER 


There has been no major flood in Boulder since 1894, when waters from a 
intense rainfall in the Rocky Mountain Front area of Boulder Creek inundate 
the areas to the north of the Creek as far as Water Street. (Fig. 5) In 195 
Congress authorized a combination levee, channel improvement, and bridg 
reconstruction project to protect sectors of the plain that had been built u 
over the years. The authorized project was to accommodate a flow equal t 
that of 1894, or about 12,000 cfs. A flow of as much as 15,000 cfs would nc 
be impossible. The City decided not to take part in the authorized projec 
which would have required a local contribution of $1,192,000 out of a total es 
timated cost of $1,707,000. So the situation stands today with no protectio 
works, and with some local interests pushing for a Department of Agricultur 
watershed survey as a possible line to an alternative solution. 

Since 1938, the area which was flooded has been the arena for vigorous ne 
construction. Single-family residential structures have increased 34%, com 
merical structures 26%, and public structures 300%. The University of Colo 
rado has built twenty eight new multiple-unit housing structures in the floo 
plain. The City has constructed its municipal building between the channe 
and the line of the proposed levee. The largest bank has its new buildin 
astride the proposed levee line. One modern hotel is under construction | 
short distance below the canyon mouth, and another along the creek betwee: 
24th and 28th streets. 

Some managers of property in the flood plain believe that the area neve 
will be flooded again. Others are taking a calculated risk, arguing that th 
benefits of daily use will more than offset the lossesfrom an occasional over 
flow, and that another flood may not come in their period of occupance. The 
are correct in thinking that Boulder has not recently added to the national bil 
of flood losses. But when a flow equivalent to the one of 1894 does come dow 
the channel—perhaps next year, perhaps not in our lifetime—the city will be o: 
the front page of every newspaper in the country as a place where propert 
losses exceeded several millions of dollars in a few hours. Hopefully, th 
pounding creek waters, laden with boulders and debris, may not cause loss a 
life. If the experience elsewhere applies here, that unpredictable but certai 
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FIG, 5.—-MAP OF BOULDER, COLO. SHOWING APPROXIMATE FLOOD PLAIN 
AND LOCATION OF SOME NEW DEVELOPMENTS. 
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flood will be followed by organization of an emergency citizens committee to 
obtain Federal flood control for Boulder. 


FLOOD-PLAIN OCCUPANCE IN DENVER 


Denver has three distinct flood situations, one of which has been largely 
solved, one of which daily becomes more acute, and one of which is in course 
of halting but imaginative solution. The Cherry Creek flood plain (Fig. 6) was 
the scene of a damaging overflow in 1933 that covered large sections of resi- 
dential and commercial occupance in the city. This area now is protected by 
a Federal detention dam capable of controlling a storm runoff of more than 
600 cfs per sq. mile. There has been argument that such protection is unduly 
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FIG, 6.—MAP OF DENVER, COLO. SHOWING APPROXIMATE FLOOD PLAIN 
AND LOCATION OF MAJOR TRANSPORTATION LINES. 


expensive, but, regardless of cost, that sector of Denver seems to be off the 
rolls for flood-loss reporting. 

Along the main stem of the South Platte River, both above and below the 
mouth of Cherry Creek, there is a broad plain that received partial protectior 
from WPA works and that has not experienced flood flows exceeding the chan- 
nel capacity of 13,000 cfs in recent years. The Congress, in 1950, authorizec 
a reservoir and levee protection project to control flows of 85,500 cfs for this 
area. The dam was never started, and, indeed, could not be started today with- 
out dislodging a unit of the Martin Aircraft Corporation and a residential sub- 
division, both of which have since invaded the proposed reservoir area. Mean- 
while, in the South Platte flood plain commercial structures have increasec 
26% since 1936, industrial structures 40% to 50%, and residential structures 
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slightly, notwithstanding extensive demolition for a new expressway inthe val- 
ley. The damage potential unquestionably has increased at a rapid rate in this 
lower valley area. 

A little progress has been made in the small tributary gulches where there 

s been a persistent tendency for residential and commercial occupance to 
press toward the dry stream beds, thus reducing the capacity of the channels 
and exposing more property to flood loss. Denver is one of the few cities that 
has begun a program of zoning, building regulation, and land acquisition to 
control further encroachment on such dry gulches, and a start has been made 
at metropolitan area planning for regulation across political boundaries. ! 

While the Cherry Creek area of flood loss has been eliminated by engineer- 
ing works, the South Platte damage potential has been growing under the stimu- 
lation of highway and partial flood-control construction. A modest startis un- 
derway to curb further encroachment on the tributary gulches. 


RANGE OF CHOICE IN ADJUSTMENT TO FLOOD HAZARD 


The traditional approach to flood losses in the United States has been to 
either control the flood flows by engineering works or to bear the losses, let- 
ting the public come to the aid of the less prosperous flood sufferers. This is 
shown by the heavy public emphasis upon investment in Federal flood-control 
works and by the Red Cross and other emergency aid programs. 

There are, of course, other possible adjustments which individual property 
managers can make. The whole range of choice may be listed as follows: 

Bear the Losses.—A good many households and firms bear the losses when 
they occur, and a few actually make financial plans to do so. 

Emergency Evacuation and Rescheduling.—Property may be removed from 
the reach of floods, and property movement and production operations may be 
rescheduled so as to avoid losses through interruption. This requires a rela- 
tively accurate system of flood forecasting and a plan for emergency action 
when the critical forecast is received. 

Prevent Flows.—In some areas it is practicable to prevent certain flows, 
oarticularly the more frequent flows, by land treatment and associated mea- 
sures. The experience with such measures is only beginning to build up. 

Elevate Land.—By land fill it is practicable in some situations to raise 
oroperty above the level of flood waters. The effect upon channel capacity 
yaries according to location with respect to flowage and pondage areas on the 
‘lood plain. 

Control Flows.—This is the conventional engineering solution, involving 
shannel improvement, levees, cut-offs, and storage or detention dams. 

Change Structures.—In advance of a flood warning, structures may be al- 
ered so as to prevent or reduce flood losses when the water rises. These 
shanges may include such alterations as packing of machinery, bricking in of 
Ow openings, cut-off valves on sewers, and rearrangement of electrical cir- 
uits. 

Change Land Use.—The use of flood-plain land may be changed so as to in- 
roduce a use that is less susceptible to flood loss. This may range from the 


4 7 “Report for Inter-county Regional Planning Commission on Storm Drainage in the 
Denver Metropolitan Area,” by Dale Rea, 1957. 
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transfer of an entire town from a riverine to an upland site, to public acquisi- 
tion of flood-hazard areas for recreational or parking purposes. It may be 
guided by public regulation, including zoning, building ordinances, subdivision 
regulation, and land acquisition. 

Insure.—Although insurance against flood losses is not generally available 
in the United States, there is an inactive program for Federal-state subsidized 
insurance under the Flood Insurance Act of 1956, andthere are a few instances 
of coverage by private companies where special structural adjustments have 
been made. 

Public Relief.—Through the medium of the Red Cross or by direct grants 
and loans under Public Roads, Civil Defense, Corps of Engineers and small 
business programs, public aid is given to individuals and local governments 
that suffer heavily from floods. 

There are concrete examples of each type of adjustment, and there are 
strong arguments that can be made for and against each type of adjustment, 
depending upon the local situation. To understand how and why particular ad- 
justments are made in a given situation it is necessary to look to the various 
elements that enter into decisions as to management of urban flood-plain prop- 
erties. 


ELEMENTS IN DECISIONS AS TO FLOOD-PLAIN USE 


As with most other resource-management situations, at least seven con- 
siderations seem to enter into decisions as to flood-plain adjustments. It may 
be helpful to examine how each one figures in the urban situations which have 
been studied. 

Estimating the Resource.—There is widespread ignorance of the flood hazard 
and a tendency to minimize it. Many people building or buying in flood plains 
are unaware of the precise hazard they are running, or grossly misinterpret 
the technical estimates. This applies even in places where there have beer 
public plans for flood protection. A man says he need not worry about floods 
because a 200-yr. flood occurred the year before. A Federal housing insur- 
ance office says it does not insure mortgages for new buildings in flood plains 
but has no map showing where the flood-hazard areas are located. 

Discounting Future Benefits and Costs.—It is common in Federal flood pre- 
vention and flood-control studies to discount future benefits and costs from the 
proposed works. These cost-benefit ratios have significance in Congressional 
decisions so long as they are below unity. There isa tendency to group all 
programs with a ratio of more than unity together and to pay little attention tc 
the means by which the ratios were calculated. So far as private planning o: 
flood-plain adjustments are concerned there is little evidence of such discount- 
ing procedure. A few managers of large industrial units make cost-benefi 
calculations but most managers do not. 

Harmonizing Two or More Uses.—In the Federal plans there is a disposi- 
tion to harmonize flood-control plans with other water uses such as irrigatiorx 
and navigation. At the level of local governments and private managers suck 
multiple-uses tend to be ignored. Few attempts are made to prevent harmfw 
encroachments, and there is little attention to combining levees with highways. 
or to harmonizing a park development with a floodway improvement. One o: 
the hopeful moves in this direction is under planning supported by the Urba1 
Renewal Administration. 
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Projecting Future Demand.—Typically, the Federal agencies assume that 
nere will be little change in the demand for flood-plain land while private 
aanagers tend to overestimate the effect of protection and to rush in at the hint 
r prospect of some control work. 
Projecting Technological Change.— Private managers tend to overestimate 
ne physical benefits from any kind of prevention or protection work. The 
ompletion of a single-purpose power dam upstream, (and even downstream), 
ne beginning of a watershed improvement, or the clearing of an upstream 
hannel is taken as assurance that floods will be abated in the future. Thus, a 
ttle protection work may encourage a great deal of channel encroachment. 

Integrating Regional Uses.—While the Federal agencies try to assess the 
pnsequences of a flood-control program upon other water uses inthe same or 
djoining basins, there is little attention to this on the part of local and private 
lanagers, and metropolitan area plans for storm water disposal are rare. 
ven less common are efforts to assess the possible use of flood plains in 
erving regional requirements for land and for riverine facilities. 

Setting Social Guides.—We come now to one element which affects all the 
est. Society, through its public agencies, definitely restricts some kinds of 
ecisions by property managers and clearly encourages some other kinds of 
ction. Federal agencies tend to guide the occupance of flood plains by pro- 
iding prevention and protection plans, by giving public relief, by issuing storm 
nd flood warnings, and by offering information on flood occurrence. Seven 
tates attempt to regulate channel encroachment. A few cities and counties 
xercise some regulation over flood-plain use. One Federal agency cooperates 
ith state planning agencies in assisting with local plans for flood-plain use.8 


BROADENING THE RANGE OF CHOICE 


In the present circumstances the social guides in the form of information, 
2gulation, and investment affecting flood-plain occupance tend to encourage 
irther encroachment upon the flood plains at the same time that they lead to 
2avier Federal expenditures for flood control. It is not entirely whimsical to 
1y that the national situation is somewhat like a local situation where investi- 
ition showed forty new houses in a flooded area could not be economically 
rotected by levee works. It was noted, however, that if twenty additional 
yuses were to be constructed protection then would be feasible. Giyen the 
tisting policies one may expect that the additional twenty houses will in time 
> built and that Federal largess will follow in due course. 

We know that the managers of flood-plain properties often underestimate the 
zard, fail to discount future benefits and costs, overestimate the demand for 
eir land, overestimate the effects of engineering and land-treatment works, 
id ignore possible combinations of flood-loss reduction with other local im- 
‘ovements or with regional plans for land and water use. We also know that 
; a practical matter most managers who face aflood threat have a choice be- 
een only two alternatives. They can bear the losses, or they can press for 
Federal project of some sort. They do not receive technical advice as to the 
sibilities of emergency evacuation, land elevation, structural changes, or 


i 
8 “Regulating Flood-Pain Development,” by Francis C. Murphy, Univ. of Chicago 
sorgraphy Research Paper, 1958, No. 56. 
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land-use changes. They do not have insurance readily available. If they suffe 
unduly heavy losses they have the prospect of public aid. If they are energeti 
enough they may receive public protection. With few exceptions they do n¢ 
know the precise character of the flood hazard, and they are not rubjectto pub 
lic curbs against any further encroachment upon the stream channel. 

If we are to break out of the present situation in which floodlosses promis 
tokeep pace with flood protection for at least several decades it seems essen 
tial to broaden the range of choice. Each manager—public or private—shoul 
be given the opportunity to choose among the whole range of possible adjust 
ments, with the public agencies strictly limiting those choices which woul 
cause damage to others or to the public safety. This calls for a new and fres 
approach to the problem of flood loss reduction inthe United States. An impor 
tant beginning in that direction was made at the Conference called by the Coun 
cil of State Governments in December, 1958.2 The new approach does not re 
quire abandonment of existing public programs. It requires that they kh 
supplemented by new or expanded efforts to broaden the choice open to all wk 
face flood hazard. 

The chief directions in which it now seems necessary to move are these: 


1. Publication and wide distribution of flood hazard maps and reports for a 
important flood hazard areas, to be made available to both property manage1 
and to public lending and construction agencies. 


2. Technical advice to property managers as to the means and costs of re 
ducing flood losses by emergency measures, by structural changes, and k 
land-use changes. 


3. A cooperative Federal-state program of flood insurance in which pre 
miums are in proportion to hazard. 


4. State regulation ofany further channel encroachments which would caus 
damage to others or to the public safety, and Federal requirement that add: 
tional Federal expenditures for flood control be contingent upon such regulz 
tion. 


5. An improved system for issuing flood warnings and for getting them int 
the hands of the property manager concerned. 


6. Technical assistance to local governments in drawing up plans for re 
ducing flood losses by whatever combination of means may seem most feat 
ible. 


EE, 


9 “A New Attack on Flood Losses,” by Gilbert F. White, State Government, Sprit 
1959. 
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APPLICATION OF SNOW HYDROLOGY TO THE COLUMBIA BASIN2 


Closure by O. A. Johnson and P. B. Boyer 


O. A. JOHNSON, F. ASCE! and P. B. BOYER.2—Mr. Miller has presented 
1 well-written, thoughtful discussion, adding several explanations for which 
he writers are grateful. They agree that all the detailed exploratory studies 
-omprising a major research effort should be made available to those inter- 
sted in the general subject covered. His statement that the original papers 
ire nowhere referred to is a little misleading, in that the referenced “Snow 
tydrology” has some 14 pages filled with references to works by writers in 
he general field of hydrology as well as those engaged in analyzing and re- 
sorting upon the results of the recent laboratory measurements. Some of 
hese were written by Mr. Miller. 

The water-balance approach to forecasting volume of streamflow from 
nelting snow may indeed replace, in time, the index methods now used. How- 
ver, much better coverage must first be established of such elements as 
now water equivalent, soil and snow priming, evaporation and other losses, 
nd precipitation. Until such time, the index methods presently in use must 
ontinue to serve. In fact, until long-term weather forecasts can be made 
7ith improved reliability, forecasts of spring snowmelt runoff will continue to 
e made with considerable uncertainty. 

Mr. Miller chooses to reverse the order of consideration of four equations 
oncerning the heat balance concept, going from a formula containing a large 
umber of terms to one containing only two temperature measurements. For 
resentation to a hydrologist completely familiar with the subject, this order 
f consideration might logically be followed, but the writers assumed that 
ome persons whose knowledge of snow hydrology was limited would read the 
aper. In this case, the usual and preferable method, where a choice exists, 
3 to start with the least-complicated situation and to develop formulas that 
onsider additional parameters, until the most complicated situation is cov- 
red. 

Mr. Miller emphasizes the sparsity of radiation observations, particularly 
yr long wave radiation. Several relatively inexpensive and self-contained in- 
‘ruments for measuring and recording solar radiation have been developed, 
1e by an English concern, based upon research by the Meteorological Office. 
Sa result, it may soon be possible to equip stations in mountainous basins 
ith instruments for measuring daily solar radiation; this information would 
2 reported by the observer to the using offices by electric means in order 
at better forecasts of runoff may be made. 


4 January, 1959, by Oliver A. Johnson and Peter B. Boyer. 


1 Hydr. Engr., North Pacific Div., Corps of Engrs., U.S. Dept. of the Army, Portland, — 


eg. 
2 Hydr. Engr., Portland Dist., Corps of Engrs., U. S. Dept. of the Army, Portland, 
eg. 
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DETERMINATION OF HYDROLOGIC FREQUENCY FACTORS?2 


Discussion by William H. Sammons 


WILLIAM H. SAMMoNSs, ! M. ASCE.—The writer is deeply indebted to the 
author for personal assistance prior to publication of this graphical method 
for the determination of K-values.2 

In April 1958, the writer was “knee-deep” in the preparation of a paper3 
for the May meeting of the American Geophysical Union. A phone call to the 
author solved the major problem which dictated a solution. This will be re- 
ferred to later as Method A and will be presented briefly for the benefit of the 
reader. Methods employed in the May, 1958 paper stimulated the writing of a 
subsequent paper4 in September, 1958. This latter paper presented theory, 
examples and K-values complied from other sources.®,6,7,8,9,10 (A few new 
applications were also presented. 

The K-values presented here as Table 2, “Theoretical Log-Probability 
Critical Values,” will be discussed under Method B, and was designed to sup- 
dlement available tables.4,8 

Method A.—Limited time demanded the expediency of this procedure which 
.§ the author’s suggested method (with a few liberties taken). The writer be- 
lieves the K-values developed (tables of the K-value were furnished by the 
uuthor and a few office copies were printed) served the purpose for which they 
were intended.3 The writer does not recommend the indiscriminate use of the 
inpublished table. Briefly, Method A is as follows: 


4 July, 1959, by Ven Te Chow. 

1 Hydraulic Engineer, U.S. Dept. of Agric., Soil Conservation Service, Central Tech- 
ical Unit, Beltsville, Maryland. 

2 “A Graphical Method for Determining the Log-Normal Hydrological Frequency Fac- 
or,” by Ven Te Chow, (Private communication with Ven Te Chow, April 5, 1958. Not 
or publication). 

3 “Agricultural Watershed Influence on Runoff Potential,” by Melvin H. Kleen, and 
Vm. H. Sammons, presented at the thirty-ninth annual meeting of the May 5-8, 1958 at 
Vashington, D. C. (Manuscript only). eats 

4 «Relationship Between the Original and Transformed Variates for the Logarithmic- 
-robability Law,” by Wm. H. Sammons, (For in-Service use and distribution at the SCS 
lydrology Meeting, Jackson, Miss., October 13-17, 1958). 

5 “On the Logarithmic-Probability Law,” by A. A. Kalinske, Trans. AGU, 1946, Vol. 
7, No. 5, pp. 709-710. 

6 Reference 5 of Proc. Paper 2084, by Ven Te Chow. 

7 “A Nomogram for Log-Normal Frequency Analysis,” by L. L. Weiss, Trans. AGU, 
957, Vol. 38, No. 1, pp. 33-37. : 

8 « Critical Values of the Log-Normal Distribution,” by Jack Moshman, J. A. Statis- 
ical Assoc., 1953, Vol. 48, No. 263 . 600-609. : 
“9 “The Log-Normal Pe oe J. Aitchison and J. A. C. Brown, Cambridge 
iv. Press, 1957. , 

5 16 SP itis a Generalized Log-Normal Distribution to Hydrologic Data,” by D. L. 
rakensiek, Trans. AGU, 1958, Vol. 39, No. 3, pp. 469-473. 
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The author’s procedure was followed verbatim. 


. Figures were plotted similar to the author’s Fig. 1. 


Probability lines were developed for twelve “key” C,-values 
(Cy = 1.1171 - 3.455). 


MINIMUM 
MEAN 
MAXIMUM 


f— O15 
' 


0.215 


0.420 


() 
O03 06609 12 15. 18. 2 24.27 30 0 03 06 09 12 15 18 21 
COEFFICIENTS OF VARIATION 


FIG, 1.—DISTRIBUTION OF THE COEFFICIENTS OF VARIATION, 


4. K-values (graphical determinations) were plotted on the logarithmic 


scale versus Cy-values on the arithmetic scale of semi-logarithmetic 
graph paper. : 


5. K-values for nine levels of probability used in Table 2 of Reference 6 


were read from these graphs as needed for use in Reference 3. 


6. These nine K-values for specific Cy-values were smoothed by com- 


puting X/X; plotting on probability paper; drawing an average line by 
eye; re-computing K-values by the author’s Eq. 1. 
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__ 7. Smoothed K-values were used to calculate X-magnitudes of the hydro- 
logic event. 
8. A table of approximately 100 Cg (or C,,) values were tabulated with K- 
values at the nine levels of probability. 


Fig. 1 is presented to illustrate the extent to which this unpublished tabula- 
tion of K-values was employed.3 

Method B.—K-values developed by Method B could have utilized an IBM- 
650 computer to a great advantage and under a different set of prevailing cir- 
cumstances, it would probably have been worth while as an official project. 
Briefly, Method B is as follows: 


1. The author’s Eq. 9 was modified as follows: 
X, = M exp K a ei Pa cS (9') 


where K-values for Cy = 0.0000 and a few selected values are tabulated 
as Table 1. 


TABLE 1.—THEORETICAL LOG-PROBABILITY CRITICAL VALUES 
FOR THE STANDARD NORMAL VARIATE EXCEEDED WITH 
SPECIFIED PROBABILITY? 


Probability 


Probability % 
in (a) 


in % 


99.9 -3.090232 +3.090232 
99.75 2.807034 2.807034 
99.6 2.652070 2.652070 
99.5 2.5758 2.9758 
99 2.326348 2.326348 
97.5 -1.959964 +1.959964 
96 -1.750868 +1.750868 
95 1.6449 1.6449 
93.5 1.514102 1.514102 
90 1.281552 1.281552 
85 1.036433 1.036433 
-0.674490 +0.674490 


0.000000 


0.000000 


. 
2 Probability in percentage equal to or greater than the given variate. 


Five levels of probability were employed to compute K-values. 
Codex Book Company, Inc., Logarithmic Normal Probability Paper No. 
31,376 was used to plot X/X versus probability at the 99, 84.13, 50, 
15.87 and 1% probability levels (used in smoothing). 

The five levels of probability, “computed” K-values, were plotted on K 
and E, 20 x 20 to the inch, coordinate paper No. 359-10 1/2 LG, versus 
Cx ; 

sindoth curves were drawn and computational errors located and cor- 
rected by computational methods. 

6. The five K-values for specific C, (or C,)-values were smoothed by the 
process described under Method A, step 6. These were later plotted on 


wn 


> 


oy 


Se = > 


110 


February, 1960 


TABLE 2.—THEORETICAL LOG-PROBABILITY CRITICAL VALUES 


HY 2 


Coefficients of 


Probability, in percentage, equal to or greater than the 


given variate 


— —r =T em | 
Skew Variation 99 AK 84.13 AK 50 AK 15.87 AK if AK 
Cs Cy ACy - jo te [ a) + ~ 
3.00 | 0.8177 -- | 1.0437 -- | .76004 -- 27621 -- | .71333 -- | 3.7787 -- 
3.01 | 0.8197 20 | 1.0420 17 | .75929 75 27648 27 | «11235 98 | 3.7804 17 
3.02 | 0.8217 20 | 1.0402 18 | .75855 74 27675 27 |.71139 96 | 3.7822 18 
3.03 | 0.8237 20 | 1.0384 18 | .75780 75 27701 26 | .71043 96 | 3.7840 18 
3.04 | 0.8257 20 | 1.0366 18 | .75696 84 27718 17 | .70948 95 | 3.7863 23 
| 
3.05 | 0.8277 20 | 1.0348 18 | .75622 74 27745 27 | .70853 95 | 3.7882 19 
3.06 | 0.8296 19°} 1.0331 17 | .75544 78 | .27765 20 | .70748 | 105 | 3.7895 13 
3.07 | 0.8316 20 | 1.0313 18 | .75470 74 27790 25 | .70654 94 | 3.7913 18 
3.08 | 0.8336 20 | 1.0296 17 | .75396 74 27816 26 | .70561 93 | 3.7932 19 
3.09 | 0.8356 20 | 1.0278 Sei eose 79 27841 25 |.70449 | 112 | 3.7946 14 
3.10 | 0.8375 19 | 1.0261 17 | .75244 73 27860 19 | .70365 84 | 3.7965 19 
3.11 | 0.8395 20 | 1.0244 17 | .75170 74 27885 25 | .70273 92 | 3.7980 15 
3.12 | 0.8414 19 | 1.0227 It, e510 69 27913 28 |.70171 | 102 | 3.7998 18 
3.13 | 0.8434 20 | 1.0209 18 | .75019 82 27929 16 | .70081 90 | 3.8014 16 
3.14 | 0.8453 19 | 1.0193 16 | .74946 73 27956 27 | .69961 | 120 | 3.8030 16 
3.15 | 0.8473 20 | 1.0176 17 | .74873 73 27981 25 | .69871 90 | 3.8046 16 
3.16 | 0.8492 19 | 1.0159 17 | .74800 73 27999 18 | .69789 82 | 3.8062 16 
3.17 | 0.8511 19 | 1.0142 17 | .74727 73 28017 18 | .69690 99 | 3.8087 16 
3.18 | 0.8531 20 | 1.0126 16 | .74650 77 28040 23 | .69601 89 | 3.8094 16 
3.19 | 0.8550 19 | 1.0109 17 74573 77 28058 18 | .69502 99 | 3.8110 16 
3.20 | 0.8569 19 | 1.0093 16 | .74505 68 | .28084 26 | .69404 98 | 3.8125 15 
3.21 | 0.8588 19 | 1.0077 16 | .74432 73 28102 18 | .69325 79 | 3.8140 15 
3.22 | 0.8608 20 | 1.0060 17 }.14356 76 28125 23 | .69219 | 106 | 3.8156 16 
3.23 | 0.8627 19 | 1.0044 16 | .74288 68 28150 25 | .69123 96 | 3.8172 16 
3.24 | 0.8646 19 | 1.0027 17 | .74212 76 28167 17 | .69026 97 | 3.8184 12 
3.25 | 0.8665 19 | 1.00113} 160 | .74140 12 28184 17 | .68948 78 | 3.8199 15) 
3.26 | 0.8684 19 .99953 | 160 | .74068 72 28210 26 | .68834 | 114 | 3.8214 15 
3.27 | 0.8703 19 .99792 | 161 | .73996 72 28226 16 | .68757 77 | 3.8229 15 
3.28 | 0.8722 19 .99631| 161 | .73921 75 28242 16 | .68662 95 | 3.8242 13 
3.29 | 0.8741 19 .99475 | 156 | .73854 67 28267 25 | .68568 94 | 3.8258 16 
3.30 | 0.8760 19 -99318 | 157 |.73786 68 | .28292 25 | .68474 94 | 3.8273 15 
3.31 | 0.8778 18 99166 | 152 | .73712 74 | .28302 10 | .68388 86 | 3.8289 16 
S327 0.8797 19 -99009 | 157 | .73645 67 | .28326 24 | .68295 93 | 3.8300 11 
3.33 | 0.8816 19 -98851| 158 | .73570 75 | .28342 16 | .68203 92 | 3.8317 17 
3.34 | 0.8835 19 -98695 | 156 | .73503 67 | .28366 24 | .68110 93 | 3.8329 12 
3.35 | 0.8854 19 -98537 | 158 | .73428 75 | .28381 15 | .68019 91 | 3.8344 15 
3.36 | 0.8872 18 .98380 | 157 | .73362 66 | .28400 19 | .67935 84 | 3.8359 15 
3.37 | 0.8891 19 -98234| 146 | .73287 75 | .28415 15 | .67844 91 | 3.8374 15 
3.38 | 0.8909 18 -98086 | 148 | .73221 66 | .28433 18 | .67761 83 | 3.8386 12 
3.39 | 0.8928 19 29'7933)| 153 || .7315% 70 | .28456 23 | .67652 | 109 | 3.8399 13 
3.40 | 0.8946 18 .97780 | 153 | .73085 66 | .28474 18 | .67570 82 | 3.8411 12 
3.41 | 0.8965 19 .97630 | 150 | .73011 74 | .28488 -67480 90 | 3.8426 15 
3.42 | 0.8983 18 -97485 | 145 | .72946 65 .67399 81 | 3.8438 12 
3.43 | 0.9002 19 .97333 | 152 | .72876 70 .67292 | 107 | 3.8451 13 
3.44 | 0.9020 18 .97188 - 72811 -67211 81 | 3.8464 13 
3.45 | 0.9039 .97042 -72737 .67123 3.8476 12 
3.46 | 0.9057 -96895 | 147 | .72672 .67043 80 | 3.8489 13 
3.47 | 0.9075 -96748 | 146 | .72603 -66946 97 | 3.8502] 13 
3.48 | 0.9093 -96603| 145 | .72538 -66866 80 | 3.8514 125 
3.49 | 0.9112 -96452 - 72465 -66780 86 | 3.8526 12 
3.50 | 0.9130 -96310 -72397 -66683 3.8539 13 
3.51 | 0.9148 .96166| 144 | .72332 -66605 3.8552 13 
3.52 | 0.9166 .96024 -72264 -66509 3.8565 | 13 — 
3.53 | 0.9184 -95883 .72199 -66431 3.8576 | 11 
| 0.9202 -95741 -72131 -66336 3.8589 13 
0.9220 .95600 - 72067 -66259 3.8601 
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TABLE 2.—THEORETICAL LOG-PROBABILITY CRITICAL VALUES 


: == 


Coefficients of Probability, in percentage, equal to or greater than the given variate 
7 aa i T T T T T T 
Skew Variation 99 84.13 50 : 
¢ Pan TAGs i AK AK AK saat! AK : AK 
s v v = - + Ay 


3.55 | 0.9220 18 -95600 | 141 | .72067 64 | .28719 16 | .66259 77 | 3.8601 12 
3.56 | 0.9238 18 -95459 | 141 | .71999 68 | .28735 16 | .66165 94 | 3.8614 13 
3.57 | 0.9256 18 -95319 | 140 | .71931 68 | .28750 15 | .66071 94 | 3.8625 11 
3.58 | 0.9274 18 -95180 | 139 | .71867 64 | .28766 16 | .65995 76 | 3.8637 12 
3.59 | 0.9292 18 -95040 | 140 | .71800 67 | .28781 15 | .65902 93 | 3.8649 12 


3.60 | 0.9310 18 -94901 | 139 | .71736 64 | .28796 15 | .65827 75 | 3.8660 11 
3.65 | 0.9399 89 -94215 | 686 | .71408 | 328 | .28865 69 | .65409 | 418 | 3.8712 52 
3.70 | 0.9487 88 -93547 | 668 | .71087 | 321 | .28935 70 | .64990 | 419 | 3.8765 53 
3.75 | 0.9574 87 -92895 | 652 | .70773 | 314 | .29004 69 | .64568 | 422 | 3.8818 53 
3.80 | 0.9661 87 -92249 | 646 | .70455 | 318 | .29064 60 | .64154 | 414 | 3.8868 50 


3.85 | 0.9747 86 -91619 | 630 | .70146 | 309 | .29124 60 | .63755 | 399 | 3.8916 48 
3.90 | 0.9832 85 | .91004/} 615 | .69844 | 302 | .29185 61 | .63352 | 403 | 3.8962 46 
3.95 | 0.9916 84 -90401 | 603 | .69544 | 300 | .29239 54 | .62963 | 389 | 3.9004 42 
4.09 | 1.0000 84 -89806 | 595 | .69248 | 296 | .29291 52 | .62580 | 383 | 3.9045 41 
4.05 | 1.0083 83 -89225 | 581 | .68951 | 297 | .29337 46 | .62194 | 386 | 3.9086 41 


4.10 | 1.0165 82 -88656 | 569 | .68663 | 288 | .29384 47 | .61820 | 374 | 3.9125 39 
4.15 | 1.0247 82 -88083 | 573 | .68363 | 300 | .29429 45 | .61450 | 370 | 3.9162 37 
4.20 | 1.0328 81 -87545 | 538 | .68095 | 268 | .29475 46 | .61065 | 385 | 3.9200 38 
4.25 | 1.0408 80 -87008 | 537 | .67819 | 276 | .29515 40 | .60706 | 359 | 3.9234 34 
4.30 | 1.0488 80 -86477 | 531 | .67545 | 274 | .29553 38 | .60352 | 354 | 3.9266 32 


4.35 | 1.0567 79 -85957 | 520 | .67270 | 275 | .29587 34 | .59993 | 359 | 3.9297 31 
4.40 | 1.0645 78 -85449 | 508 | .67005 | 265 | .29621 34 | .59646 | 347 | 3.9327 30 
4.45 | 1.0723 78 -84947 | 502 | .66741 | 264 | .29654 33 | .59304 | 342 | 3.9356 29 
4,50 | 1.0800 are .84457 | 490 | .66478 | 263 | .29682 28 | .58956 | 348 | 3.9383 27 
4.55 | 1.0877 AG -83971 | 486 | .66221 | 257 | .29715 33 | .58614 | 342 | 3.9408 25 


4.60 | 1.0953 76 .83497 | 474 | .65968 | 253 | .29744 29 | .58282 | 332 | 3.9433 25 
4.65 | 1.1029 76 -83026 | 471 | .65711 | 257 | .29764 20 | .57955 | 327 | 3.9457 24 
4.70 | 1.1104 75 .82567 | 459 | .65465 | 246 | .29793 29 | .57623 | 332 | 3.9480 23 
4.75 | 1.1178 74 .82117| 450 | .65219 | 246 | .29811 18 | .57315 | 308 | 3.9503 23 
4.80 | 1.1252 74 -81673 | 444 | .64977 | 242 | .29833 22 | .56996 | 319 | 3.9523 20 


4.85 | 1.1325 73 .81239 | 434 | .64740 | 237 | .29858 25 | .56672 | 324 | 3.9542 19 
4.90 | 1.1398 73 .80808 | 431 | .64501 | 239 | .29875 17 | .56352 | 320 | 3.9559 17 
4.95 | 1.1470 72 .80387| 421 | .64266 | 235 | .29889 14 | .56056 | 296 | 3.9575 16 
5.00 | 1.1542 72 -79970 | 417 | .64035 | 231 | .29907 18 | .55749 | 307 | 3.9591 16 
5.10 | 1.1684 | 142 .79159| 811 | .63583 | 452 | .29934 27 | .55156 | 593 | 3.9622 31 


5.20 | 1.1824 | 140 .78373 | 786 | .63137 | 446 | .29958 24 | .54572 | 584 | 3.9650 28 
5.30 | 1.1962 | 138 .77612 | 771 | .62705 | 432 | .29978 20 | .54011 | 561 | 3.9671 21 
5.40 | 1.2098 | 136 .76875 | 737 | .62281 | 424 | .29995 17 | .53445 | 566 | 3.9689 18 
5.50 | 1.2232 | 134 .76161| 714 | .61869 | 412 | .30009 14 | .52900 | 545 | 3.9706 17 
5.60 | 1.2365 | 133 .75463 | 698 | .61463 | 406 | .30018 9 | .52371 | 529 | 3.9718 12 


1.2496 | 131 .61067 | 396 | .30024 6 | .51848 
1.2625 | 129 -60682 | 385 | .30028 4 | .51345 
1.2753 | 128 -60302 | 380 | .30028 0 | .50855 
2 
2 
8 


1.2879 | 126 .59932 | 370 | .30026 -50371 
1.3004 | 125 .59567 | 365 | .30024 .49874 


49419 


.59212 | 355 | .30016 


6.20 | 1.3127 

6.30 | 1.3249 | 122 .58863 | 349 | .30008 8 | .48964 
6.40 | 1.3369 | 120 .58522 | 341 | .29995 13 | .48524 
6.50 | 1.3488 .58187 | 335 | .29983 12 | .48084 


.29969 14 | .47658 


6.60 | 1.3605 .57861 | 326 
6.70 | 1.3722 .57537 129955 47228 
6.80 | 1.3837 57223 .29939 .46812 
6.90 | 1.3950 56913 (29919 46409 
7.00 | 1.4063 .56611 4.29902 46012 


.29879 .45629 


1.4174 -56313 
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TABLE 2.—THEORETICAL LOG-PROBABILITY CRITICAL VALUES (Continued) 


HY 2 


Coefficients of 


Probability, in percentage, equal to or gre 


ater than the given variate 


Skew | Variation 99 ak |84=18 |) ag feo AK | 15.87 AK 1 AK 
Ge Cy, [ACy - - eee Sz. ret 

7.20 | 1.4284 | 110 | .66554 | 458 | .56022 | 291 | .29859 | 20 | .45243 | 386 | 3.9670 | 14 
7.30 | 1.4394 | 110 | .66102 | 452 | .55732 | 290 | .29833 | 26 | .44874 | 369 | 3.9656 | 14 
7.40 | 1.4501 | 107 | .65667 | 435 | .55453 | 279 | .29813 | 20 | .44498 | 376 | 3.9640 | 16 
7.50 | 1.4608 | 107 | .65238 | 429 | .55175 | 278 | .29785 | 28 | .44149 | 349 | 3.9623 | 17 
7,60 | 1.4718 | 106 | .64817 | 421 | .54904 | 271 | .29762 | 23 | .43787 | 362 | 3.9603 | 20 
7,70 | 1.4819 | 105 | .64407 | 410 | .54635 | 269 | .29734 | 28 | .43438 | 349 | 3.9585 | 18 
7.70 | 1.4819 | 105 | .64407 | 410 | .54635 | 269 | .29734 | 28 | .43438 | 349 | 3.9585 | 18 
7.80 | 1.4923 | 104 | .64003 | 404 | .54373 | 262 | .29709 | 25 | .43096 | 342 | 3.9564] 21 
7.90 | 1.5026 | 103 | .63608 | 395 | .54114 | 259 | .29679 | 30 | .42768 | 328 | 3.9543 | 21 
8.00 | 1.5127 | 101 | .63226 | 382 | .53862 | 252 | .29652 | 27 | .42428 | 340 | 3.9521 | 22 
8.10 | 1.5228 | 101 | .62848 | 378 | .53613 | 249 | .29622 | 30 | .42115 | 313 | 3.9499 | 22 
8.20 | 1.5328 | 100 | .62479 | 369 | .53368 | 245 | .29593 | 29 | .41797 | 318 | 3.9476 | 23 
8.30 | 1.5427 | 99 | .62115 | 364 | .53127 | 241 | .29562 | 31 | .41486 | 311 | 3.9453 | 23 
8.40 | 1.5526 | 99 | .61757 | 358 | 52886 | 241 | .29531 | 31 | .41169 | 317 | 3.9430 | 23 
8.50 | 1.5623 | 97 | .61409 | 348 | .52655 | 231 | .29503 | 28 | .40861 | 308 | 3.9407 | 23 
8.60 | 1.5719 | 96 | .61068 | 341 | .52425 | 230 | .29470 | 33 | .40569 | 292 | 3.9384 | 23 
8.70 | 1.5815 | 96 | .60731 | 337 | .52198 | 227 | .29437 | 33 | .40282 | 287 | 3.9359 | 25 
8.80 | 1.5910 | 95 | .60401 | 330 | .51975 | 223 | .29405 | 32 | .39990 | 292 | 3.9335 | 24 
8.90 | 1.6004 | 94 | .60079 | 322 | .51757 | 218 | .29375 | 30 | .39704 | 286 | 3.9310 | 25 
9.00 | 1.6097 | 93 | .59760 | 319 | .51541 | 216 | .29340 | 35 | .39434 | 270 | 3.9284 | 26 
9.10 | 1.6189 | 92 | .59451 | 309 | .51328 | 213 | .29307 | 33 | .39159 | 275 | 3.9259 | 25 
9,20 | 1.6281 | 92 | .59143 | 308 | .51119 | 209 | .29276 | 31 | .38878| 281 | 3.9233 | 26 
9.30 | 1.6372 | 91 | .58842 | 301 | .50913 | 206 | .29241 | 35 | .38622 | 256 | 3.9205 | 28 
9.40 | 1.6462 | 90 | .58549 | 293 | .50710 | 203 | .29208 | 33 | .38361| 261 | 3.9178 | 27 
9.50 | 1.6552 | 90 | .58256 | 293 | .50509 | 201 | .29174 | 34 | .38103 | 258 | 3.9151 | 27 
9.60 | 1.6640 | 88 | .57974 | 282 | .50313 | 196 | .29141 | 33 | .37853 | 250 | 3.9123 | 28 
9.70 | 1.6729 | 89 | .57688 | 286 | .50116 | 197 | .29106 | 35 | .37603 | 250 | 3.9095 | 28 
9.80 | 1.6816 | 87 | .57414 | 274 | .49925 | 191 | .29071 | 35 | .37360 | 243 | 3.9067] 28 
9.90 | 1.6903 | 87 | .57143 | 271 | .49736 | 189 | .29039 | 32 | .37110 | 250 | 3.9040 | 27 
10.00 | 1.6989 | 86 | .56875 | 268 | .49549 | 187 | .29003 | 36 | .36874 | 236 | 3.9010 | 30 
10.10 | 1.7074 | 85 | .56615 | 260 | .49366 | 183 | .28968 | 35 | .36643 | 231 | 3.8980 | 30 
10.20 | 1.7159 | 85 | .56355 | 260 | .49183 | 183 | .28933 | 35 | .36404 239 | 3.8951 | 29 
10.30 | 1.7243 | 84 | .56101 | 254 | .49006 | 177 | .28900 | 33 | .36180 | 224 | 3.8921 | 30 
10.40 | 1.7327 | 84 | .55851 | 250 | .48827 | 179 | .28864 | 36 | .35948 | 232 | 3.8890 | 31 
10.50 | 1.7410 | 83 | .55604 | 247 | .48654 | 173 | .28830 | 34 | .35730 | 218 | 3.8860 | 30 
10.60 | 1.7493 | 83 | .55359 | 245 | .48480 | 174 | .28795 | 35 | .35505 | 225 | 3.8830 | 30 
10.70 | 1.7574 | 81 | .55123 | 236 | .48313 | 167 | .28761 | 34 | .35295 | 210 | 3.8800 | 30 
10.80 | 1.7655 | 81 | .54887 | 236 | .48146 | 167 | .28727 | 34 | .35079 | 216 | 3.8769 | 31 
10.90 | 1.7736 | 81 | .54654 | 233 | .47979 | 167 | .28690 | 37 | .34873 | 206 | 3.8738] 31 
11,00 | 1.7816 | 80 | .54427 | 227 | .47816 | 163 | .28655 | 35 | .34671 | 202 | 3.8707 | 31 
11.10 | 1.7896 | 80 | .54199 | 228 | .47654 | 162 | .28622 | 33 | .34460 | 211 | 3.8676 | 31 
11,20 | 1.7975 | 79 | .53978 | 221 | .47495 | 159 | .28586 | 36 | .34256 | 204 | 3.8646 | 30 
11.30 | 1.8053 | 78 | .53760 | 218 | .47339 | 156 | .28551 | 35 | .34063 | 193 | 3.8615 | 31 
11.40 | 1.8131 | 78 | .53546 | 214 | .47183 | 156 | .28516 | 35 | .33863 | 200 | 3.8584] 31 
11.50 | 1.8209 | 78 | .53331 | 215 | .47028 | 155 | .28482 | 34 | .33665 | 198 | 3.8552 | 32 
11.60 | 1.8286 | 77 | .53121 | 210 | .46877 | 151 | .28448 | 34 | .33470| 195 | 3.8522 | 30 
11.70 | 1.8362 | 76 | .52916 | 205 | .46728 | 149 | .28414 | 34 | .33279 | 191 | 3.8491 | 31 
11.80 | 1.8439 | 77 | .52710 | 206 | .46578 | 150 | .28378 | 36 | .33096 | 183 | 3.8460 | 31 
11,90 | 1.8514 | 75 | .51510 | 200 | .46433 | 145 | .28343 | 35 | .32919 | 177 | 3.8429 | 31 
12.00 | 1.8589 | 75 | .52312 | 198 | .46289 | 144 | .28310 | 33 | .32726 | 193 | 3.8398 | 31 
12.10 | 1.8663 | 74 | .52117 | 195 | .46148 | 141 | .28277 | 33 | .32553 | 173 | 3.8368 | 30 
12.20 | 1.8738 | 75 | .51921 | 196 | .46004 | 144 | .28242 | 35 | .32371 | 182 | 3.8337 | 31 
12.30 | 1.8811 | 73 | .51733 | 188 | .45866 | 138 | .28207 | 35 | .32203 | 168 | 3.8306] 31 
12.40 | 1.8884 | 73 | .51546 | 187 | .45728 | 138 | .28172 | 35 | .32027| 176 | 3.8274 | 32 
12.50 | 1.8957 | 73 | .51358 | 188 | 45591 | 137 | .28138 | 34 | .31862] 165 | 3.8243] 31 
12,60 | 1.9030 | 73 | .51174 | 184 | .45455 | 136 | .28103 | 35 | .31689 | 173 | 3.8213 
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TABLE 2.—THEORETICAL LOG-PROBABILITY CRITICAL VALUES (Continued) 


Coefficients of Probability, in percentage, equal to or greater than the given variate 
Variation 99] 84.13 7" 50 15.87 Re 
2 AK | = AK : AK + AK + AK 


- 
& 


a ht Pe 


5 =. 
T 50993 | 181 | .45321 | 134 |.28070 | 33 | .31519 170 | 3.8182 | 31 
-50815 | 178 | .45191 | 130 | .28036 | 34 | .31361 | 158 | 3.8151 | 31 
50626 | 189 | .45047 | 144 | .27995 | 41 | .31187 | 174 | 3.8120 | 31 
-50463 | 163 | .44930 | 117 | .27969 | 26 | .31030 | 157 | 3.8089 | 31 
-50122 | 341 | .44676 | 254 | .27902 | 67 | .30717 | 313 | 3.8029 | 60 


1.9593 | 138 -49788 | 334 | .44429 | 247 | .27836 66 | .30411 | 306 | 3.7967 62 
1.9730 | 137 -49462 | 326 | .44185 | 244 | .27771 65 | .30103 | 308 | 3.7905 62 
1.9866 | 136 -49142 | 320 | .43946 | 239 | .27705 66 | .29804 | 299 | 3.7844 61 
2.0000 | 134 -48831 | 311 | .43712 | 234 | .27639 66 | .29513 | 291 | 3.7783 61 
2.0133 | 133 -48526 | 305 | .43482 | 230 | .27574 65 | .29231 | 282 | 3.7721 62 
2.0264 | 131 -48228 | 298 | .43258 | 224 | .27509 65 | .28956 | 275 | 3.7659 62 


2.0264 | 131 | .48228 | 298 | .43258 | 224 | .27509 65 | .28956 | 275 | 3.7659 62 
2.0394 | 130 -47938 | 290 | .43037 | 221 | .27447 62 | .28678 | 278 | 3.7598 61 
2.0522 | 128 -47654 | 284 | .42822 | 215 | .27382 65 | .28422 | 256 | 3.7536 62 
2.0650 | 128 -47373 | 281 | .42609 | 213 | .27321 61 | .28154 | 268 | 3.7475 61 
2.0776 | 126 -47100 | 273 | .42401 | 208 | .27257 64 | .27900 | 254 | 3.7413 62 


2.0900 | 124 -46835 | 265 | .42198 |; 203 | .27195 62 | .27651 | 249 | 3.7351 62 
2.1020 | 120 -46580 | 255 | .42003 | 195 | .27136 59 | .27411 | 240 | 3.7290 61 
2.1146 | 126 -46316 | 264 | .41801 | 202 | .27074 62 | .27165 | 246 | 3.7228 62 
2.1268 | 122 -46063 | 253 | .41606 | 195 | .27013 61 | .26928 | 237 | 3.7167 61 
2.1388 | 120 -45816 | 247 | .41416 | 190 | .26952 61 | .26696 | 232 | 3.7104 63 


2.1507 | 119 -45575 | 241 | .41230 | 186 | .26894 58 | .26468 | 228 | 3.7045 59 
2.1625 | 118 -45337 | 238 | .41046 | 184 | .26834 60 | .26244 | 224 | 3.6986 58 
2.1742 | 117 -45104 | 233 | .40865 | 181 | .26775 59 | .26025 | 219 | 3.6928 57 
2.1858 | 116 -44876 | 228 | .40688 | 177 | .26717 58 | .25814 | 211 | 3.6871 56 
-44652 .40513 | 175 | .26658 59 | .25600 | 214 | 3.6815 57 


-44433 -40343 | 170 | .26602 56 | .25399 | 201 | 3.6758 55 
113 .44216 | 217 | .40173 | 170 | .26545 57 | .25192 | 207 | 3.6703 55 
112 -44004 | 212 | .40007 | 166 | .26488 57 | .24996 | 196 | 3.6648 52 
111 -43795 | 209 | .39844 | 163 | .26434 54 | .24796 | 200 | 3.6596 52 
43588 -39681 | 163 | .26377 57 | .24598 | 198 | 3.6538 58 


-43387 .39523 | 158 | .26322 55 | .24411 | 187 | 3.6483 55 
108 -43190 | 197 | .39367 | 156 | .26268 54 | .24220 | 191 | 3.6427 56 
108 42994 | 196 | .39213 | 154 | .26214 54 | .24038 | 182 | 3.6372 55 
107 -42801 | 193 | .39061 | 152 | .26160 54 | .23858 , 180 | 3.6317 55 
-42614 -38913 | 148 | .26107 53 | .23689 | 169 | 3.6261 56 


42429 .38766 | 147 | .26055 52 | .23509 | 180 | 3.6206 55 


-42236 | 193 | .38613 | 153 | .26000 55 | .23334 | 175 | 3.6152 54 

42065 .38477 | 136 | .25950 .23167 | 167 | 3.6096 56 

-41888 .38336 -25898 -23000 56 

.37661 -25647 +22205 263 

-37029 -25406 .21465 241 

: -36436 -25173 .20781 250 
2.5397 -35878 24949 -20134 
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probability paper, as noted under step 3 above, asa check on linearity. 
Forty Cg (or Cy) values proportionally spaced throughout Table 2 con- 
stituted the sample tested. 

7, First differences, AK are shown in Table 2 as an aid to prospective 
users. 


Needless to say, the five computed K-values for each entry of Cg will al- 
low further probability levels to be obtained by the author’s Fig. 1 method 
utilizing tabulated K-values for specified Cy values and solving the author’s 
Eq. 1. (If publication space were available the 24 columns, 32 lines and 8 
pages of tabulations of the individual steps could be made available). Table 2 
does not contain an additional forty-eight computed Cg—K-values which have 
been computed to date for specific problems and to satisfy the curiosity of the 
writer. A brief summary of similar work by the writer as noted above is 
available. 

Precaution Necessary.—As pointed out in detail and in examples,4 K-val- 
ues may be used with the original variate or with a log-transformation vari- 
ate. The latter is more efficient, statistically speaking, than the scalar trans- 
formation of the original variate. H. Alden Foster suggests!1 a method for 
making theoretical corrections to Cy and Cg. Similar work by the writer also 
supports this generalized approach. 


ie a nT ee eee pe SS ea a ee ee a eS 


11 Discussion of “The Log-Probability Law and Its Engi i i 
gineering Applications,” : 
Alden Foster, 1959, Proc. No. 709, Vol. 81, pp. 14-19. ee a 


HY 2 115 


ROLL WAVES AND SLUG FLOWS IN INCLINED OPEN CHANNELS? 


Discussion by R. Hugh Taylor and John F. Kennedy 


R. HUGH TAYLOR,! A.M. ASCE and JOHN F. KENNEDY,2 A.M. ASCE.— 
The author has explored an obscure byway of hydraulics and has presented an 
interesting body of experimental data resulting from what was obviously a 
painstaking and thorough laboratory investigation. While knowledge is valua- 
ble per se, the fact that the phenomena dealt with here are so little known 
confers an additional value upon the results, they offer the engineer a means 
of broadening his intuition in an area he is not likely to have experienced di- 
rectly. In this respect, then, the author has done the profession a definite 
service. 

It is to be regretted, however, that the author’s treatment of his results 
leads one to conclude that they are more directly useful than they actually 
are. This misunderstanding seems to have two different sources. First, the 
choice of nomenclature is unfortunate in that the term “roll waves” already 
has a generally accepted meaning, if not a precisely defined one. In spite of 
this, however, the author has chosen to apply the term exclusively to insta- 
bilities of laminar flow, a phenomenon which is of predominately academic 
interest to the hydraulic engineer. Furthermore, the term is so defined in the 
paper as to exclude the meanings already well established in the literature, 
notably in the writings of Thomas? and of Dressler.5 Second, the confusion is 
enhanced by the reference to impact and freeboard problems in channel de- 
sign caused by instabilities of flow. Certainly, the laminar instabilities the 
author calls “roll waves” are not of this magnitude, if they do not occur at 
Reynolds numbers greater than 415. 

The author has likewise limited himself by choosing to include within his 
term “slug flow” only those unsteady flows which have arisen more or less 
directly from laminar flow. This may well be a reservation without a reason, 
if the properties of the limiting wave are not influenced by the presence or 
absence of prior laminar flow. Unless this is true, however, the similarity of 
these waves to the roll waves reported by Cornish13 or those discussed by 
Dressler® is questionable. It seems improbable that the flow with roll waves 
which Dressler attributes to Cornish, or the flow with roll waves observed by 
V. A. Vanoni and the writers in the Santa Anita Flood Control Channel, near 
Arcadia, Calif., could ever have been laminar. This is not to object to the 
scope of the paper, but merely to question the relation between the actual 


scope and that which is suggested by the title and introduction. 


ES ee a ee 
4 July, 1959, by Paul G. Mayer. ; 

1 Dept. of Civ. Engrg., Calif. Inst. of Tech., Pasadena, Calif. 

2 Dept. of Civ. Engrg., Calif. Inst. of Tech., Pasadena, Calif. 
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Purpose.—The purpose of this discussion is three-fold and may be sug- 
gested by a passage from the paper itself: “Rigorous mathematical treatment 
of complex hydraulic phenomena encounters great difficulties. Convenient ap- 
proximations often lead to equations of questionable veracity. The background 
of experimental reality provides important insight into physical phenomena.” 
Specifically, the writers’ intent is to call attention to some theoretical argu- 
ments in a recent English paper, to enter an objection to what the author pre- 
sents as the theory of “ roll waves” and to the presentation and analysis of his 
experimental results, and to mention some recently reported experiments 
which supplement the findings of the author. The objections will be entered 
first. 

On the Theory of "Roll Waves" and the Presentation and Analysis of Data.— 
The first objection is that the theoretical development is weak. If oneassumes 
irrotational flow and surface disturbances small in amplitude with respect to 
their wave length, then the linear analysis of Airy2 follows. It is well known 
that in this analysis the celerity, c, of a wave is given by 


c2 = (a r+ & ) tanh (Dl ee ee eee (i) 


where J" = - = “kinematic surface tension,” so called by analogy with kine- 


SS aos ; 27 
matic viscosity, anda = res wave number 


Under the usual simplifications, tanh (a D) is taken equal to a D whena D 
is small, and it is taken as unity when a D is large (largerthan \/2, usually). 
In the former case one speaks of shallow water waves, the expression for the 
celerity of which has no useful minimum; that is, 


Cshallow = VT a2 D+ gD 
from which 


dco" aD 
ae Vl a2D+geD 


The latter vanishes only for a = 0, which requires an infinite wave length or 
no disturbance at all. Of the two simplifications, it is only in the latter one, 
where one speaks of deep water waves, that one finds a minimum clerity and 
this case is certainly irrelevant to the problem at hand. The analysis in Eqs. 
(1) through (6) is for the deep water case, and since the author makes almost 
no use of his discussion of wave celerity, perhaps it would have been better to 
omit it completely. 

The author’s criterion for the formation of progressive waves is very hard 
to accept in the present situation. If the progressive waves associated with 
“roll waves” and “slug flow” are treated as discontinuous surge fronts, and 


the velocity is assumed constant with depth, the celerity of the waves can be 
shown to be 


1D D | 
Cr=ayi en oe (+3) lisse toy Open Te OLE 
where c is the velocity of the wave relative to the undisturbed flow, D repre- 


sents the depth of undisturbed flow, and Do is the depth upstream from the 
disturbance. If the disturbance is vanishingly small, that is, if D=Dpo, then 
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the celerity becomes that given by the author in Eq. 11, c = Vg D. The same 
limiting value is given by Eq. i for shallow water waves if surface tension is 
neglected and if the disturbance is very small compared to the wave length 
and the depth, that is, in the deep water case. 

For Eq. ii to be valid, the velocity must not vary with depth, upstream and 
downstream from the disturbance. If some other velocity profile, such as the 
parabolic one used by the author, exists upstream and downstream from the 
disturbance, the wave celerity is further modified, as was shown by 
Boussinesq,+° to 


c= VED (143 2)h EB (n2 «2 wp) BANG setae (iii) 


where 7) is the amplitude of the disturbance, Ng represents the Froude Num- 
D 
ber of the undisturbed flow, B is the momentum coefficient = sf (v /V)2ay 
fe) 


and c denotes the velocity of the disturbance relative to the average velocity 
of the undisturbed flow. 

The limiting case is again c = Vg D, but only if 7 =0 and B =1 (a rectangu- 
lar velocity profile). Acomparison of Figs. 17 and 23 indicates that in the ex- 
periments reported by the author, the depth of flow and the wave heights were 
of the same order of magnitude. Noticing also that the Froude numbers were 
quite high, and recalling that the momentum coefficient has the value 1.20 for 
a parabolic velocity distribution, one must conclude that the use of vgD as 
the wave velocity relative to the average velocity of flow is incorrect. 

Actually, taking the Froude number as 1.7, the amplitude of the disturb- 
ance as half the undisturbed depth, and a parabolic velocity distribution, the 
resulting disturbance velocity is only about half that given by the limiting 
case. Similarly, the celerity of oscillatory waves (Eq. i) is greatly modified 
if the velocity profile departs appreciably from rectangular. Thus Eq. 11 and 
the criteria stated in Eqs. 12a and 17 must be rejected. 

The kinematic argument of Lighthill and Whitham’ properly predicts a 
celerity of very long, “kinematic” waves (in which dynamic considerations 
are unimportant) equal to twice the mean stream velocity. But to convert this 
to the statement that the Froude number equals 2 is equivalent to stating that 
the wave celerity is given by Vg D, which was shown above to be incorrect for 
the problem at hand. Lighthill and Whitham themselves said: “It must be re- 
marked that the result Np < 2 has been deduced under the assumption that the 
disturbance is small. More generally,..., the value of Np which must be 
exceeded if a bore of constant strength is to be maintained, depends on the 
strength; this value of Nf is always less than 2 and tends to 2 as the strength 
approaches zero.” 

The second objection to be entered is that the author’s analysis of the ex- 
perimental data is somewhat misleading. Where dimensionless parameters 
are used, it is generally desirable to use the standard ones for which the 
reader already has an intuitive “feel,” and to resist the strong temptation to 
invent new ones unless the gain of greater simplicity of the new ones out- 
weighs theadvantages of familiarity with the established ones. The parameter 
presented in Eq. 29 which he calls the critical flow number Fey, is actually a 


modified Weber number: 


; 18 «Rssai sur la Théorie des Eaux Courantes,” Mémoires Divers Savants a l’Acade- 


mie des Sciences, 1877, Vol. 23. 


: 
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This results from the tacit assumption that the Froude number is unity in the 
development of the critical flow number in Eq. 29. 

The critical depth number, Np, which is given in Eq. 31, can also be con- 
sidered as a shear velocity Reynolds number, 


where ux = shear velocity = Vg DS. The Darcy-Weisbach friction factor, f, is 
defined by 


S°= <p Aid a iMe ea a) cana eee eens (vi) 
and can be expressed as 

ie a( 3)” EP Se MEA ‘eer eee (vii) 
or 

u, = VV£/8 
from which 

Nit Regi! Sr S49) stag wine ae (viii) 


In the laminar regime, where a parabolic velocity profile is applicable, Eqs. 
10a and (vi) yield 


fie: APRs). to Jndlnvanr Stee (ix) 
Substituting this value of f in Eq. viii gives 
Np STV SAR D Pb d > A SE SS ee (x) 


In the present experiments, a single fluid was used and thus the only vari- 
ables in the critical flow number are V, D, and S. Thus Fig. 16, in which Np 
is plotted against F,,, essentially relates V, S, and D, and would be expected 
to give a linear relation in the laminar region ona logarithmic plot. With a 
suitable manipulation of variables, Fig. 16 can be put in the form of the usual 
pipe friction diagram (except for the constants involving p and o) which is 
found in any elementary text on fluid mechanics.!9 For flow of this type, the 
laminar region ends at approximately Rp = 600 to 1,000 (although laminaz 
flow can exist at Reynolds numbers up to 2,500 under appropriate conditions) 
which corresponds to Np = 42.5 to 54.8. In Fig. 16, this is about the region i1 
which the straight line relation ends, as would be expected. It is also the up- 
per limit of the Reynolds number for which roll waves were observed and this 
confirms the conclusion that “roll waves” are confined to the laminar regime 
of flow. 


19 «Fluid Mechanics,” by R. L. Daugherty, and A. C. Ingersoll, McGraw-Hill, 195! 
p. 182. ; 
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: In a similar manner, Figs. 12, 13, 14, 15, and 18 are also essentially pipe- 
Iriction diagrams for this type of flow. This becomes more apparent when it 


is considered that 
1 


NF ( 8 i 
LSet 
Indeed, if use is made of Eq. ix the result for the laminar regime is 


1 
Ne (° 2). 
dGue)-\) 24 


Np = VS/3 VRp 


or 


which is exactly the relation the author presents as Eq. 19 and in Fig. 18. The 
runs conducted outside of the laminar regime are in the transition regime 
from laminar to turbulent flow, thus the friction factor cannot be obtained di- 
rectly from a pipe friction diagram, since the relation between Rp and f in 
this region is not well defined. Actualiy, it depends on the manner in which the 
transition occurs, which in turn depends on the manner in which the experi- 
ments were conducted. Had the experiments progressed into the fully turbu- 
lent regime, the relation between S, V, and D could be obtained directly from 
a pipe friction diagram. Thus it appears that this phase of the author’s analy- 
sis consists primarily of presenting the well known pipe friction diagram in 
several alternate forms. 

The author’s statement that the minimum slope for ‘ roll wave” formation 
is 3°/o was not actually verified by his experiments since no runs were made 
with slopes in the neighborhood of 39/o. Actually Figs. 16, 17, and 18 would 
indicate that a better criterion for the formation of “roll waves” and slug flow 
as defined by the author would be that “ roll waves” are associated with lami- 
nar flow, while slug flow occurs when the flow is in the transitional regime 
between laminar and turbulent flow. The writers also question the applica- 
bility of the statement: “All Froude numbers leading to the formation of slug 
flows were above the theoretical minimum NF = 2'', The criterion for slug 
flow formation mentioned above would indicate that slug flow could result even 
if Np < 2 if the flow were transitional or turbulent. No experiment was con- 
ducted in this range to either prove or disprove the necessity of turbulence 
for slug flow formation regardless of the size of the Froude number. Such an 
experiment would be very enlightening. 

Comparison with Other Recent Papers.—There are two papers (by A. M. 
Binnie20 and T. B. Benjamin21) which appeared at about the same time* as 
Re nah ek, gl tps ol pees hee ee a a ee 


20 «Instability in a Slightly Inclined Water Channel,” by A. M. Binnie, J. Fluid Mech., 


1959, vol. 5, pp. 561-570. F 

21 «wave Formation in Laminar Flow Down an Inclined Plane,” J. Fluid Mech., by T. 
B. Benjamin, 1957, vol. 2, pp. 554-574; Numerical corrections ibid. 1958, vol. 3, p. 657. 
Zz The actual precedence is somewhat confusing. The theoretical paper of Benjamin 
was published before the present paper (but after the thesis upon which it is based). 
Similarly, the experimental results of Binnie appeared after the oral presentation of the 
present paper but before it was published. As an aside, it might be suggested that the 
date of submission of a paper be shown in case questions of priority ever should arise. 


This is common procedure in many technical journals. 


‘7 
é 


a 


120 February, 1960 HY 2 


the present one, and which are relevant to this discussion. They are men- 
tioned at some length here because they appeared in a relatively new publica- 
tion, the Journal of Fluid Mechanics, which is published in England, and may 
not be familiar to some readers. The two papers are essentially contempo- 
rary with the author’s; it is interesting to see if they are in agreement. 

The English papers supplement the present one in that (a) they are con- 
cerned with the lower threshold of waves in laminar flow, (b) Binnie reports 
having explored the spectrum of Reynolds numbers from before occurrence of 
waves up into the turbulent patch or “slug flow” condition, and (c) Benjamin 
has developed a simplified stability analysis, not dependent upon the Airy as- 
sumptions, which Binnie has applied to laminar flow in slightly inclined open 
channels. 

A brief review of Benjamin’s stability analysis is presented by Binnie, as 
adapted for small slopes. It is essentially as follows: Define the following di- 
mensionless quantities: 


Celerity C =c/V = Cy, +i Cj, Cr and C; both being real 
Wave number a =a D 

Distance X = x/D, Y = y/D 

Time T= t.D/¥. 


Now if one imposes a small disturbance Y = 6 exp fa i(x.- C T)} on the 
laminar flow represented by Eqs. 7 through 10, there results an equation 
known as the Orr-Somerfeld equation. 22 This equation can be solved in power 
series in Rp which converges fairly rapidly for small Rp, although the alge- 
bra is quite tedious. For very small waves one finds Cy = 2 (thus confirming 
Lighthill and Whitham under the given assumptions), and 

Cy = 5 Rp(F a na as Ale ee eee ) 

Dv2 vy2 

where @ is the inclination of the channel. Within the small-wave assumption, 
all this equation can show is the relation between Rp, D, and @ under which 
the waves would start to grow exponentially. The condition for stability is that 
Gj < 0, so that the disturbance will have no part which grows exponentially 
with time. Thus, for instability, we must have 


azl 
D v2 


D 8 
f= COSeO <a 
v2 Sie 


Of interest is the wave number, Qy, for the fastest growing instability. This 
is obtained by maximizing a Cj with respect to a and is given by 


If this last result, together with the primary flow Eqs. 8 and 10 be substituted 
into Eq. xiii, the maximum value of Cj is seen to be 

eee 
a “The Theory of Hydrodynamic Stability,” by C. C. Lin, Cambridge Univ. Press. 
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ay adr2/ 3 
max Cj = 0.25( se | ~F5 (e sin 0 
£5 


cot ¢ 1/6, 11/6 
pase Ta 


5 sta (Ghee 


[tis seen that the minimum Rp for which instability can exist is obtained by 
setting max C; = 0. This gives 


Rp = 2 COteUe me apes See eiee sae ces (xv) 


Because of the linearizations involved in its derivation, the development 
mentioned above cannot lead to definite conclusions concerning fully-developed 


TABLE 1.—THRESHOLD OF INSTABILITY. COMPARISON OF CALCULATED 
WITH OBSERVED RESULTS, IN TERMS OF DEPTH REYNOLDS NUMBERS 


é 
(degrees) 


Note: Rp' is given by Rp’ =3 cot 9, Eq. xv. Ry’ is the threshold reported by Binnie. 


‘roll waves.” However, it can indicate a method of attack, and it can also be 
used to give a lower bound to the values of Rp below which unstable waves, 
and the “roll waves” arising from them, cannot occur. For this purpose the 
Reynolds numbers corresponding to neutral stability (that is, to Cj = 0) were 
salculated, and the values superimposed on the author’s Fig. 17. Binnie’s ob- 
served thresholds are also shown, for 9 = 1°, 2°, and 23/4°. (See Fig. 27.) 
These values are well below the author’s observed values and therefore can 
9€ said to be consistent. 


—s i nana 
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FIG, 27.—COMPARISON OF MAYER’S RESULTS WITH LOWER THRESHOLDS 
REPORTED BY BINNIE (CF, FIG, 17). 


FLOOD CONTROL ASPECTS OF CAUCA VALLEY DEVELOPMENT? 
Discussion by Gordon R. Williams 


GORDON R. WILLIAMS, ! F. ASCE.-—The planning and final design of water 
resources developments should be preceded by a thorough study of the basic 
meteorologic and climatologic characteristics of the general region. Such a 
study should first lead to a classification of the climate according to type 
along the lines suggested by Képpen, Thornwaite and others. It is evident that 
such a study was made for the Cauca Valley, and although not specifically 
mentioned, the climate falls into Képpen’s “tropical wet and dry” or “savanna” 
classification. Such climates extend along the entire western coast of Central 
America and into Mexico. 

The purpose in classifying the climate as to type is to determine in a gen- 
eral way the potentialities of the climate, including the possible occurrence of 
different storm types, the probable severity and duration of droughts, and the 
possibility of killing frosts and damaging winds. For example, the tropical 
wet and dry climate can be expected to havea total annual precipitation which 
will be enough to supply the annual evapo-transpiration to be expected from 
consideration of the regional temperature and humidity, but the distribution of 
the precipitation month by month will vary widely and there will be seasonal 
deficiencies which will result in crop failures without irrigation. The Cauca 
Valley definitely has this type of climate and, therefore, irrigation works are 
necessary to develop the agricultural potential of the region. 

The paper reveals that the storm rainfall potentials in the region are low, 
relative to many other parts of the world, because of the absence of hurri- 
canes or cylconic disturbances associated with cold fronts. As a measure of 
this characteristic, it is pointed out that the maximum 24-hr rainfall has 
never exceeded 6% of the average annual rainfall. This is in contrast to the 
east coast of the United States where 24-hr storm rainfalls equal to 15% to 
20% of the average annual rainfall have occurred over areas of several thou- 
sand square miles. Presumably most of the precipitation in the Cauca Valley 
results from convective cells, which are characteristics of the doldrum re- 
gion. The peak of this convective activity occurs in Colombia, in the northern 
hemisphere’s winter, spring or late fall seasons, when the solar energy in the 
vicinity of the equator is greatest. The seasonal rainfall patterns result from 
an alternate north and south migration of the centers of convective action in 
unison with the causative solar influences. No mention is made of monsoon 
conditions, but undoubtedly the mountain barriers to the west would prevent 
the development of any strong moisture-bearing monsoon winds from the 


Pacific Ocean. 


a September, 1959, by Phillip Z. Kirpich and Carlos 8. Ospina. 
1 prof, of Hydr. Engrg., Mass. Inst. of Tech., Cambridge, Mass. 
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The method of determining the spillway design storm and flood is unusual 
and is dictated bya lack of data and by the meteorologic characteristics of the 
storm rainfall. The method has considerable justification but some of the de- 
tails of the analysis might be questioned. The attempt to evaluate the 1,000- 
yr frequency for the 5-day precipitation from the 18-yr record is mathe- 
matically unrealistic. The rainfall record is so short that the slope of the 
frequency curve is very uncertain and values taken from a long extrapolation 
may not necessarily be conservative. In spillway design storms the object is 
to be conservative and to resort to the factor of safety principle may be the 
only rational approach. The method of expressing the rainfall depths as a 
percentage of the annual rainfall could be adopted for all durations from 1 to 
5 days. Then comparisons could be made with corresponding percentages that 
have been experienced in similar climates in other parts of the world, where 
the rainfall records are much longer. If the climate of a river basin has been 
correctly classified, the writer believes that statistical analogies may be 
made even from one continent to another. 

The procedure of plotting the known 5-day rainfall totals against the max- 
imum daily flow is practical and expedient in view of the lack of records and 
the non uniform character of the rainfall. However, the position of the en- 
velope curve is subject to individual judgment. It is questionable whether a 
straight line on semi-log paper is proper. If the plot is transferred to an 
arithmetic scale, the envelope curve will be found to be concave upward in the 
general form of a parabola with a vertical axis. If the infiltration losses for 
the storms are essentially constant, the envelope curve certainly violates the 
unity-hydrograph theory which shows that hydrograph peaks vary about as the 
first power of the volume of rainfall excess. If natural storage, of either 
channel, valley or lake type, is present there will be a tendency for the en- 
velope curve to be concave downward on an arithmetic scale. On the basis of 
the above discussion, and assuming that 186 mm represents the correct fac- 
tor of safety for the rainfall it is difficult to justify a daily flow in excess of 
2,500 cum per sec. The computed peak discharge of 6,000 cu m per sec with 
a Myers rating of 4,550, appears to place the river in a climatic province in 
which it obviously is not. As one schooled in the United States method of de- 
termining the maximum possible flood for spillways for major dams, the 
writer would be the last one to deplore the use of factors of safety in spillway 
design, but, in this case, it appears that the planners of the project do not 
have faith in their own very exhaustive meteorologic and hydrologic appraisal 
of the region. 

The adopted plan for flood protection and drainage of the Aguablanca Proj- 
ect, benefits from the application of hydrologic and hydraulic ingeniuity as 
well as from the sound economic analysis. The intercepting canal on the west, 
which carries the three mountain streams around the protected area, is a good 
solution which is undoubtedly aided by favorable land slopes ina southerly di- 
rection. The use of an excavated pond to reduce pumping capacities is a good 
practice which has been applied to the drainage of leveed areas in the United 
States, and, in some cases, to reduce the size of storm sewers in ordinary 
urban drainage problems. ~ 
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FRICTION FACTORS IN CORRUGATED METAL PIPES# 


Discussion by M. H. Diskin 


M. H. DISKIN.!—This paper is a valuable contribution to the knowledge of 
pipe flow in providing accurate experimental dataon loss of head in corrugated 
pipes of large size. It is to be regretted that the experiments could not be ex- 
tended to cover a larger range of Reynolds numbers, so as to show the trend 
of the f vs. Re curves before and after the maximum point observed, and in 
particular to show if the friction coefficient becomes constant at higher 
Reynolds numbers. 

With reference to Fig. 17 of the original paper, it is interesting to note 
that, if the recommended values of friction coefficient (corresponding to the 
maximum observed values) are considered as representing the values for 
rough turbulent flow, the equivalent pipe roughness, based on these values, is 
nearly the same for the three pipes and agreés with that computed for the 
smaller pipes reported. Taking Nikuradse’s equation: 


1/Vf = 1.74 + 2 log (D/2e) 


for the definition of the equivalent roughness (€), the values computed for the 
various pipes reported are given in Table 1. 


TABLE 1 


+ 


a September, 1959, by Marvin J. Webster and Laurence R. Metcalf. 
1 Hydr. Lab., Technion, Haifa, Israel. 
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REVISED COMPUTATION OF A VELOCITY 


HEAD WEIGHTED VALUE 
Discussion by Manuel A. Benson 


MANUEL A. BENSON, ! M. ASCE.—The revised computational procedure 
»roposed by the authors is apparently motivated by the difficulty of handling 
he extremely large numbers encountered in computing the velocity-head co- 
sfficient a. This difficulty is easily resolved by a coding process. 

Table 1 shows the slide-rule computation of a for a cross-section which 
1as been subdivided into five parts on the basis of hydraulic radius and rough- 
ness. Cols. 2 and 3 show the areas (a) and conveyances (k) of the separate 
subsections. 

Lines are drawn arbitrarily through the a and k columns, in effect shifting 
the decimal point, so that k/a-values are, in general, small numbers greater 
than 1. 

Col. 4 shows the slide-rule computations of k3/a2, All the significant fig- 
ires in Cols. 2 and 3 are used for this computation. Cols. 5, 6, and 7 are 
rapid mental computations based on only the figures to the left of the lines in 
the a and k columns. These computations are made in order to locate the 
lecimal points in the slide-rule computations of Col. 4. The correctly po- 
sitioned figures are shown in Col. 8. Cols. 4 through 7 need not be recorded; 
chey are shown only to illustrate the mental processes involved in locating the 
Jecimal points. 

The values in Col. 8 are coded values and do not represent the actual val- 
1es of the cubed conveyance divided by the squared area. However, a, com- 
outed as 


> (k3/a2) 
K3/A2 


s dimensionless and is correct regardless of the coding used. 
If the geometric mean of the conveyances at both ends of a reach is used in 
1 


he discharge formula Q = Ky, S”, it is possible to develop a formula for Q 
which makes unnecessary a trial-and-error procedure such as the authors 


use. The formula is: 


f September, 1959, by J. M. Lara and K. B. Schroeder. 
~ Hydr. Engr., Research, U. S. Geol. Survey, Washington 25, D. C. 


= 
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in which Q is the discharge in cubic feet per second; K; and Kg are the tota 
conveyances in section1 (upstream) and section 2 (downstream) respectively- 


Ke = a8 A R2/3 


1, .9 is the length of reach, in feet, between sections 1 and 2; Ah denotes th 
water-surface drop, in feet, between sections 1 and 2; g is the constant o 
gravity, 32.2 fps2; @ and @o represent the velocity-head coefficients a 
sections 1 and 2; and Ej 9 is the coefficient representing the percentage o 


TABLE 1 
----- Approximate values of----- 
Area Conveyance ee 2 2 3 /a2 
Subsection value of k/a (k/a)“ (k/a)*k k®/a' 
a b 3/42 
k? /a 
(1) (2) (3) (4) (5) (6) (7) (8) 
a 
z 25,15 33!1,000 573 1 1 33° 57.3 
| I 
2 1,2163 29!8,000 1660 2 4 120 166.0 
| 
1 
3 3113 63,000 268 2 4 24 26.8 
i} 
| 
4 1,11 42 241 2,000 1085 2 4 96 108.5 
| 
5 1,6!12 16} 0,000 1570 1 1 16 —-:15.7 
| 
Total 6,8/45(A) 109! 4,000(K) 281 2 4 400 281.0 


recovery of velocity head. For expanding reaches, Ej 2 is commonly as. 
sumed as 0.5; for contracting reaches, Ej 2 is assumed as 1.0 or 0.9. £ 
reach is contracting if the bracketed quantity in the denominator is positive 
expanding if it is negative. 

This method, like that of the author’s, permits computation of velocitie: 
and discharges in each subdivision. It is necessary only to proportion th 
total discharge on the basis of the relative conveyance in each subdivision. 

Using the base data of the author’s sample computations (areas, convey 
ances, and length of reach), the writer took 17 min to compute the two value 
of a and a discharge of 10,700 cfs by means of the suggested formula for di 
rect computation of discharge. 

The U.S. Geological Survey has developed similar formulas for multiple 
section reaches covering rather complex conditions, such as alternate con 
traction and expansion or abrupt contractions at a bridge opening. Thes: 
formulas provide direct solutions and eliminate trial-and-error procedures. 
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REVISED COMPUTATION OF A VELOCITY HEAD WEIGHTED VALUE? 
Discussion by Marcel Bitoun 


MARCEL BITOUN,! M. ASCE.—Just to keep the record straight, Ber- 
nouilli’s theorem, in the Appendix to Messrs. Lara & Schroeder’s paper, is 
not correctly applied. The eddy loss hy is always a subtractive term: 


hy =m | hy, - hy, | 


where the signs || indicate that only the absolute value of the difference 
ayy - hy, must be considered, irrespective of its algebraic sign. 


Then, for a contracting reach, and if m = 0.10: 
hg = f + (by, - hyy) - 0.10 | hy, - hy, | 
or, since hy» = hy: 


hg =! f.+4 (hy, — hy») ah 0.10 (hy, as hy, ) 


=f + 1.10 (hy, - hy,). 


a September, 1959, by J. M. Lara and K. B. Schroeder. 


_ 1 Chf. Design Branch, Div. of Flood Control, Pa. Dept. of Forests & Waters, Harris- 
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DISCHARGE FORMULA FOR STRAIGHT ALLUVIAL CHANNELS?@ 


Discussion by T. Blench 


4 i BLENCH,! F. ASCE.—The writer is interested in this paper as an ex- 
ample of the extreme swing away from formulas based mainly on theoretical 
speculations and idealizations, to empirical fitting curves that, apart from 
avoiding dimensional absurdity, seem to be unrelated to any dynamical con- 
siderations; however, they do seem to be biased in favor of the belief that 
Manning’s formula is of the proper type, for most alluvial channels with 
dunes. A neat fit has been obtained to a miscellany of slightly incongruous 
data by the patient selection of sufficient indices of non-dimensional groups of 
variables; the various graphs, combined with the references to sources of 
data, should be most useful to laboratory research workers. 

However, for various reasons, the writer would not use the various curves 
in practice. First, the data aremainly from laboratory flumes in which V b/v 
lies between about 105 and 106 and the breadth-to-depth ratio, b/d, is prob- 
ably round 3.0 or less; in the field V b/v usually lies between 106 and 108 
with b/d from about 5 to 30. Some of the data (Gilbert, Figs. 3,4) are for un- 
naturally uniformized bed material; some are for bed-conditions that had not 
become steady (Ref. 29); some are for mixed sediment with indeterminate but 
probably significant amounts in suspension (Refs. 23,29). The writer has al- 
ready pointed out2 that extrapolation of such data leads to results differing 
greatly from field experience in certain respects. Second, the curves and in- 
dices seem to have been adjusted to favor Manning’s formula whereas3;4 there 
isa large body of field evidence that Manning’s index of 2/3 should be re- 
placed by 3/4 to overcome the difficulty that his n is a function of hydraulic 
radius; the writer has experience of a canal system where, for apparently 
identical bed, n had to be specified at values from 0.018 for the main channels 
(10,000 cfs) up to 0.0225 for the small laterals (10 cfs)—these values, based 
on experience, correspond very closely to the indicial adjustment just men- 
tioned. Third, no cognizance is taken of bed-load charge which is important 
in the Gilbert experiments, in many river problems, and in occasional special 
canal ones; for example, one test problem’ showed that n altered from about 
0.021 down to 0.017 for a given channel breadth and 0.25 mm sand as a con- 
sequence of changing bed-load charge from about zero to 0.04% by weight, 
and increasing the slope to permit carriage of the load. (Note that Manning’s 
n may keep fairly constant for various discharges at one site in one channel 
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because the effects of varying charge and side shear at varying stages coun- 
teract each other to some extent and the incorrect index, 2/3, acts to reduce 
the variation of n; these are not the conditions of the problem just mentioned.) 
Fourth, there are simple and dynamically satisfactory formulas for mobile- 
bed channels,* based on an enormous collection of field data, that solve all 
the practical problems of the present paper, plus others involving bed-load 
charge. 

Although the so-called regime formulas? have the practical advantages of 
simplicity and as good accuracy as any in the field—since they were dis- 
covered by analyzing major field data—they do also have the fundamental ad- 
vantage of being consistent with the friction-factor diagram of rigid-boundary 
hydraulics, which they generalize. In fact, the regime flow formula can be 
written for study (it is written differently and very simply for design pur- 
poses) as: 


which is a generalized form of Blasius’ equation, as it should be since the 
boundary is formed from the water-sediment complex just as the so-called 
“smooth” rigid pipe boundary is formed from the pure water. (C is bed-load 
charge and a is a constant.) Moreover, when we consider that a bed-factor 
can be defined as Fp = v2/d and a side-factor as F,.= v3/b, the above equa- 
tion can be converted into: 


L 
Me wes Bo Oik ad Pall onthe it gen ee (2) 
in which x measures an equivalent roughness height and is equal to 


1 

(vF,)?/[Fp ad 2. C)?}. So now the formula is of “rough boundary” form as it 
should be since the bed is covered with dunes. In fact, regime channels are 
rough and smooth at the same time, drawing attention to the artificiality of 
these adjectives, as used in hydraulics, and to the fact that there should be a 
universal flow formula which breaks into special forms for special boundary 
phases merely by expressing the equivalent roughness height properly; this is 
consistent with the belief in a universal velocity distribution, independent of 
boundary phase, as held by authorities on rigid-boundary hydraulics. Revert- 
ing to rigid-boundary speculations, one of them converts Eq. 2 into the form 
of the Blasius equation by replacing x by laminar film thickness.° If a uni- 
versal flow formula exists, and if Eq. 2 is correct, then Manning for rigid 
boundary conditions must be amended to have 3/4 as an index instead of 2/3. 

Statistically the writer would like to see confidence bands on all charts to 
draw attention to the amount of variation that might be reasonably expected. 
Although this necessary statistical refinement is essential for a proper under- 
standing—its omission is comparable with describing the rainfall of a place 
by its long-term average without any indication of the scatter about the aver- 
age—it is seldom used; a few hydrologists have begun to expound its value in 
order to prevent the usual misinterpretation of the significance of flood fre- 
quency curves. In theabsence of confidence bands the writer can only give his 
opinion, based on experience of plotting hydraulic data, that quite a different 
looking set of indices would produce just as good fitting curves and, in par- 


ticular, so would a set based on the belief that Manning’s 2/3 ought to be re- 
placed by 3/4. 
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